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FOREWORD 

The  study  reported  herein  was  sponsored  by  the  U.  S.  Array  Elec- 
tronics Proving  Ground  (AEPG)  of  the  U.  S.  Army  Test  and  Evaluation  Com- 
mand under  Project  9-CO-002-000-010,  "Sensor  Surveillance  System."  It 
was  conducted  by  personnel  of  the  U.  S.  Army  Engineer  Waterways  Experi- 
ment Station  (WES)  during  the  period  August  1971  through  December  1972 
under  the  direct  supervision  of  Mr.  Bob  0.  Berm,  Acting  Chief,  Environ- 
mental Research  Branch  (ERB) , Mobility  and  Environmental  Systems  Labo- 
ratory (MESL),  and  under  the  general  supervision  of'  Messrs.  W.  G. 

Shockley,  Chief,  MESL,  and  W,  E.  Grabau,  Chief,  Environmental  Systems 
Division.  Field  tests  were  performed  in  August  1971  (wet  season)  and 
March  1972  (dry  season).  Personnel  making  significant  contributions  were 
Messrs,  H.  W.  West  (Project  Engineer),  of  the  Environmental  Simulation 
Branch,  M.  A.  Zappi,  Environmental  Characterization  Branch,  and  C.  Miller, 
B.  T.  Helmuth,  and  C.  Lebron,  SP4  T.  Engdahl,  and  SPi(  B.  Shumard  of  ERB; 
Messrs.  J.  R.  Curro,  R.  E.  Leach,  M.  Carlson,  D.  H.  Douglas,  E.  Perry, 
and  L.  0.  Steen,  and  Dr.  B,  Rohani  of  the  Soil  Dynamics  Division,  Soils 
and  Pavements  Laboratory  (SPL);  Mr.  J.  D.  Broughton,  Geology  Branch, 

SPL;  and  Messrs.  M.  Sava.ge  and  S.  W.  Guy  of  the  Operations  Branch,  In- 
st rument at icn  Services  Division. 

The  reduction  of  seismic  and  environmental  data  and  mappiing  of  the 
seismic  intrusion  detector  test  areas  were  accomplished  by  Mr.  Curro. 

The  theoretical  depth-ol’-penetra.tion  predictions  of  the  air-delivered 
sensors  were  produced  by  Dr.  Rohani.  The  report  was  prepared  by 
Mr.  West  and  Dr.  Rohani. 

Acknowledgment  is  made  to  Messrs.  L.  Jay,  Project  Coordinator, 

Colin  Giorgi , and  Merle  Wittraeyer,  and  LT  Gary  Wills  of  the  AEPG,  Fort 
Huachuca,  Arizona,  for  their  many  contributions  to  the  field 
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Huachuca,  Arizona,  for  their  many  contributions  to  the  field,  experi- 
mental program,  including  the  establishment  of  meteorological  sta- 
tions within  the  Fort  Huachuca  East  Range,  acquiring  test  vehicles, 
obtaining  high-speed  photographs  for  the  determination  of  impact  veloc- 
ities of  the  ADSID/S  sensors,  and  providing  personnel  support  for  con- 
ducting special  seismic  tests;  and  to  SP4  W,  Richardson  and  SP4  K. 
Gilliland  of  the  AEPG  for  their  assistance  and  cooperation  in  obtaining 
environmental  data  and  conducting  seismic  response  tests.  Acknowledg- 
ment is  also  made  to  CPT  E.  R.  Bedard  and  SGT  Serf ass  of  the  Combat 
Surveillance  Electronic  Warfare  School,  Unattended  Ground  Sensor  Divi- 
sion, Fort  Huachuca,  for  their  assistance  and  cooperation  in  obtaining 
penetration  data  on  the  ADSID/S  sensors . 

Director  of  WES  during  this  study  and  preparation  of  this  report 
was  COL  Ernest  D.  Peixotto,  CE.  Technical  Director  was  Mr.  F.  R.  Brown 
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CONVERSION  FACTORS,  METRIC  TO  BRITISH  AND  BRITISH  TO 
METRIC  UNITS  OF  MEASUREMENT 


XI 


Units  of  measurement  used  in  this  report  can  be  converted  as  follows : 
Multiply By  To  Obtain 


Metric  to  British 


millimeters 

0.03947 

inches 

centimeters 

0.0328 

feet 

centimeters  per  second 

1.968 

feet  per  minute 

square  centimeters 

0.1550 

square  inches 

cubic  centimeters 

0.06102 

cubic  inches 

meters 

3.281 

f<  et 

meters  per  second 

2.237 

miles  per  hour 

kilometers 

0.621*4 

miles 

square  kilometers 

0.3861 

square  miles 

kilograms 

2.205 

pounds 

kilograms  per  square  centime  ter 

1 ).  O 

. r_e_  _) 

pounds  per  square  inch 

grams 

0.0022 

pounds 

grams  per  cubic  centimeter 

0.0361 

pounds  per  cubic  inch.. 

kilopascals  (kilonewtons  per 
square  meter) 

0.1*450 

pounds  per  square  inch 

Celsius  or  Kelvin  degrees 

9/5 

Fahrenheit  degrees* 

Briti 

.sh  to  Metric 

inches 

0.0254 

meters 

pounds  per  square  inch 

6.8948 

kilopascals  (=  kilonewt 
pier  square  meter) 

cubic  feet 

O.02O3 

cubic  meters 

miles  (U.  S.  statute) 

1609.344 

.meters 

feci 

0 . 3048 

meters 

* To  obtain  Fahrenheit  (F)  readings  from  Celsius  (c)  readings,  use  the 
following  equation:  F = 9/5(C)  + 32.  To  obtain  Fahrenheit  from 

Kelvin  (k),  use:  F = 9/5  (K  - 273.15)  + 32. 
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SUMMARY 


To  obtain  a quantitative  basis  for  experiment  design  decisions 
regarding  evaluation  of  seismic  intrusion  detector  (SID)  systems,  test 
areas  were  selected  (one  for  hand-emplaced  STD's  and  six  for  air- 
implanted  STD's)  and  special  seismic  and  environmental  field  experiments 
were  conducted  at  Port  Huachuea,  Arizona,  during  August  3973  and  March 
1972.  Data  obtained  were  assumed  to  be  representative  of  wet-  and  dry- 
terrain  conditions  at  the  test  sites.  This  report  discusses  the  selec- 
tion of  the  SID  test  areas,  their  wet-  and  dry-terrain  characteristics, 
and  the  procedures  for  collecting  environmental  data  and  conducting 
seismic  response  tests  (i.e.  man-walking,  drop-hammer , M151  wheeled  and 
M577  tracked  vehicles,  ambient  and  induced  noises).  A3  so  presented  are 
experimental  detection  distance  data  on  the  hand-emplaced  miniaturized 
SID  (MINIS ID),  and  experimental  data  on  the  depth  of  penetration  and 
angle  of  impact  of  the  air-delivered  seismic  intrusion  detector/short 
(ADSID/S). 

This  report  also  presents  a theoretical  system  for  modeling  the 
quantitative  effects  of  the  terrain  on  SID  detection  performance  and 
methods  by  which  SID  test  results  can  be  extrapolated  from  one  site  con- 
dition to  another.  This  system  has  been  implemented  on  the  U.  3.  Army 
Engineer  Waterways  Experiment  Station  computer  and  involves  (a)  acquisi- 
tion of  terrain  and  seismic  data  for  input  to  the  model,  (b)  exercising 
the  model  for  each  set  of  terrain  factor  data,  and  (c)  portraying  the 
performance  predictions  for  each  set  of  terrain  and  seismic  factc^  data. 

A system  of  equations  is  presented  for  determining  soil  penetra- 
tion and  deceleration  of  air -delivered  SID's.  The  system  employs  a 
series  of  computer  routines  that  treat  the  soil  as  an  elastic-plastic, 
strain-hardening,  and  compressible  medium.  The  inpvuts  for  the  pene- 
tration equations  are  the  dynamic  stress-strain  data  obtained  from 
laboratory  triaxial  compression  and  uniaxial  strain  tests. 

The  hand-emplacement  area  was  found  to  provide  near-ideal  terrain 
conditions  for  performing  engineering  tesus  of  hand-emplaced  SID's. 
Surface  and  subsurface  terrain  lata  of  the  area,  show  more  variation 
from  site  to  site  than  between  wet-  and  dry-season  conditions.  The 
seismic  responses  measured  at  the  10  test  sites  within  the  hand- 
emplacement  area  show  some  site-to-site  variation,  and  this  variation 
should  be  considered  in  SID  engineering  testing;  however,  it  is  believed 
that  this  variation  can  be  accounted  for  analytically  and,  therefore, 


t 
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will  not  necessarily  result  in  any  erratic  or  inconsistent  results  in 
engineering  tests  performed  on  the  hand-emplaced  SID's.  The  variations 
in  seismic  response  parameters  portrayed  by  airdrop  areas  3 and  4-6, 
measured  in  the  dry  season,  generally  indicate  that  the  detection  per- 
formance of  air -implanted  SID's  can  be  evaluated  for  a large  range  of 
terrain  conditions . The  variation  in  these  conditions  allows  the 
testing  of  air-implanted  SID's  under  varying  stress  conditions  (i.e. 
high  to  low  rigid-body  deceleration  values)  and  under  a variety  of  soil 
strength  conditions . 

The  theoretical  detection  prediction  scheme  and/or  concepts  pro- 
vide an  analytical  framework  for  predicting  SID  performance  for  man- 
walking targets  in  various  world  terrains.  The  scheme  requires  specific 
terrain  and  seismic  input  data  for  execution,  and  as  a result,  four 
types  of  maps  were  prepared.  These  include  (a)  terrain  factor  complex 
maps  for  subsurface  layers  1 and  2 to  provide  inputs  to  the  theoretical 
seismic  wave  propagation  model,  (b)  surface  condition  maps  to  be  used 
with  the  forcing  function  model  under  development,  (c)  seismic  response 
factor  complex  maps  to  show  the  relative  seismic  response  of  the  hand- 
emplacement  area,  and  (d)  SID  one-man  walking  detection  performance 
map.  The  predicted  depths  of  penetration  of  the  ADS1D/S  compared 
reasonably  well  with  the  very  limited  experimental  penetration  data; 
they  show  the  depth  of  penetrations  to  be  very  dependent  upon  near- 
surface moisture  conditions. 

Appendix  A describes  the  operating  principles  of  hand-emplaced 
miniaturized  SID's  and  air-delivered  SID's  used  in  special  tests  to 
yield  data  for  aiding  in  the  evaluation  of  the  techniques  for  predicting 
SID  performance.  Appendix  B describes  the  test  sites  and  contains  sup- 
plementary refraction  seismic  data.  Appendix  C presents  the  equations 
for  predicting  penetration  depths  and  deceleration  values  for  air- 
implanted  SID's.  Appendix  D describes  the  computer  models  used  to  pre- 
dict SID  detection  performance. 
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EFFECTS  OF  TERRAIN  ON  THE  PROPAGATION  OF  MICROSETSMIC 
WAVES  AND  IMPLANTATION  CHARACTERISTICS  OF 
AIR-DELIVERED  SENSORS  AT  FORT  HIJACHUCA,  ARIZONA 
WET-  AND  DRY-SEASON  CONDITIONS 


PART  I:  INTRODUCTION 

Background 

1.  The  U.  S.  Army  Electronics  Proving  Gr5fend  (AEPG),  Fort 
Huachuca,  Arizona,  has  the  responsibility  for  developing  methodology  and 
conducting  engineering  tests  for  the  U.  S.  Army  Test  and  Evaluation  Com- 
mand (T.ECOM)  for  the  purpose  of  evaluating  seismic  intrusion  detection 
svstems . 

2.  The  development  of  a capability  for  testing  the  engineering 
adequacy  and  operational  suitability  of  seismic  intrusion  detection  de- 
vices must  necessarily  rest  on  an  unde rst and:. ng  of  the  physical  phenom- 
ena involved  in  the  detection  jjrocess,  Generally,  a seismic  intrusion 
detector  (SID)  is  intended  to  provide,  at  a minimum,  the  information 
that  something  (hereafter  referred  to  as  a target)  is  mov'ng  across  thc- 
terrain  within  some  specified  distance  from  the  SID.  It  would,  of 
course,  be  far  better  if  more  information  could  be  obtained,  and  thus 
efforts  ria-ve  been  made  to  design  SID's  that  will  provide  information  on 
target  characteristics  (i.e.  will  discriminate  among  types  of  targets), 
target  rate  of  movement,  target  direction,  and  target  range . 

3-  The  detection  process  alone  (i.e.  determination  that  there  is 
a target  within  a given  area)  is  relatively  complex,  and  with  the  addi- 
tion of  other  types  of  information,  the  process  becomes  very  complex 
indeed.  In  view  of  this,  it  seems  appropriate  to  provide  a brief  dis- 
cussion of  the  detection  process  and  of  the  physical  phenomena  involved 
therein . 


1 


Transfer  of  energy  from 
seismic  source  to  substrate 

It.  The  defection  process  consists  of  a sequence  of  interrelated 
physical  phenomena  (fig.  l).  The  sequence  starts  with  the  target . 


"W*8" 

SEISMIC  WAVE  PROPAGATION 


Fig.  1,  Signal  generation  and  propagation  of  microseismic  waves 


Anything  that  moves  on  the  surface  of  the  ground  applies  a force  to  the 
substrate.  The  thing  may  be  a man,  a vehicle,  the  pressure  pulse  of  an 
acoustic  wave,  an  animal,  or  even  trees  being  moved  slightly  by  the  wind. 
If  the  force  is  raised  to  a level  such  that  tne  substrate  is  stressed 
enough  to  deform,  even  if  only  minutely,  the  energy  is  carried  away  from 
the  point  of  deformation  by  seismic  waves.  The  implications  of  the 
statement  that  anything  that  applies  a force  to  the  substrate  of  suffi- 
cient magnitude  to  stress  and  deform  it  will  produce  a seismic  wave 
train  should  be  thoroughly  understood.  For  example,  wind  blowing 
through  trees  and  shrubs  transmits  forces  down  the  stems  and  into  the 
substrate,  and  the  deformations  result  in  the  generation  of  seismic 
wave  trains.  Those  waves  are  not  qualitatively  different  from  the  ones 
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produced  by  a moving  tank.  A passing  train  produces  seismic  signals; 
so  does  traffic  on  highways.  Factories  produce  seismic  signals . In 
fact,  the  active  crust  of  the  earth  produces  seismic  wave  trains 
throughout  the  surface  of  the  globe.  No  area  on  the  planet  is  entirely 
free  of  such  extraneous  wave  trains,  and  those  waves  will  inevitably  be 
sensed  by  military  SID's,  just  as  will  the  waves  from  a military  target. 
The  point  is  that  if  one  is  interested  in  detecting  a.  specific  type  of 
target  "false  alarms"  will  always  occur;  the  only  point  of  issue  is 
their  frequency. 

5.  The  amplitude  and,  to  some  degree,  the  type  of  seismic  wave 

generated  depends  upon  its  stress  history.  For  example,  a walking  man 
applies  a .localized  force  at  the  point  of  each  footfall,  and  that  force 
and  the  resultant  substrate  stress  are  characterized  by  a particular 
relationship  between  magnitude  and  time.  That  is,  the  force  is  not  ap- 
plied instantaneously;  rather,  it  rises  and  falls  according  to  a partic- 
ular pattern,  even  though  the  entire  time  involved  is  only  a fraction  of 
a second.  This  "force-time  history"  can  vary  remarkably.  The  major 
variations  can  be  classed  as:  (a)  variations  in  the  nature  of  the  tar- 

get itself,  and  (b)  variations  in  the  nature  of  the  substrate  (soil, 
rock,  pavement,  snow,  or  any  other  material  composing  the  medium  on 
which  the  target  is  moving). 

6.  Variations  in  the  nature  of  the  target  are  important  because 
they  result  in  differing  rates  and  magnitudes  of  force  application  to 
the  substrate.  That  is,  they  affect  the  force-time  history  by  which 
the  available  energy  is  applied  to  the  substrate.  The  major  potential 
variations  are  size  of  the  mass  and  rate  of  force  application. 

7.  Other  things  being  equal,  a large  mass  will  result  in  greater 
forces  than  a small  mass,  and  these  will,  of  course,  be  reflected  in 
the  force-time  history.  Thus , a large  man  tends  to  produce  a force-time 
history  different  from  that  of  a small  man;  this  also  is  true  for  large 
and  small  vehicles. 

8.  Again  other  things  being  equal,  a force  applied  rapidly  will 
produce  a force-time  history  different  from  that  of  the  same  force 
applied  slowly.  Thus,  a man  placing  his  feet  slowly  and  carefully  on 


the  ground  will  produce  a force-time  history  significantly  different 
from  that  produced  by  the  same  man  walking  normally.  Similar  differ- 
ences occur  in  the  force-time  histories  of  vehicles.  For  example,  a 
wheeled  vehicle  moving  on  a smooth , even  surface  applies  force  rela- 
tively smoothly,  without  sharp  discontinuities  in  force  level.  On  the 
other  hand,  the  track  pads  of  a tracked  vehicle  strike  the  ground  se- 
quentially, and  thus  produce  a force-time  history  consisting  of  a series 
of  high  and  low  force  levels.  Obviously,  the  rate  at  which  those  cycles 
occur  is  a function  of  vehicle  speed.  This  situation  becomes  even  more 
complex  if  the  surface  over  which  the  vehicle  is  moving  is  irregu'ar, 
since  in  this  instance  the  dynamic  response  of  the  vehicle,  regare j ess 
of  whether  it  is  wheeled  or  tracked,  will  result  in  maxima  and  minima  in 
the  force  applied  to  the  substrate. 

9*  Variations  in  the  nature  of  the  substrate  may  also  affect  the 
force-time  history,  since  that  history  is  a record  of  the  ielation  be- 
tween force  and  time  as  measured  at  a specific  point.  If  the  surface  is 
soft  and  spongy,  the  thing  applying  the  force  (i.e.  an  infantryman 'c 
foot,  the  track  pad  of  a tank,  etc.)  requires  a longer  time  between  ini- 
tial contact  and  achievement  of  maximum  force  level  than  would  be  the 
case  on  a firm,  hard  surface.  Thus,  the  force-time  histories  of  the  two 
situations  will  not  be  identical,  even  though  no  intrinsic  changes  in 
the  character  of  the  target  are  present.  One  implication  of  this  is 
that  the  force-time  history  of  a man  walking  in  a sod-covered  area  will 
be  different  from  that  of  the  same  man  walking  at  the  same  pace  in  an 
area  of  bare  ground. 

10.  The  force  applied  to  the  ground  surface  creates  a set  of 
stresses  in  the  substrate  material,  and  if  the  stress  level  is  high 
enough,  measurable  deformation  occurs.  The  energy  of  deformation  is 
then  carried  away  from  the  point  of  deformation  by  seismic  waves.  Thus 
it  follows  that  the  stress-time  history  and  the  characteristics  of  the 
substrate  materials  control  the  nature  and  magnitude  of  the  seismic 
waves . 

11.  There  are  three  major  seismic  wave  inodes,  all  of  which  move 
radially  outward  from  a point  of  substrate  deformation  (i.e.  the  point 
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at  which  stress  is  applied):  (a)  compression  waves  (P-waves),  in  which 

the  principal  particle  motion  is  along  a radial;  (b)  shear  waves  (S- 
waves),  in  which  the  principal  particle  motion  is  at  right  angles  to  a 
radial;  and  (c)  surface  waves  (Rayleigh  waves),  in  which  the  principal 
particle  motion  is  elliptical,  with  the  plane  of  motion  chiefly  perpen- 
dicular to  the  surface  of  the  propagating  medium  and  the  major  axis  of 
the  ellipse  usually  subparallel  to  the  surface.  Compression  waves  move 
in  all  directions , so  that  the  advancing  wave  front  is  approximately  a 
hemisphere.  Shear  waves  also  move  in  all  directions,  so  that  the  ad- 
vancing wave  front  is  also  approximately  a hemisphere.  However,  they  do 
not  necessarily  move  at  the  same  velocity  as  the  P-wave , even  in  the 
same  material.  Unlike  P-waves  and  S-waves , Rayleigh  waves  move  only  at 
shallow  depths  and  parallel  to  the  surface , so  that  at  a point  a few 
meters  from  the  target,  the  advancing  wave  front  is  cylindrical,  with 
the  vertical  axis  of  the  cylinder  at  the  target.  For  all  practical  pur- 
poses, Rayleigh  waves  (assuming  very  low  amplitudes  such  as  are  of  in- 
terest in  SID  utilization  and  design)  affect  a substrate  zone  equal  to 
about  one-half  wavelength.  Amplitudes  of  Rayleigh  waves  are  at  a maxi- 
mum at  the  surface  and  decrease  progressively  with  depth.  The  propaga- 
tion velocities  are  not  the  same  as  those  of  either  P-waves  or  the 
S-waves . 

12,  Generally,  both  P-waves  and  S-waves  are  diffracted  in  the 
direction  of  greater  substrate  densities,  and  since  nearly  all  natural 
substrates  increase  in  density  with  depth,  these  waves  tend  to  diffract 
downward.  The  result  is  that  they  appear  to  attenuate  very  rapidly  at 
the  surface.  They  may,  however,  reflect  from  a subsurface  discontinuity 
and  reappear  at  the  surface  at  some  distance  from  the  target.  Thus 
there  may  be  zones  in  which  the  P-waves  and  S-vraves  cannot  be  detected 
at  the  surface.  Conversely,  the  Rayleigh  wave  propagates  only  along  the 
surface;  thus  there  can  be  no  discontinuities  in  the  area  over  which  the 
signal,  can  be  detected.  This,  and  the  fact  that  about  two-thirds  of  the 
energy  at  the  source  is  carried  away  in  the  Rayleigh  wave , is  the  prin- 
cipal reason  for  choosing  the  Rayleigh  wave  as  the  mode  upon  which  to 
base  SID  designs . 
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13.  All  seismic  waves  are  affected  by  the  nature  of  the  materials 
through  which  they  propagate.  However,  since  the  Rayleigh  wave  is  the 
principal,  mode  of  interest,  the  following  discussion  will  be  restricted 
to  matters  of  concern  to  Rayleigh  wave  propagation. 

lU.  The  "ideal"  situation  for  the  propagation  of  Rayleigh  waves 
would  consist  of  a completely  homogeneous  elastic  half-space.  Anything 
that  departs  from  these  ideal  conditions  serves  to  attenuate  the  wave 
form  more  rapidly  than  would  be  the  case  by  purely  geometric  attenuation 
(i.e.  the  attenuation  resulting  from  the  same  amount  of  energy  being  ap- 
plied over  a longer  wave  front  as  the  wave  moves  outward  away  from  the 
source).  There  are  several  basic  types  of  nonhomogeneities  commonly 
found  in  nature.  The  dramatic  effects  of  nonhomogeneities  on  wave  prop- 
agation stem  primarily  from  the  fact  that  each  type  of  material  tends  to 
act  as  a specific  medium  for  a specific  suite  of  wave  frequencies.  That 
is,  a specific  type  of  material  tends  to  propagate  certain  frequencies 
more  efficiently  than  others.  The  result  is  that  substrate  materials 
act  as  selective  filters.  For  example,  most  targets  generate  a seismic 
wave  train  (or  signal)  consisting  of  a complex  of  frequencies  (or  wave- 
lengths) ranging  from  very  low  (i.e.  very  long  waves ) to  very  high  (i.e, 
very  short  waves).  Generally,  the  signal  contains  a broad  spectrum  of 
frequencies  as  it  emerges  from  the  source.  However,  as  the  wave  train 
noves  radially  away  from  the  source,  two  things  happen  to  it.  First, 
since  each  frequency  tends  to  propagate  at  a slightly  different  speed, 
the  wave  train  tends  to  separate  into  sections,  each  having  a character- 
istic frequency.  This  effect  is  usually  not  obvious,  since  the  process 
rarely  has  time  to  produce  complete  frequency  separation  before  the  wave 
train  dies  completely  away. 

15.  Second,  some  frequencies  are  propagated  efficiently  over 
rather  long  distances  and  others  die  out  quickly.  The  effect  of  this 
phenomenon  is  to  filter  out  some  of  the  original  frequencies,  leaving 
a signal  characterized  only  by  those  frequencies  that  are  efficiently 
propagated.  In  practical  terms , the  implication  is  that , in  some  ter- 
rains and  at  long  detection  distances , all  targets  tend  to  be  character- 
ized by  wave  trains  exhibiting  the  same  frequencies.  It  is  obvious  that 
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in  such  situations,  virtually  all  targets  will  look  the  same  insofar  as 
the  frequency  composition  of  their  signatures  is  concerned. 

16.  There  are  several  basic  types  of  variations  in  substrate 
characteristics,  each  of  which  interacts  with  the  seismic  wave  trains 
in  specific  ways.  One  such  basic  type  is  stratification.  Nearly  all 
substrates  are  stratified  to  seme  degree,  and  many  are  divided  very 
sharply  into  distinct  layers.  For  example,  most  agricultural  areas  are 
characterized  by  soils  having  a minimum  of  three  more-or-less  distinct 
strata;  (a)  a disturbed  layer,  down  to  plow  depth,  consisting  of  rela- 
tively low-density  materials;  (b ) a layer  in  which  the  parent  material 
is  only  partially  modified  by  soil-forming  processes;  and  (c)  a layer 
consisting  of  the  original  parent  material,  which  is  often  rock.  Each 
such  layer  tends  to  be  characterized  bv  somewhat  different  propagation 
characteristics,  with  the  result  that  the  seismic  waves  try  to  propagate 
independently  through  each.  The  interference  that  occurs  at  the  inter- 
faces absorbs  energy,  and  the  wave  as  a whole  therefore  tends  to  atten- 
uate rapidly.  Many  other  conditions  of  stratification  also  occur;  the 
hand  of  man  is  not  required  to  produce  them. 

17.  It  should  be  noted  that  stratification  is  not  a constant  with 
respect  to  time.  The  most  common  cause  of  a change  in  stratification  is 
a change  in  moisture  content.  The  presence  of  water  in  a soil  affects 
the  overall  density,  the  intraparticle  adhesion,  and  in  some  cases  even 
particle  orientation,  all  of  which  change  the  elasticity  and/or  viscos- 
ity of  the  soil,  and  therefore  the  wave  propagation  characteristics. 

The  effects  of  moisture  are  often  dramatic.  For  example,  the  deep, 
nearly  homogeneous,  silty  sand  soils  that  occur  along  the  Gulf  Coast  of 
the  United  States  seem  superficially  to  approach  the  ideal  as  a propaga- 
tion medium  and,  indeed,  do  approach  the  ideal  during  periods  when  mois- 
ture content  is  uniform  throughout  the  soil  mass.  However,  during  and 
for  a short  period  after  rain,  the  near-surface  soil  stratum  is  near 
saturation,  whereas  the  soil  at  depth  has  a much  lower  amount  of  con- 
tained water.  The  result  is  an  effective,  if  temporary,  stratification 
that  may  result  in  a more  rapid  attenuation  of  the  seismic  wave  than 
would  normally  be  the  case.  In  practical  terms,  this  means  that  the 


7 


effective  detection  range  of  a SID  might  be  significantly  reduced  during 
the  period  of  stratification.  Somewhat  similar  situations  can  occur  and 
persist  for  longer  periods;  for  example,  detection  distances  in  some 
tropical  situations  are  characteristically  less  during  the  wet  season 
chaii  during  the  dry. 

18.  A change  in  stratification  of  the  kind  described  above  may 
also  result  in  a change  in  the  frequency  characteristics  of  the  signal 
as  a function  of  time.  The  reason  for  this  is  that  a wet  soil  may  not 
efficiently  propagate  the  same  frequencies  as  the  same  soil  in  a dry 
state.  The  practical  effect  of  this  phenomenon  (see  paragraphs  l4  and 
15)  is  that  the  signals  arriving  from  the  same  target  at  the  same  dis- 
tance may  be  different  in  wet  and  dry  seasons. 

19.  A second  basic  type  of  variation  in  substrate  characteris- 
tics is  facies  changes.  Nearly  all  substrates  change  laterally  to  some 
extent , and  in  many  geographic  situations  such  changes  occur  abruptly 
and  within  short  distances.  This  phenomenon  may  have  important  effects 
on  seismic  wave  propagation  if  an  interface  occurs  between  the  target 
and  the  SID.  For  example,  consider  a situation  in  which  the  first  half 
of  the  propagation  medium  consists  of  a material  that  efficiently  propa- 
gates low-frequency  waves , but  quickly  attenuates  high-frequency  waves , 
and  the  second  half  consists  of  a material  that  efficiently  transmits 
high-frequency  waves,  but  not  low-frequency  waves.  The  signal  leaves 
the  target  with  a normal  spectral  composition  of  both  high  and  low  fre- 
quencies. During  the  fii'st  half  of  the  path,  all  of  the  high  frequen- 
cies are  attenuated  and  lost,  so  that  only  low  frequencies  remain  as  the 
signal  crosses  the  interface  between  material,  types.  Since  the  second 
type  will  efficiently  transmit  only  high  frequencies,  the  low  frequen- 
cies will  be  quickly  attenuated  and  lost,  leaving  nothing.  The  prac- 
tical implication  of  this  is  that  a seismic  wave  train  may  be  eliminated 
in  a distance  less  than  would  be  the  effective  detection  distance  if  the 
entire  propagation  path  were  in  either  one  of  the  substrate  types. 

20.  Nonhomogeneities  in  the  soil  mass  are  a third  basic  type  of 
variation  in  substrate  characteristics.  Many  substrate  materials  con- 
sist of  particles  that  are  both  small  with  respect  to  the  size  (i.e. 
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wavelength)  of  the  seismic  waves  that  propagate  through  the  substrate 
and  also  materials  with  particles  roughly  equal  to  the  size  of  the  prop- 
agating wave.  On  the  other  hand,  many  soils  consist  of  aggregates  of 
different  materials,  some  of  which  contain  relatively  large-size  parti- 
cles such  as  boulders  in  sand.  The  effect  is  somewhat  like  plums  in  a 
pudding.  In  addition  to  the  size  disparity,  there  are  many  situations 
in  which  the  "plums"  exhibit  elastic  and/or  viscous  properties  markedly 
different;  from  those  of  the  matrix  material. 

21.  Such  internal  nonhomogeneities  may  arise  from  a number  of 
causes.  For  example,  the  widespread  boulder-clays  of  the  northeastern 
United  States  are  materials  in  which  the  matrix  is  basically  silty  clay, 
and  the  plums  are  cobbles  and  boulders  of  solid  rock.  In  forested 
areas,  the  nonhomogeneity  may  be  caused  by  the  roots  of  trees.  In  arid 
and  subarid  climates,  discontinuous  masses  of  caliche  (a  form  of  calcium 
carbonate)  may  form  in  the  soil.  In  the  Arctic  the  plums  may  be  ice 
wedges  buried  in  the  soil. 

22.  Generally,  substrate  materials  exhibiting  this  type  of  non- 
homogeneity attenuate  all  seismic  waves  much  more  rapidly  than  a homoge- 
neous material  of  comparable  elastic  and  viscous  properties.  However, 
the  effect  is  usually  more  dramatic  on  high  frequencies  than  on  low. 

The  reasons  are  to  be  found  in  the  frequency-dependence  of  propagation 
efficiencies  (paragraph  l4).  For  example,  consider  a material  such  as  a 
fine  sand  with  large  rocks  in  which  the  matrix  efficiently  propagates 
low  frequencies.  A broad-spectrum  signal  propagating  through  such  a ma- 
terial will  be  rapidly  attenuated  because  of  a complex  of  phenomena: 

The  high-frequency  components  are  quickly  filtered  out  because  of  the 
properties  of  the  matrix,  the  low  frequencies  are  attenuated  because  the 
basic  wave  forms  are  distorted  as  they  pass  around  the  plums,  and  so  on. 

23.  Yet  another  basic  type  of  variation  occurs  as  changes  in  sur- 
face geometry.  Very  few  terrains  exhibit  perfectly  planar  surfaces, 
which  is  the  ideal  situation  for  the  propagation  of  Rayleigh  waves.  In- 
stead, nearly  all  are  irregular  to  some  degree.  Generally,  the  more  ir- 
regular the  topographic  surface,  the  less  efficient  is  the  propagation 
of  the  Rayleigh  wave.  There  is,  however,  one  important  proviso;  namely, 


9 


that  the  irregularities  he  large  enough  to  significantly  interfere  with 
the  wave  form,  but  not  so  large  that  one  or  more  wave  forms  can  he  ac- 
commodated on  the  feature.  Since  the  wavelengths  that  are  of  primary 
interest  from  the  point  of  view  of  SID  design  and  utilization  vary  in 
length  from  about  2 to  50  m,*  it  is  obvious  that  surface  geometry  fea- 
ture size  is  also  frequency-dependent.  That  is,  a feature  (such  as  a 
ditch)  that  is  a meter  across  may  well  interfere  significantly  with  the 
propagation  of  a high-frequency  wave  (i.e.  a wave  having  a length  of 
only  a few  meters);  whereas  it  might  not  interfere  significantly  with  a 
low-frequency  wave  (i.e.  a wave  having  a length  of  tens  of  meters).  The 
converse  is  also  possible  ; that  is , a high-frequency  wave  may  not  be 
significantly  affected  by  a feature  tens  of  meters  across,  because  sev- 
eral complete  wave  forms  of  the  high-frequency  wave  can  be  accommodated 
on  it . 

Transfer  of  seismic 
energy  from  substrate  to  SID 

24,  The  previous  paragraphs  have  dealt  with  the  transfer  of 
energy  from  a seismic  source  to  the  substrate  and  the  propagation  of 
seismic  waves  within  the  substrate.  The  following  paragraphs  present 

a brief  discussion  of  the  phenomena  dealing  with  the  transfer  of  seismic 
energy  from  the  vibrating  medium  (i.e.  substrate)  to  the  SID. 

25.  As  the  propagating  seismic  wave  passes  the  geophone  of  the 

SID,  the  geophone  is  carried  along  by  the  motion  of  the  substrate  par- 
ticles, if  the  geophone  is  properly  placed  in  contact  with  the  ground, 
and  if  certain  other  conditions  are  met.  In  rhe  case  cf  hand-emplaced 
SID's,  the  various  conditions  are  relatively  easy  to  meet,  but  this  is 
not  necessarily  the  case  with  air-delivered  or  artillery-delivered  de- 
vices. There  are  three  major  conditions  that  can  affect  SID  perfor- 
mance: (a)  depth  of  penetration  (or  emplacement),  (b)  geophone  atti- 

tude, and  (c)  coupling  of  the  geophone  with  the  substrate  material. 

2.6,  In  t.he  context  of  penetration  depth,  it  must  be  recalled  that 

* A table  of  factors  for  converting  metric  units  of  measurement  to 
British  units  and  British  units  of  measurement  to  metric  units  is 
presented  on  page  xi. 
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the  particle  motion  produced  by  a Rayleigh  wave  is  at  a maximum  at  the 
surface,  and  that  the  amplitude  decreases  as  a function  of  the  ratio 
depth/wavelength  of  the  substrate  material  and  becomes  effectively  un- 
detectable at  a ratio  of  3,  .2  in  almost  all  substrate  materials  3 The 
practical  effect  of  this  characteristic  of  the  wave  mode  is  that  maximum 
detection  distance  (or  range ) is  achieved  when  the  geophone  is  at  or 
very  close  to  the  surface.  Thus,  if  a geoplione  is  emplaced  at  a depth 
of  1 m,  it  may  not  detect  a target  that  is  readily  detected  by  a geo- 
phone placed  at  the  surface,  because  the  amplitude  of  the  particle  mo- 
tion at  that  depth  will  be  somewhat  lower.  This  effect  becomes  more  and 
more  pronounced  with  depth,  even  in  nonstratified  substrates. 

2".  If  the  substrate  is  strongly  stratified,  the  effect  may  be 
even  more  pronounced.  This  situation  can  occur  if  the  second  layer  has 
propagation  characteristics  markedly  different  from  those  of  the  surface 
layer.  In  this  situation,  it  is  common  to  find  that  very  little  of  the 
available  energy  is  coupled  into  the  second  layer,  and  therefore  the 
particle  motion  is  very  small.  Thus,  a geophone  emplaced  in  the  second 
layer  may  be  scarcely  affected  by  a wave  passing  above  it  and  which  is 
propagating  almost  entirely  in  the  surface  layer. 

28.  From  these  considerations,  it  is  evident  that  depth  of  pene- 
tration of  the  geophone  is  critical.  If  it  is  placed  too  deep,  the  geo- 
phone may  be  below  the  zone  in  which  active  particle  motion  is  taking 
place , 

29.  Ail  existing  SID’s,  and  nearly  all  those  presently  in  design 
stages,  utilize  geophones  that  operate  at  maximum  efficiency  when  the 
axis  of  the  geophone  is  parallel  to  the  pull  of  gravity.  The  point  is 
that  the  geophone  should  be  vertical  with  respect  to  sea  surface, 

not  vertical  with  respect  to  the  substrate  surface  (except  when  sea 
level  surface  and  .he  substrate  surfaces  are  parallel).  The  geophones 
become  less  efficient  with  increasing  angles  of  inclination,  and  most 
become  completely  inoperative  at  angles  from  the  vertical  of  between 
^5  and  60  deg,  depending  upon  type.  The  practical  effect  of  this  phe- 
nomenon is  that  two  SID's  placed  virtually  side  by  side  will  not  neces- 
sarily detect  the  same  target.  If  one  geophone  is  vertical  and  the 
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other  is  at  some  angle,  the  vertical  geophone  will  characteristically 
achieve  greater  detection  ranges  than  the  one  emplaced  at  an  angle.  It 
is  also  appax'ent  that  the  probable  detection  range  of  an  air-  or 
artillery-delivered  SID  cannot  be  reliably  estimated  unless  the  attitude 
of  the  SID  geopnone  is  known. 

30.  In  order  for  the  particle  motion  produced  by  the  wave  to  be 
sensed  by  the  geophone,  the  case  of  the  geophone  must  move  as  if  it  were 
a substrate  particle.  Generally,  this  can  only  be  achieved  if  the  geo- 
phone is  in  intimate  and  solid  contact  with  the  substrate . If  it  is 
not,  only  a fraction  of  the  real  particle  motion  may  be  transferred  to 
the  geophone , with  the  result  that  the  geophone  will  record  the  passing 
of  a wave  having  an  apparent  amplitude  of  much  less  than  that  of  the 
real  wave.  The  practical  effect  of  this  phenomenon  is,  again,  a reduc- 
tion in  the  detection  range, 

31.  Failure  to  achieve  adequate  coupling  can  arise  from  several 
mechanisms.  For  example,  the  impact  of  an  air-delivered  SID  may 
"crater"  the  substrate , leaving  only  the  very  tip  of  the  projectile  em- 
bedded in  the  ground.  This  limited  surface  contact  may  be  inadequate  to 
transfer  all  of  the  motion  from  the  substrate  to  the  SID,  The  impact  of 
the  SID  may,  in  some  soil  types,  cause  the  soil  to  surge  away  from  the 
projectile  body,  resulting  in  a tubular  cavity.  The  walls  of  the  cavity 
may  break  and  collapse  into  the  space  between  the  cavity  wall  and  the 
projectile,  which  results  in  a very  loose  material  actually  in  contact 
with  the  SID.  The  loose  material  may  act  as  an  energy-absorbing  system, 
effectively  preventing  any  significant  amount  of  seismic  motion  from 
penetrating  to  the  SID.  The  practical  effect,  of  course,  is  to  reduce 
the  effective  detection  range  of  the  device. 

32.  Finally,  the  passing  seismic  wave  produces  appropriate  motion 
in  the  geophone,  and  the  geophone  then  sends  an  electrical  analog  signal 
to  the  SID  logic,  which  the  logic  then  processes  in  some  way.  The  pre- 
cise operation  of  the  logic  is  of  major  importance  in  any  attempt  to 
predict  detection  range.  One  example  may  suffice. 

33.  Let  it  be  assumed  that  a particular  SID  has  a logic  system 
that  will  transmit  a "detection”  alarm  only  when  a particular  wave 
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amplitude  is  exceeded  four  times  within  a 5-sec  period.  Thus,  in  ef- 
fect, the  logic  "listens"  for  a particular  wave  amplitude.  When  it 
"hears"  it,  it  stores  the  memory  away  for  future  reference,  and  con- 
tinues to  store  away  similar  memories  for  a period  of  5 sec.  It  then 
adds  the  memories,  and  if  the  total  is  four  or  more,  it  transmits  an 
alarm. 

34.  Now  let  it  be  assumed  that  the  SID  is  in  a place  where  there 
is  a cyclic  background  "noise"  (i.e.  a wave  train  generated  by  some  sit- 
uation in  the  ambient  environment),  which  occurs  in  such  a way  that  the 
threshold  amplitude  of  the  SID  is  exceeded  once  every  2 sec.  The  SID, 
of  course,  hears  the  background  noise  signals,  but  since  only  two  such 
events  occur  in  each  5-sec  period,  the  SID  does  not  alarm.  Now  let  it 
be  assumed  that  an  intruder  approaches  the  position  of  the  SID.  At  some 
point,  while  the  intruder  is  still  some  distance  away,  the  SID  will  be- 
gin to  hear  some  of  the  intrude r's  footsteps.  It  will  not  hear  them 
all,  at  first,  because  of  natural  variations  in  the  intruder's  pace,  be- 
cause of  slight  changes  in  the  character  of  the  surface,  and  so  on. 
However,  the  SID  adds  the  record  of  those  footsteps  that  it  does  hear  to 
the  background  noise  events,  and  when  the  total  exceeds  four  in  a 5-sec 
period,  the  SID  alarms.  In  this  situation,  the  practical  effect  of  the 
background  was  to  increase  the  apparent  sensitivity  of  the  SID,  since  it 
achieved  a detection  at  a distance  that  it  would  not  have  achieved  had 
there  been  no  background  seismic  events. 

35.  Now  let  it  be  assumed  that  the  same  SID  is  emplaced  in  a re- 
gion in  which  the  background  noise  exceeds  the  threshold  amplitude  on  a 
random  schedule,  rather  than  on  a regular  cyclic  schedule.  Given  this 
situation,  a number  of  things  may  occur.  First,  a chance  concentration 
of  noise  events  may  cause  the  SID  to  alarm  in  the  complete  absence  of  an 
intruder.  A different  concentration  may  produce  the  same  effect  as  that 
described  in  paragraph  3^+  above,  thus  producing  an  abnormally  long  de- 
tection distance.  A still  different  random  concentration  of  noise 
events  might  also  saturate  the  SID  logic,  in  which  case  an  intruder 
could  pass  through  the  area  without  detection,  since  the  SID  would  be 
alarming  constantly;  the  signals  added  by  the  intruder  could  not  be 


differentiated  from  those  of  the  background  noise.  The  point  is  that 
the  effects  of  background  seismic  events  are  not  necessarily  either  con- 
sistent or  constant.  At  the  extremes,  they  can  produce  what  appears  to 
be  abnormal  sensitivity  in  a.  sensor,  and  they  can  prevent  detection 
entirely . 

Scope  of  engineering  tests 

36.  From  the  preceding  discussion,  it  is  apparent  that  any  objec- 
tive test  of  a SID  should,  if  possible,  account  quantitatively  for  all 
of  the  potential  variables  that  have  been  described.  Only  in  this  way 
is  it  possible  to  achieve  confidence  that  the  system  is  operating  as  in- 
tended, and  only  in  this  way  can  its  true  performance  characteristics  be 
established. 

Objectives 

37-  The  overall  objective  of  the  seismic  and  environmental  study 
reported  herein  was  to  improve  procedures  for  conducting  engineering 
tests  to  evaluate  various  SID  systems.  The  specific  objectives  were  to:. 

a.  Determine  the  significant  terrain  and/or  seismic  response 
factors  required  for  interpreting  the  engineering  test 
results  and  extrapolating  the  results  to  other  test 
conditions . 

b.  Determine  the  variation  (i.e.  site  to  site  and  season  to 
season)  in  the  significant  environmental  and  raicroseismic 
response  factors  measured  at  Fort  Huachuca, 

o_.  Develop  and/or  assemble  field  data  collection  procedures 
for  the  acquisition  of  the  significant  terrain  arid  seis- 
mic response  data. 

d.  Map  the  distribution  of  the  terrain  and  microseismic  re- 
sponse factors  of  the  hand-emplacement  test  area  consid- 
ered significant  to  the  engineering  evaluation  (i.e. 
those  factors  related  to  the  prediction  of  SID  detection 
range)  of  a typical  Phase  III  SID  system. 

e_.  Examine  WES  analytical  systems  to  determine  their  useful- 
ness in  extrapolating  SID  test  results. 
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f_.  Determine  theoretical  estimates  of  the  depth  of  penetra- 
tion and  deceleration  for  a typical  air-delivered  SID  in 
selected  airdrop  test  areas  under  varying  terrain  condi- 
tions (i.e.  under  two  different  moisture  conditions  that 
were  assumed  to  be  typical  of  wet-  and  dry-terrain  condi- 
tions , respectively ) . 


Approach 

38.  To  provide  a quantitative  basis  for  making  experiment  design 
decisions  in  regard  to  evaluating  SID  systems,  test  areas  were  selected 
(one  area  for  hand-emplaced  SID's  and  six  areas  for  air-implanted  SID's) 
and  special  seismic  and  environmental  field  experiments  were  conducted 
by  WES  personnel  at  Fort  Huachuca,  Arizona,  during  August  1971  and  March 
1972,  respectively.  Data  obtained  during  these  periods  were  assumed  to 
be  representative  of  the  respective  wet-  and  dry-terrain  conditions  at 
Fort  Huachuca.  Also , special  tests  with  air— implanted  SID  s were  con — 
ducted  at  Fort  Huachuca  by  Combat  Surveillance  Electronic  Warfare  School 
(CSEWS)  personnel  during  August  1972  to  provide  experimental  data  on  the 
depth  of  penetration  of  SID's  in  the  Fort  Huachuca  terrain. 

Sensor  detection  performance 

39.  To  provide  a framework  for  studying  the  quantitative  effects 
of  the  terrain  on  SID  detection  performance , the  detection  problem  was 
divided  into  efforts  to  develop  procedures  or  models  to  describe  the 
following  phenomena: 

a_.  The  transfer  of  energy  from  the  target  to  the  substrate 
material . 

b_.  The  propagation  of  the  microsei.smic  waves  within  the 
substrate  materials. 

c_.  The  transfer  of  energy  from  the  vibrating  medium  (sub- 
strate) to  a SID. 

d.  The  response  of  the  sensor  logic  to  the  vibrating  medium. 
These  four  facets  of  the  seismic  detection  problem,  are  related  as  shown 
in  fig.  2.  Each  facet  'was  studied  within  the  framework  of  theoretical 
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mechanics;  and,  where  possible,  mathematical  models  that  closely  simu- 
late the  phenomena  involved  in  each  facet  were  formulated.  The  input 
requirements  of  the  models  define  the  data  that  must  be  measured  cr  in- 
ferred to  make  a prediction. 

1)0 . The  performance  prediction  procedures  form  the  keystone  for 
extrapolating  SID  test  results  from  one  site  condition  to  another.  The 
concept  for  extrapolation  or  comparison  of  SID  performance  in  different 
areas  is  shown  graphically  in  fig.  3.  The  system  involves  (a)  acquisi- 
tion and  mapping  of  the  terrain  data  required  by  the  various  models , 

(b)  exercising  the  forcing  function  and  wave  propagation  submodels  for 
each  set  of  terrain  factor  data,  and  (c)  portraying  the  performance  pre- 
diction for  each  set  of  terrain  factor  data. 

Sensor  air-implantation  performance 

1)1.  The  development  of  quantitative  tools  for  predicting  the  per- 
formance of  air-implanted  sensors  has  been  restricted  to  studi.es  of  air- 
delivered  SID  (ADSID/S)  soil  penetration  and  deceleration.  The  approach 
to  the  soil  penetration  problem  employs  a series  of  computer  routines 
that  treat  the  soil  as  an  elastic-plastic,  strain-hardening,  compress- 
ible medium.  The  dynamic  stress-strain  data  required  as  input  to  the 
computer  routines  are  obtained  from  laboratory  triaxial  compression  and 
uniaxial  strain  tests. 

Program  structure 

1)2.  In  practice,  the  entire  research  effort  was  divided  into 
three  subprograms,  and  those  subprograms  are  used  as  the  basis  for  the 
structure  of  this  report.  Thus,  each  subprogram  is  described  in  a part, 
as  follows : 

a.  Part  II : Sensor  Test  Areas , Test  Procedures , and  Envi- 

ronmental Characterization.  This  includes  a description 
of  the  process  of  test  area  selection,  the  test  program, 
the  data  collection  and  processing  methods  employed,  and 
related  matters, 

b_.  Fart  III:  Analysis  of  Data.  This  includes  an  analysis 

of  both  environmental  and  seismic  response  data,  a dis- 
cussion of  seismic  propagation  models  and  their 
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utilization , a discussion  of  sensor  implantation  data  and 
models , and  related  matters . 

Part  IV : Mapping  of  Terrain  and  Seismic  Response  Factors 

and  SID  Performance  Values.  This  includes  a discussion 
of  techniques  developed  for  mapping  terrain  factors  re- 
lated to  seismic  sensor  utilization,  end  related  matters. 

Definitions 


Certain  terms  used  in  this  report  are  defined  as  follows : 
a..  Coefficient  of  friction  (p):  The  coefficient  of  friction 

between  the  projectile  (or  sensor)  and  the  soil, 
dimensionless , 

b_.  Compressibility  parameter  (pp/p):  The  ratio  of  the  cur- 

rent  wet  density  (p  ) of  the  substrate  material,  to  the 

P 

initial  wet  density  (p)  of  the  substrate  material.  Sym- 
bolically the  compressibility  parameter , pp/p  — exp  a , 


where 


is  the  vertical  or  axial  strain. 


Compression  wave  velocity  (Vp):  The  speed  of  a compres- 

sion (primary)  wave  through  a medium.  Compression  waves 
have  the  greatest  velocity  of  any  elastic  wave  in  the 
same  medium.  The  motion  of  the  particles  is  parallel  to 
the  direction  of  propagation,  V^  is  defined  mathemat- 
ically as 


- ("f-)1 


where 


Vp  = compression  wave  velocity,  LT 

O 

A = Lames  constant,  FL~ 

_2 

G - shear  modulus,  FL 

p = mass  density  or  soil  wet  density  divided  by 

-1*  2 

gravitational  constant,  g,  FL  T* 

Cone  penetration  resistance  (CPR):  The  force , in  kilo- 

pascals,  required  to  make  a 30-deg  right  circular  cone 


2 

Of  3.23  cm  base  (or  projected  area)  penetrate  a soil 
mass.  The  CPR  is  commonly  accepted,  as  a dimensionless 
index  of  the  shear  resistance  of  the  soil.  The  CPR  of 
the  surface  soil  and  the  average  CPR  of  the  0-  to  1+5-cm 
layer  are  used  in  this  report. 

Drag  coefficient  (C^):  The  coefficient  representing  the 

geometry  of  the  frontal  face  of  the  projectile  and  the 
frictional  forces  acting  on  the  projectile  during  a pene- 
tration event,  dimensionless. 

Moisture  (or  water)  content  (w):  The  ratio  of  the  weight 

of  water  (W  ) in  a sample  of  soil  to  the  dry  weight  of 

V 

solids  (W^)  in  the  sample,  expressed  as  a percentage.  It 
can  be  written 

W 

W = — 

W 

s 

Peak  particle  velocity  (PPV ) : The  maximum  positive  am- 

plitude of  the  time  rate  of  change  of  a particle  of  the 
medium  with  respect  to  a specified  reference  frame.  The 
particle  velocity  was  measured  in  the  time  domain  by  the 
geophones  used  in  this  study. 

Peak  summed  particle  velocity  (PSPV ) : The  peak  or  maxi- 

mum value  of  SPV. 

Rayleigh  wave  velocity  (VR):  The  speed  of  a Rayleigh 

wave (particle  motion  is  elliptically  retrograde  and  par- 
allel to  the  direction  of  propagation)  along  the  free 
surface  of  a medium;  depends  on  Poisson’s  ratio  (v)  of 
the  medium.  For  values  of  Poisson's  ratio  0 < v < 0.5  , 
the  Rayleigh  wave  velocity  has  the  range  0.875Vq  < V0 

< 0.955Vr,  , where  V = shear  wave  velocity,  LT“^  . 
b b 

Root  mean  square  (ms)  velocity:  The  velocity  corre- 

sponding to  the  square  root  of  the  sum  of  the  squares  of 
the  individual  velocities  (i.e.  the  particle  velocity  at 
each  frequency  component)  within  the  seismic  wave. 
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k.  Shear  wave  velocity  (Vg):  The  speed  of  a shear  wave 

(particle  motion  of  the  medium  is  perpendicular  to  the 
direction  of  propagation)  through  a medium.  Mathemat- 
ically, Vg  = (G/p)1/2 

1_.  Span ; Absolute  difference  between  the  maximum  positive 
and  maximum  negative  amplitudes  (or  velocities)  within 
the  seismic  wave. 

m.  Sum-frequency  ( f ) : The  frequency  at  which  the  maximum 

value  of  FSPV  occurs  on  the  SPV -versus- frequency 
relations . 

n.  Summed  particle  velocity  (SPV);  The  particle  velocity 
summed  in  the  frequency  domain  over  a 7-Hz  bandwidth. 

The  summation  is  accomplished  by  performing  a Fourier 
transform  to  the  analog  particle  velocity  recording  (time 
domain)  and  converting  it  to  the  frequency  domain. 

o.  Thickness  of  refraction  soil  layers  (T^  and  T,^)t  The 
thickness  of  layer  1 (T^)  is  the  vertical  depth  to  the 
interface  between  the  surface  and  the  next  deeper  layer 
as  distinguished  by  their  differing  compression  (primary) 
wave  velocities.  The  thickness  of  layer  2 (T,,)  is  the 
vertical  depth  from  the  bottom  of  the  first  layer,  as  de- 
termined above,  to  the  surface  of  the  next  layer.  The 
compression  wave  velocities  of  these  two  layers  are  de- 
termined by  techniques  of  refraction  seismology.  Tne 
above-defined  layers  often,  but  do  not  necessarily, 
correspond  to  soil  layers  as  defined  in  soil  moo  .anics 
studies.  Mathematically , the  thickness  of  layer  .1  is 


where 


T = thickness  cf  first  refraction  layer,  m 

X,  = distance  from  the  seismic  source  to  point  at 
x which  first  change  in  compression  wave  veloc- 
ity occurs , in 
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Vp  = compression  wave  velocity  in  first  layer , 
m/sec 

Vp2  = compression  wave  velocity  in  second  layer, 
m/  sec 


The  thickness  of  layer  2 is  determined  by: 


T .5*  P. 


V - V 
P3  P2 

V + V 
F3  P2 


where 


T2  = thickness  of  the  second  layer,  m 

X = distance  from  the  seismic  source  to  the  point 
^ at  which  second  change  in  compression  wave 


VP3 


velocity  occurs , m 

= compression  wave  velocity  in  third  layer , 


m/sec 

It  should  be  noted  that  a subsurface  layer  will  not  be 
detected  if  the  substrate  material  is  an  infinite , homo- 
geneous, isotropic,  elastic  medium. 

Ultimate .shear  strength  (c^):  The  maximum  or  ultimate 

shear  stress  of  the  substrate  material.  Mathematically, 

(“a  - °r) 


0 = 
u 


max. 


V? 


2 

cm 


where 


o = axiax  stress,  xg/um 
a 

a = radial  stress  (also  confining  pressure), 

r , , ? 

kg/cnr- 

Wet  density  (p):  Wet  ui.it  weight;  the  total  weight  (WT) 

of  the  wet  soil  (solids  plus  moisture)  from  a soil  sample 
per  unit  of  total  soil  volume  (VT)  of  the  sample  (volume 
of  solids  plus  voids).  Symbolically  this  is 


W„ 


P = 


cm 


Young's  modulus  of  elasticity  (E):  Young's  modulus  is 


the  stress-strain  ratio  in  simple  tension  or  compression. 
Mathematically 


E = 


de 

a 


where 

d = differential 
o = axial  stress,  kg/cm^ 

a = radial  stress  (also  confining  pressure), 
kg/ cm^ 

ea  = axial  strain,  percent 
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PART  II:  SENSOR  TEST  AREAS,  TEST  PROCEDURES,  AND 

ENVIRONMENTAL  CHARACTERIZATION 


1*1*.  This  part  of  the  report  discusses  the  selection  of  the  SID 
test  areas,  their  terrain  characteristics  and  the  procedures  used  in 
collecting  environmental  data  and  in  conducting  seismic  response  tests 
(i.e.  man-walking,  drop-hammer,  vehicle,  and  ambient  or  induced  noise), 
hand-emplaced  miniaturized  SID  (MINISID)  performance  tests,  and.  ADSID/S 
air-implantation  tests . 

1*5 . The  field  program  consisted  of  (a)  field  reconnaissance  and 
site  selection,  (b)  environmental  data  collection  to  provide  input  to 
the  sensor  detection  and  implantation  performance  models,  (c)  seismic 
response  tests  to  establish  the  seismic  response  characteristics  of  the 
various  test  sites  at  Fort  Huachuca,  and  (d)  data  collection  for  verifi- 
cation of  the  sensor  detection  and  implantation  performance  models. 

1*6.  The  data  collected  for  the  environmental  characterization  are 
listed  below. 

ia.  Input  requirements  for  detection  performance  model; 

(1)  Transfer  of  energy  from  source  to  substrate: 

(a)  Surface  soil  moisture  content,  percent 

(b)  Surface  soil  wet  density,  g/cm 

(c)  Surface  soil  type  (Unified.  Soil  Classification 
System,  USCS) 

(d)  Surface  soil  strength  (CPR),  kPa 

(e)  Surface  vegetation  ground  cover  height,  cm;  and 
density  of  ground  cover 

(2)  Propagation  of  seismic  waves: 

(a)  Thickness  of  refraction  layer,  cm  (each  layer) 

(b)  Compression  wave  velocity  m/sec  (ea.ch  layer) 

(c)  Shear  wave  velocity,  m/sec  (each  layer) 

(d)  Soil  wet  density,  g/cm  (each  layer) 

(e)  Surface  terrain  profiles,  xyz  coordinates 

b_.  Input  requirements  for  sensor  implantation  performance 
model : 
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(l)  Initial  Young's  modulus  of  elasticity  (E) , kg/cm 

2 

- (2)  Ultimate  shear  strength  (o^),  kg/cm 

(3)  Compressibility  parameter  (p^/p) , dimensionless 

(4 ) Soil  wet  density  (p),  g/cm^ 

(5)  Coefficient  of  friction  between  the  projectile  and 
the  soil  (p),  dimensionless 

(6)  Drag  coefficient  (Cp),  dimensionless 

Note:  These  data  are  determined  from  laboratory  tests  on 

bulk  samples  obtained  in  the  field. 
c_.  Background  meteorological  data: 

(1)  Air  temperature,  deg  Celsius 

(2)  Wind  speed,  m/s ec 

(3)  Wind  direction,  deg 

(4)  Precipitation  (rainfall,  snowfall),  cm 

U7.  The  tests  conducted  for  seismic  responses  are  listed  below. 

a. .  Basic  tests  for  model  development: 

(1)  Man  walking 

(2)  Drop  hammer,  or  controlled  source 

(3)  Vehicle 

(4)  Ambient  or  induced  noise 

b.  Tests  to  verify  prediction  models : 

(1)  Sensor  detection  performance  (MlNISID  tests) 

(2)  Sensor  implantation  performance  (ADSID/S  tests) 

Location  and  Climate  of  Fort  Huachuca  Reservation 


48.  The  Fort  Huachuca  Reservation  is  in  Cochise  County,  Arizona, 
about  76  km  north  of  the  Mexican  border.  The  reservation  is  on  the 
northeast  flank  of  the  Huachuca  Mountains  and  on  the  southern  part  of 
the  San  Pedro  River  drainage  basin.  The  climate  is  mild  and  sunny  and 
the  annual  precipitation  is  about  36  cm  (15-yr  record,  1956-1970)  of 
which  approximately  1.6  cm  is  the  water  equivalent  of  16  cm  of  snow. 

The  precipitation  distributed  by  months  is  as  follows,  (For  a detailed 
summary  of  the  meteorological  data  at  Fort  Huachuca,  see  reference  2.) 
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Month 

Precipitation 

cm 

Month 

Precipitation 

cm 

January 

1.52 

(0.17)* 

July 

11.30 

February 

1.U5 

CO 

0 

August 

9.22 

March 

1.1+5 

(0.4-1) 

September 

4.l4 

April 

0.1+8 

(0.07) 

October 

1.70 

(0.02) 

May 

0.18 

November 

1.24 

(0.20) 

June 

1.02 

December 

2.62 

(0.31) 

Total 

36.32 

(1.6) 

* Numbers  in  parentheses  indicate  snow . 

49-  The  mean  daily  maximum  temperature  varies  between  and  19  C 
(Celsius)  during  the  winter  months  (November-March ) and  between  23  and 
32  C during  the  summer  months  (April-October ) ; whereas  the  mean  daily 
minimum  temperature  varies  between  1 and  8 C during  the  winter  months 
and  between  11  and  19  C during  the  summer  months.  The  prevailing  -winds 
during  the  months  of  September-May  are  in  a southwesterly  direction 
(225  deg)  and  during  the  months  of  June-August  are  in  a westerly  direc- 
tion (270-285  deg).  The  mean  wind  speed  during  all  months  is  approxi- 
mately 3 to  4 m/s ec,  although  peak  wind  gusts  with  speeds  up  to  approxi- 
mately 19  to  26  m/sec  occur  during  all  times  of  the  year. 

50,  The  field  program  was  conducted  at  the  following  selected 
areas  on  the  East  Range  of  the  Fort  Huachuca  Reservation  and  at  the 
Willcox  Playa  near  Willcox,  Arizona : 

a..  A hond-empl anemone  area  (trapezoidal-shaped  with  an  area 
of  approximately  3 sq  km). 

b_.  Six  airdrop  (i.s.  air-implantation)  areas  (rectangular- 
shaped, each  100  m wide  and  1200  m long).  Five  areas 
were  located  on  the  Fort  Kuachuca  East  Range  and  one  at 
Willcox  Playa. 

£.  Special  airdrop  area  used  by  CSEWS  for  testing  ADSID/S. 
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Test  Areas 


51.  All  test  areas  were  in  Fort  Iiuachuca's  East  Range  (fig.  h), 
except  airdrop  area  6,  which  was  in  a flat  lake  bed  (Willcox  Playa), 
near  Willcox,  Arizona  (fig.  5).  The  East  Range  covers  approximately 
117  sq  km  and  lies,  in  general,  on  an  alluvial  apron  along  the  northeast 
flank  of  the  Huachuca  Mountains . Near  the  eastern  boundary  of  the 
range,  the  alluvial  apron  merges  into  the  floodplain  of  the  San  Pedro 
River.  The  East  Range  is  drained  by  intermittent  streams  (arroyos), 
which  flow  either  northerly  to  the  Babocomari  River  or  easterly  to  the 
San  Pedro  River.  Most  of  the  arroyos  are  well  incised,  contain  steep, 
near -vertical  banks,  and  vary  in  width  from  1 to  2 m to  as  much  as 

200  in. 

52.  The  genei'al  terrain  is  relatively  flat,  except  for  the  north- 
east corner,  which  is  flanked  by  hills  some  of  which  reach  heights  of 
100  m above  the  surrounding  terrain.  The  predominantly  flat  areas  vary 
between  el  1230  and  13o0.*  The  soils  are  predominantly  fine-to-coarse 
sands  (classified  SC  by  the  USC3)  with  a small  mixture  of  silt  and/or 
clay.  Gravels,  cobbles,  and  boulders  'up  to  30  cm  in  diameter  are  also 
quite  prevalent  within  the  substrate  and  usually  occur  within  23  cm  of 
the  surface.  The  upper  banks  of  the  arroyos  normally  contain  gravels 
and  large  stones  embedded  in  the  fine-co-coarse  sands.  Most  of  the 
soils  on  the  East  Range  are  strongly  calcareous , and  hardpans  of  nodular 
caliche  are  present.  Caliche  deposits  were  also  noted  in  a few  scattered 
areas,  but  were  normally  found  at  depths  greater  than  1 m below  the  sur- 
face. The  caliche,  in  most  cases,  seems  to  develop  as  a thin,  layered, 
local  deposit, 

33.  Vegetation  at  the  East  Range  includes  desert  grass,  mesqu.ite , 
creosote  bush,  acacia,  and  an  unidentified  bush  (probably  burro  bush). 
Stands  of  yucca  and  a few  cacti  were  also  present,  mostly  on  the  slopes 
of  the  arroyos,  but  were  not  very  abundant'  (i.e.  they  were  widely 


* All  elevations  (el)  cited  herein  are  in  meters  referred  to  mean  sea 
level . ... 
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spaced).  The  vegetation  in  general  was  most  dense  in  the  bottoms  and 
along  the  upper  banks  of  the  arroyos. 

Hand-emplacement  test  area 

5k,  The  trapezoidal-shaped  area  selected  for  con'  ting  te:  ts 
with  hand-emplaced  SID's  was  in  the  south-southwest  part  of  the  East 
Range  (see  fig.  4).  The  north  boundary  was  2880  m long,  the  east  bound- 
ary 20t0  m,  the  south  boundary  2970  m.  and  the  west  boundary  l68o  m. 

The  geographical  coordinates  of  the  four  corners  of  the  area  are  as 
follows : 


HE  corner: 
SE  corner: 
NW  corner: 
SW  corner: 


lat.  31°36’52"  N 
long.  110°15'50"  W 

lat.  31°35'53"  N 
long.  110  It '59"  W 

lat.  31°35't3"  N 
long.  110°17'17"  W 

lat.  31° 34' 55"  N 

“I  i ^ « O -»  Z'’  • « fM|  t t 

lung.  -L-lu  xo  ’ ,50  W 


55.  This  area  was  selected  as  being  representative  of  most  of  the 
physiographic,  geologic,  vegetation,  and  soils  conditions  found  on  the 
East  Range ; it  contained  numerous  homogeneous , relatively  flat  areas 
suitable  for  man-walking  tests.  Also,  its  network  of  unimproved  roads 
and  trails  provided  sites  for  vehicle  tests,  and  it  was  sufficiently  re- 
moved from  the  high-activity  (i.e.  noise)  areas.  The  east-west  orienta- 
tion of  its  long  axis  (see  fig.  4)  afforded  the  use  of  the  area  for 
testing  with  minimum  dissection  by  drainage  channels  (arroyos). 

Airdrop  test  areas 

56 . The  airdrop  areas  were  selected  to  provide  a variety  of  soils 
materials  ranging  in  penetration  difficulty  from  relatively  easy  to  ex- 
tremely difficult.  The  areas  were  tentatively  selected  on  the  basis  of 
available  published  data,  topographic  and  geologic  maps,  and  by  a ste- 
reoscopic examination  of  the  air  photos,  whereby  photo  patterns  were 
isolated  on  the  basis  of  their  tone,  texture,  and  geometry.  The  final 
selection  of  the  airdrop  areas  was  made  after  a detailed  ground  recon- 
naissance of  the  East  Range.  The  airdrop  area  at  Willcox  Playa  (see 
fig.  5)  was  near  the  center  of  the  large  playa  and  was  selected  to 


provide  an  easily  penetrable,  deep,  near-homogeneous  material  for  test- 
ing air-implanted  SID's. 

57.  Airdrop  areas  1-6  were  rectangular  in  shape  (each  100  m by 
1200  m)  and  were  positioned  at  the  following  geographical  locations  (co- 
ordinates  are  in  reference  to  the  SE  corner  of  each  area). 


a. 

Airdrop 

area 

1: 

lat.  31°38'13"  N 
long.  110°l6'40" 

b. 

Airdrop 

area 

2: 

lat.  31°39 ' 57"  N 
long.  110°l4 ' 12" 

c_. 

Airdrop 

area 

3: 

lat.  31°36' 36"  N 
long.  110°l8'15" 

d. 

Airdrop 

are  3 

4: 

lat.  31°37'57"  N 
long.  110°l4'10" 

e . 

Airdrop 

area 

5 = 

lat.  31°38'19"  N 
long.  U0°12'31" 

f_. 

Airdrop 

area 

6: 

lat.  32°09'47"  N 

long.  110°51'06" 

58.  A special  airdrop  area  selected  by  CSEWS  for  conducting 
drop  tests  with  an  ADSID/S  had  approximately  the  same  dimensions  as 
the  other  airdrop  areas.  It  was  located  at.  lat.  31°35'27"  N and 
long.  110°13'36"  W.  The  terrain  .description  of  the  special  airdrop  area 
and  the  results  of  the  experimental  tests  with  the  ADSIS/S  are  included 
herein  for  comparison  with  the  other  airdrop  areas  and  for  evaluating 
the  theoretical  depth-of-penetration  predictions  presented  herein. 

59-  Terrain  conditions  of  the  various  areas  are  discussed  below: 
a..  Airdrop  area  I , positioned  at  cl  4270  on  the  east  side  of 
and  parallel  to  an  unimproved  airstrip,  had  a general 
slope  of  approximately  3 deg  from  south  to  north  and  very 
few  irregularities  along  the  ground  surface.  The  sub- 
strate material  varied  from  a loose  sandy  silt  (classi- 
fied SC  by  the  USCS)  to  a reddish-brown  clayey  sand  (SC), 
except  for  the  northernmost  end  (last  20  m) . Here  the 
substrate  materials  consisted  of  a coarse  sand  with  cali- 
che nodules  and.  considerable  amounts  of  gravel,  cobbles, 
and  rocks  up  to  10  to  15  cm  in  diameter. 
b_.  Airdrop  area  2 was  situated  at  el  435C  on  the  west  slope 
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of  an  igneous  hill  near  the  northeast  corner  of  the  East 
Range  (see  fig.  *0 . This  area  was  selected  because  it 
represented  the  upper  range  of  soil  strengths  that  might 
be  encountered  on  the  East,  Range.  The  ground  surface 
along  the  area,  was  very  irregular  and  consisted  of  a 
coarse  gravelly  sand  with  some  cobbles  and  boulders  up 
to  50  cm  in  diameter.  Patches  of  bedrock  were  exposed 
within  the  area. 

Airdrop  area  3 was  a rather  large  arroyo  whose  top  width 
was  approximately  200  m and  bottom  width  was  approxi- 
mately 150  m.  The  ground  surface  along  the  bed  of  the 
arroyo  was  relatively  flat  and  was  composed  of  a fine-to- 
coarse  sand  with  some  mixture  of  silt  and  clay  (SC  to  SM 
by  USCS).  A small,  but  deeply  entrenched,  wash  channel 
(2  m wide  and  5-7  m deep)  was  also  in  the  main  bed  of 
the  arroyo.  This  wash  channel  contained  an  assortment  of 
sand,  gravel,  pebbles,  and  a few  small  be  J.ders.  The  bed 
of  the  arroyo  also  contained  a very  dense  stand  of  reeds 
(a  species  of  grass)  1 to  2 m in  height. 

Airdrop  area  4 was  on  an  alluvial  plain  of  the  San  Pedro 
River  at  el  i+150  (fig.  4).  The  ground  surface  along  the 
center-line  length  was  relatively  flat  (less  than  2-deg 
slope),  but  along  the  southwest-northwest  boundary  it  was 
very  irregular  as  a result  of  numerous  small,  well- 
defined  wash  channels.  The  substrate  material  consisted 
of  a loose,  gray,  sandy  silt  with  pebbles  in  the  surface 
layer  (0-10  cm)  and  a firm,  light-brown , sandy  clay  with 
various  amounts  of  friable  caliche  down  to  a depth  of  1 m. 
A' layer  of  caliche  approximately  30  cm  thick  was  encoun- 
tered at  a depth  of  approximately  1 m near  the  southwest 
corner  of  the  area.  The  dominant  vegetation  consisted 
of  creosote  bush  approximately  1.5  m in  height,  with  a 
few  very  widely  spaced  mesquite  trees  2-3  m in  height. 
Airdrop  area  5 was  also  in  an  alluvial  plai a of  the 
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San  Pedro  River  at  el  4100-1+150  (fig.  4).  The  topography 
was  slightly  undulating  as  a result  of  small  well-defined 
wash  channels  that  drain  the  area  in  the  east-northeast 
direction.  The  substrate  materials  in  general  consisted 
of  a sandy  silt  (SC)  in  the  surface  layer  (0-10  cm)  and  a 
very  hard,  compact,  sandy  silt  with  inclusions  of  pebbles 
and  caliche  in  the  subsurface  layer  (10-120  cm).  A layer 
of  white  caliche  was  also  detected  at  a depth  of  approxi- 
mately 100  cm  in  the  southeast  part  of  the  area.  Vege- 
tation covering  the  area  was  composed  of  a relatively 
dense  mixture  of  creosote  bush  (height  1-1.5  m)  and 
acacia  (height  1-1.5  m). 

.f . Airdrop  area  6 was  in  a rather  large  lake  bed  (Vfillcox 
Playa)  (see  fig.  5).  The  topography  was  extremely  flat 
and  very  smooth.  The  substrate  material  in  the  0-  to 
20-cm  layer  was  a compact,  brown,  silty  clay  (CL),  and  in 
the  20-  to  150-m  layer  it  was  a homogeneous , blue-green 
to  black  clay  (CH).  The  playa  contained  no  vegetation 
and  is  inundated  during  certain  periods  of  the  year 
(July-September ) . 

£L*  Special.  G'SEWS  airdrop  area  was  on  an  unimproved  airstrip 
located  approximately  2 km  west  of  the  hand-emplacement 
area  (see  fig.  4)  at  el  4300,  The  substrate  material 
varied  from  sandy  silt  to  clayey  sand  (SC)  with  small 
amounts  of  gravel  and  small  pebbles.  The  ground  surface 
was  relatively  smooth  and  contained  a small  grass  cover 
along  the  outer  edges  of  the  airstrip. 


Seismic  Response  Tests 


Test  site  selection 

60.  Within  each  sensor  test  area,  one  or  more  sites  were  selected 
for  conducting  seismic  response  tests  and  collecting  environmental  data. 
In  the  hand-emplacement  area,  10  sites  (sites  335-342,  348,  and  349) 
were  selected  for  conducting  man-walking  and  drop-hammer  seismic 
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response  tests,  and.  two  sites  (cross-country  sites  344  and  345)  were  se- 
lected for  conducting  vehicle  seismic  response  tests.  The  locations  of 
all  sites  are  shown  on  the  aerial  mosaic  in  fig.  6.  Within  each  airdrop 
test  area,  except  for  the  special  CSEWS  area,  one  site  was  selected  for 
conducting  man-walking  and  drop-hammer  seismic  response  tests.  All  man- 
walking  and  drop-hammer  test  sites  were,  in  most  cases,  in  areas  portray- 
ing the  least  ground  surface  irregularity  within  the  general  area. 

Test  site  layout 

6l.  Sites  were  prepared  in  a specified  manner  for  the  tests  (drop- 
hammer,  man-walking,  and  moving-vehicle  tests).  For  the  drop-hammer  and 
man-walking  tests , a reference  line  110  m long  was  established  by  plac- 
ing wood  stakes  at  5-m  intervals.  For  the  two  vehicle  cross-country 
test  sites,  a reference  line  1100  m long  was  established  by  placing  wood 
stakes  at  50-ra  intervals.  Typical  site  layouts  are  shown  in  fig.  7» 

Test  procedures 

The  seismic  responses,  or  analog  signatures  , resulting  from  a. 
walking  man  and  a controlled  source  (drop-hammer)  were  measured  at  each 
test  site.  Also,  the  ambient  background  noise  and  any  induced  noises 
that  occurred  during  testing  were  recorded.  Special  vehicle  tests  were 
also  conducted  in  the  hand-emplacement  area  with  two  military  vehicles 
(M151  Jeep  and  M57T  .Armored  Command  Post).  The  following  paragraphs  dis- 
cuss the  instrumentation  used  and  the  procedures  for  each  type  of  test. 

63.  Instrumentation . The  analog  signatures  resulting  from  a 
walking  man,  a hammer  drop,  or  a moving  vehicle  were  measured  with  the 
following  scientific  geophones:  a Mark  Products  Model  L-4-3D  (three- 

dimensional  geophone)  positioned  at  sta  0 (see  fig.  7)  and  a Mark 
Products  L-4-1D  (one-dimensional,  or  verticil.,  geophone)  positioned  at 
. sta  -5  for  the  man-walking  and  drop-hammer  tests , and.  an  1-4-3D  geophone 
positioned  at  sta  0 for  the  vehicle  tests.  The  natural  frequency  of  the 
geophones  was  1 Hz.  The  L-4-3D  geophone  was  emplaced  by  excavating  a 
hole  approximately  30  cm  in  diameter  and  22  cm  deep  and  orienting  the 
three  internal  geophones  as  follows : The  vertical  geophone  was  normal 

to  the  ground  surface,  the  radial  geophone  was  parallel  to  the  reference 
line,  and  the  transverse  geophone  was  perpendicular  to  the  line.  The 
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100  m 


1000  m 


4-  900  m 


400  m 


300  m 


200  m 


Triaxis  Geophone v 

(Mark  Products  L-4-3D h1 

Vertical.  Geophone 
(Mark  Products  L-4-XD) 


100  m 


Drop-hammer  and  man-walkin 


itc  layout  (di, stance  markers 


at  5~m  intervals ) 


1 -100  m 

b.  Vehicle  site  layout 


distance  markers  at 


50-m  intervals 


Fig.  7.  Site  layouts  for  seismic  response  tests 

material  removed  from  the  hole  was  then  tightly  packed  around  and  also 
spread  over  the  top  of  the  geophone.  The  same  emplacement  procedure  was 
used  for  the  L-4-1D  geophone  at  sta  -5,  except  that  a hole  approximately 
20  cm  in  dimeter  was  excavated.  Fig.  8 shows  two  geophones  in  place. 

64.  Seismic  analog  signatures  were  recorded  on  analog  M magnetic 
tape  and  on  a direct-print  oscillograph  (fig.  9),  which  were  housed  in 
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a mobile,  self-contained  vehicle  that  was  normally  positioned  approxi- 
mately 50-75  m from  the  site  (or  walk  line).  A vide-band- frequency 
(2-300  Hz)  battery-powered  recording  system  was  used  to  obtain 

f 

continuous  high-quality  data  and  an  adequate  dynamic  range.  All  instru- 
ments used  in  the  seismic  tests  were  calibrated  in  accordance  with  the 
National  Bureau  of  Standards  procedures,  and  no  filtering  of  any  type 
was  used.  Seven  channels  were  used  to  record  the  following: 

a.  Four  particle  velocity  measurements  (one  three- 

dimensional  geophone  array  at  sta  0 and  one  vertical 
geophone  at  sta  -5 ) . 
b_.  One  time  code. 

c_.  One  voice  comment  or  remote  acoustic  data  (microphone), 
d.  One  location  of  energy  source  (position  indicator). 

65 . Man-walking  tests.  In  a man-walking  test , a man  walked  at  a 
constant  rate  along  the  prescribed  110-m  path,  starting  from  sta  +100. 
During  the  test,  radio  communication  was  maintained  between  the  man  and 
the  recorder  for  positioning  each  footstep  (i.e.  each  foot  encounter 
with  the  ground)  along  the  line.  Prior  to  the  actual  test,  a practice 
walk  was  made  to  establish  a step  of  uniform  length  and  rate. 

66.  Drop-hammer  tests.  A special  drop-hammer  apparatus  (fig.  1C) , 
designed  by  U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES)  was  a 
controlled  energy  source  whose  seismic  response  could  be  compared  in  am- 
plitude and  frequency  with  that  resulting  from  an  average  footstep.  A 
sample  analog  signature  from  the  drop  hammer  is  compared  with  that  from 

a footstep  in  fig.  11.  The  hammer,  which  consists  of  a 9-9~kg  annular 
weight  with  a Teflon  sleeve,  is  connected  by  a metal  guidepost  to  the 
footing,  which  is  a plywood  base  topped  by  two  layers  of  0.2-cro-thick 
hard  rubber  and  two  layers  of  1 .9-cni-thick  soft  foam  rubber.  In  a drop- 
hammer  test  the  hammer  was  dropped  two  times  at  each  5-m  station  along 
the  reference  line  starting  at  sta  0 and  continuing  out  to  the  distance  at 
which  the  seismic  signal  could  not  be  detected  over  the  system  and  ambient 
noise  level.  Before  the  hammer  was  planed  on  the  soil  at  each  5-m  sta- 
tion, surface  litter  was  removed,  and  the  hammer  was  positioned  so  that 
it  would  drop  vertically  and  have  uniform  contact  with  the  ground. 
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a,  drop-hammer  signatures  for  five  drops 


NOTE:  THESE  RECORDS  WERE  OBTAINED  FROM  THE  VERTICAL 
GEOPHONE  FOR  THE  10-METER  DISTANCE. 


SCALE,  TIME 

O 10  *0  MSEC 

l_j ! i__l _ I 


Fig.  11.  Comparison  of  drop-hammer  and  footstep  signature 
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Vehicle  tests  ■were.-  conducted  in  the  hand- 


67  Vehicle  tests. 
emplacement  area  with  two  vehicles  (MIS 1 and  M57Y ) in  the  same  manner  as 
the  man-va.lkiug  tests,  A vehicle  was  positioned  at  some  arbitrary  dis- 
tance beyond  sta  +1000  so  that  it  could  attain  a preselected  speed  upon 
reaching  sta  +1000,  Speeds  of  8 and  32  km/hr  were  used,  and  two  tests 
were  conducted  at  each  speed. 

MINISID  and  ADSID/S  Tests 

MIN I SID  tests 

68.  Man-walking  tests  were  also  conducted  with  four  Phase  III 
MINISID 's  at  each  of  the  seismic  test  sites  to  yield  data  for  aiding 
in  the  evaluation  of  the  techniques  for  predicting  Sip  performance. 

(For  a discussion  of  MINIS ID  operating  principles  see  Appendix  A.) 

The  MJNISID  test  procedures  are  discussed  below. 

6,9.  A small  trench  approximately  5 cm  wide,  10  cm  long,  and  5 cm 
deep  was  excavat.ed  at  sta  0 (fig.  7)  for  emplacement  of  the  four  MINISID 
geophor.es.  Each  geophone  was  placed  in  the  trench  by  hand  and  forced, 
down  until  its  spike  was  completely  below  the  bottom  of  the  small  trench. 
After  the  geophone  was  in  place  ana  as  nearly  vertical  as  possible  (the 
MINISID  geophone  senses  the  vertical  component  of  the  seismic  wave), 
slight  pressure  was  applied  to  the  side  of  the  geophone  to  ensure  that 
it  did  not  move  and  was  in  firm  contact  with  the  soil.  Emplacement  of 
the  four  MINISID' s is  shown  in  fig.  12.  Once  the  geophones  were 


a.  Geophones  and  detectors  b.  Geop’none  in  trench 

Fig.  12.  Emplacement  of  MINISID  geophones 
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emplaced,  the  material  removed  from  the  trench  vas  tightly  packed  around 
and  loosely  spread  over  the  top  of  each  of  the  geophones. 

70.  The  triggering  of  the  MINISID's  was  monitored  by  using  a por- 
tatale  that  was  positioned  in  the  instrumentation  van  located  approxi- 
mately 50-75  m from  the  walk  line.  The  four  MINISID's  were  set  on  the 
following  gains . 

an  MINISID  A:  Low  gain. 

b_.  MINISID's  B and  C:  Medium  gain. 

c_.  MINISID  D:  High  gain. 

71.  Individual  walks  were  made  for  each  of  the  four  MINISID's. 

For  the  first  test  (i.e..  MINISID  A,  which  was  set  on  low  gain)  a man 
walked  from  sta  +100  down  the  110-m  line  until  an  activation  was  indi- 
cated on  the  portatale . At  the  point  of  activation  the  man  was  signaled, 
by  radio,  to  stop  the  walk  and  to  dete."mine  his  position  (i  e.  distance) 


with  reference  to  the  MINISID.  This  distance  was  rec 


..  the  detec  ■ 


tion  distance  of  the  MIN£pID.  The  man  then  walked  back  to  sta  +100,  or 
at  least  15  m back  fron/the  detection  distance,  paused  for  approximately 

/ 30-45  sec,  and  then  started  walk- 
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Fig.  13.  ADS ID /S  sensor 


ing  down  the  line  again.  Walk 
tests  were  conducted  until  the 
MINISID's  had  triggered  several 
times  for  each  gain  setting, 
respectively. 


72.  Airdrop  tests  were  con- 
ducted by  personnel  of  CSFWS  at 
Fort  Huachuca  during  August  1972 
with  a dummy  (i.e,  not  instru- 
mented) Phase  III  ADSID/C,  MA.-37 , 
at  the  special  airdrop  area  (unim- 
proved airstrip)  described  above. 
Fig.  13  is  a drawing  of  the 
ADSID/S.  The  results  of  these  ex- 
perimental tests  have  been 
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included  herein  to  aid  ir.  the  evaluation  of  the  techniques  for  predict- 
ing the  penetration  performance  (i.e.  depths  of  penetration)  of  the 
ADSID/S's.  These  tests  are  discussed  below. 

73*  There  were  120  ADSID/S's  delivered  from  a fixed-wing  aircraft 
by  dropping  three  censors  on  each  of  several  passes  over  a marked  target 
area.  The  drop  interval  for  the  three  sensors  was  regulated  by  an  in- 
tervalometer  mechanism  and  normally  resulted  in  a ground  impact  spacing 
of  100-150  m.  After  each  completed  aircraft  flyby,  a CSEWS  field  team 
assisted  by  one  man  from  WES  -visited  the  impact  area  to  determine  each 
sensor's  depth  of  penetration  (i.e.  penetration  depth  measured  along  the 
longitudinal  axis  of  the  impact  trajectory),  angle  of  impact  with  re- 
spect to  the  vertical,  and  data  on  the  terrain  conditions  (i.e.  moisture 
content,  density,  cone  penetration  resistance).  In  addition,  the;  AEPG 
provided  a team  to  take  high-speed  photographs  for  the  determination  of 
the  near-surface  impact  velocities  for  11  of  the  sensor  tests  (tests  5, 

./V.  , 

9,  11,  12,  3U,  36,  75,  76,  and  78-BC).  ' ■■ 


Ambient-  and  Induced-Noise  Tests 


74.  The  seismic  motion  resulting  from  both  the  ambient  or  natural 
noise  and  any  induced  noises  created  by  aircraft  (fixed-  and  rotary- 
wing), vehicles  traveling  along  roads  near  the  areas,  wind  gusts,  rain, 
etc.,  were  recorded  at  each  test  site.  Approximately  60  sec  of  dataware 
recorded  at  each  site  as  follows;  The  ambient  seismic  motion  was  mea- 
sured at  a time  when  members  of  the  test  team  were  quiet  and  motionless 
and  normally  before  the  start  of  a drop-hammer  test;  in  the  case  of  in- 
duced noises,  the  seismic  motions  were  measured  as  they  occurred  during 
testing. 

Environmental  Characterization 

75.  To  obtain  information  that  allows  the.  SID  test  results  to  be 
evaluated  within  a rigorous  theoretical  framework , specific  environmen- 
tal data  are  required.  For  convenience,  the  environmental  data 
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collected  can  be  considered  in  support  of:  (a)  the  seismic  response 

tests  and/or  the  prediction  of  sensor  detection  performance,  and  (t>) 
sensor  implantation  performance.  The  data  requirements  for  each  cate- 
gory are  listed  in  paragraphs  U6a_  and  b_. 

76.  In  general,  the  terrain  factors  required  (except  soil  type, 
Atterberg  limits,  and  soil  particle  sizes)  are  sensitive  in  varying  de- 
grees to  changes  in  moisture  content  and  temperature.  For  this  reason 
the  field  data  collection  program  was  conducted  under  conditions  assumed 
to  be  representative  of  wet-  and  dry-season  conditions.  It  was  recog- 
nized that  the  probability  of  sampling  the  extreme  ranges  of  moisture 
content  during  only  one  sample  period  each  within  the  dry  and  wet  sea- 
sons was  very  remote.  Therefore,  it  was  decided  that  meteorological 
data  from  which  soil  moisture  content  can  be  derived  should  be  collected 

at  specific  locations  for  at  least  a complete  yearly  cycle.  The  meteo- 

3 

rological  data  collected  are  compatible  with  the  WES  computer  model  for 
predicting  soil  moisture  in  the  0-  to  15-  and  15-  to  30-cm  layers  and 


include : 

a. 

Air  temperature,  deg  C. 

b. 

Wind  speed,  m/sec. 

£.• 

Wind  direction,  deg. 

d. 

Precipitation  (rainfall,  snowfall),  cm. 

77-  The 

meteorological,  data  were  measured  by  AEPG  personnel  at 

three  locations  within  the  East  Range  (see  fig.  U),  and  were  recorded 
continuously  during  the  period  July  197-1  to  the  time  of  this  study.  The 
data  were  recorded  on  paper  strip  charts  and ; therefore,  must  be  manu- 
ally extracted  and  put  in  a computer-compatible  format  before  soil  mois- 
ture data  can  be  obtained.  This  work  is  not  within  the  scope  of  this 
study;  however,  data  were  extracted  from  the  records  and  tabulated  for 
the  times  of  the  seismic  response  tests  for  each  of  the  sites  and  are 
included  in  table  1. 

78.  As  stated  above,  the  meteorological  data  were  obtained  to 
provide  a means  for  estimating  the  range  of  soil  moisture  contents  at 
the  East  Range  SID  test  sites.  A secondary,  but  important,  use  of  the 
meteorological  data  is  the  interpretation  of  causes  for  high  ambient 
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seismic  noise  caused  by  meteorological  phenomena.  For  example,  as 
stated  in  paragraphs  4 and  74,  wind  blowing  against  vegetation  can  cause 
considerable  seismic  noise,  and  for  this  reason,  a cursory  description 
of  vegetation  was  made  in  the  field  program.  The  vegetation  factors  for 
which  data  have  been  obtained  include : 
a_.  Tree  type  (species). 

b_.  Stem  base  diameter(s),  cm;  and  number  of  stem  components 
per  plant  assemblage. 
c_.  Tree  height,  cm. 
d..  Average  tree  spacing,  cm. 

e_.  Tree  crown  characteristics  (diameter  and  depth),  cm. 
f_.  Depth,  cm;  and  density  of  ground  cover. 

Data  collection  procedures  (sensor 
detection  performance  requirements ) 

79 . Soil  moisture-density  and  USCS  classification.  In  support 

of  seismic  response  testing,  soil  samples  for  describing  the  substrate 
were  obtained  by  digging  a pit  to  a depth  of  1 m near  sta  +5  at  sites 
335-342,  347-349,  and  4lQ-4l3.  Also,  for  three  sites  (sites  33b,  338, 
and  340)  a hole  approximately  10  cm  in  diameter  was  hand-augered  to  a 
depth  of  3 m.  Compact  soil  prevented  sampling  beyond  the  1-m  depth  at 
the  other  sites.  Soil  samples  from  which  moisture  content  and  wet  den- 
sity were  determined  were  obtained  at  the  surface  and  at  depths  of 
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tional  depths  of  150,  200,  and  300  cm,  respectively.  Bulk  samples  were 

taken  at  various  depths  (as  specified  by  the  field  geologist)  to  obtain 
USCS  classification  for  eat  i .identifiable  soil  layer.  Soil  moisture 
content,  density,  and  classification  data  for  the  wet  and  dry  seasons 
for  the  surface  layer  and  layers  1 and  2 of  the  various  sites  (335-342, 
347-349,  and  410-413)  are  listed  in  table  2.  A comparison  of  the  mois- 
ture data  for  the  wet  and  dry  seasons  is  given  in  table  3.  The  surface 
soil  plasticity  and  liquid  limit  data  (Atterberg  limits)  for  the  sites 
within  the  hand-emplacement  test  area  are  compared  in  fig.  l4. 

80.  CPR  measurements . Cone  penetrometer  readings  (a  measure  of 
soil  strength)  were  obtained  at  each  site,  at  six  stations  (0,  +10,  +20, 


Fig.  l4.  Surface  coil  plasticity  and  liquid  limit  data  for 
10  sites  within  the  hand-emplacement  test  area 

+30,  +40 , and  +50)  along  the  line  at  the  surface  and  at  2.5-cm  incre- 
ments to  a depth  of  15  cm,  and  then  at  7-5-cm  increments  to  a depth  of 
45  cm  (when  possible ) . The  cone  penetrometer  data  for  the  wet  and  dry 
seasons  for  all  sites  except  the  CSEWS  airdrop  area  are  given  in 
tables  2 and  4.  The  CPK  data  for  the  CSEWS  area  for  each  sensor's  loca- 
tion are  given  in  table  5- 

8l.  Compression  wave  velocity  and  thickness  of  layers.  To  obtain 
a measure  of  the  compression  wave  velocity  and  thickness  of  the  refrac- 
tion layers , a refraction  seismic  survey  was  conducted  at  the  various 
sites  with  the  WES  modified  Geo-Space  Corporation  GT-2B,  12-channel 
portable  seismograph,  (fig.  15).  The  refraction  survey  consisted  of  em- 
placing 11  geophones  at  specified  intervals  (3m  for  the  dry  season  and 
5 m for  the  wet  season)  between  two  stations  (0  and  +30,  and  0 and  +50) 
and  recording  the  arrival  times  of  the  initial  seismic  impulse  at  each 
respective  geophone.  The  seismic  signal  was  generated  by  a hammer  blow 
at  each  end  of  the  refraction  line  (i.e.  sta  0 and  +30,  or  sta  0 and 
+50).  Refraction  data  (compression  wave  velocities  and  thicknesses  of 
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Fig.  15.  Instrumentation  used  in  seismic  refraction  tests 


layers)  for  the  wet-  and  dry-season  conditions  are  summarised  for  the 
surface  and  first  and  second  refraction  layers  in  table  2.  The  refrac- 
tion data  portrayed  in  a time-versus-distance  format  for  the  sites  are 
included  in  Appendix  B,  figs.  B12-B39 . 

82.  Shear  wave  velocity.  For  this  study  the  surface  shear  wave 
velocity  was  considered  equal  to  the  measured  surface  (Rayleigh)  wave 
velocity  (see  paragraph  h3ij  . The  shear  wave  velocity  was  determined  by 
comparing  similar  distinct  wave  forms  from  the  time-displacement  refrac- 
tion records  of  six  successive  geophones.  The  shear  wave  velocity  was 
calculated  by  dividing  the  distance  (3m)  between  successive  geophones 
by  the  time  lag  between  similar  wave  shapes.  The  depth  at  which  this 
calculated  velocity  applies  was  computed  as  one-half  the  wavelength  of 
the  wave  used  in  determining  the  shear  wave  velocity.  (The  wavelength 
is  the  wave  period  multiplied  by  the  wave  velocity.)  By  using  several 
successive  geophones,  various  depths  and  their  corresponding  velocities 
were  determined.  From  these  data  the  shear  wave  velocity  for  various 
selected  depths,  such  as  the  first-layer  refraction  depth  (or 
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thickness),  was  computed.  The  shear  wave  velocity  data  for  the  various 
sites  are  contained  in  table  2. 

83.  Topography . Ground  surface  configuration  data  were  obtained 
by  surveying  a detailed  profile  along  each  110-m  drop-hammer  and  man- 
walking test  site  and  each  1100-m  vehicle  test  site.  Since  the  drop- 
hammer  and  man-walking  sites  selected  were  in  relatively  flat  terrain, 
there  were  very  few  terrain  surface  changes  along  the  110-m  line.  In 
the  cross-country  vehicle  test  sites  (sites  3^3  and  344),  a few  minor 
elevation  changes  occurred  along  the  1100-m  lines,  but  these  changes 
were  not  abrupt;  therefore,  the  effect  on  the  seismic  response  would 
probably  be  quite  small.  The  surface  configuration  data  for  the  drop- 
hammer  and  man-walking  sites  and  the  vehicle  sites  have  not  been  in- 
cluded in  this  report. 

84.  Vegetation.  The  vegetation  surrounding  each  man-walking  and 
drop-hammer  site  was  described  by  determining  the  plant  type  (species), 
stem  base  diameters,  heights,  spacing,  number  of  stems  pei-  plant,  and 
crown  characteristics  (diameter  auu  depth)  of  the  dominant  and  codomi- 
nant plants  in  the  population,  and  determining  the  amount  of  grass  or 
other  low-growing  plants  that  occupied  the  site.  Also,  the  height  of 
ground  cover  along  the  110-m  line  was  measured  and  recorded.  The  vege- 
tation data  describing  the  various  sites  are  listed  in  table  6.  General 
terrain  descriptions  of  the  various  sites  are  given  in  paragraphs  2-13, 
Appendix  B. 


Data  collection  procedures 
(sensor-implantation  requirements ) 

85.  In  support  of  providing  input  data  for  the  theoretical,  equa- 
tions that  predict  th<  depth  of  penetration  of  air-delivered  SID's, 
special  soil  samples  were  obtained  at  depths  of  30,  90,  and  120  cm  from 
two  locations  within  each  of  airdrop  areas  1,  4,  and  5,  respectively. 
The  locations  of  the  soil  pits  in  each  of  the  three  areas  are  shown  in 
fig.  l6.  The  equations  used  to  predict  the  depths  of  penetration  are 
described  in  Appendix  C.  The  input  requirements  for  the  sensor- 
implantation  equations  are  listed  in  paragraph  46b_. 


86.  Data  on  E 


Pp/P  , 


and  Cp  were  obtained  from 


48 


the  results  of  special  uniaxial  strain  and.  triaxial  compression  tests. 
These  special  laboratory  tests  were  performed  on  a soil  sample  from  the 
30-m  depth  at  each  site,  at  two  moisture  contents  believed  to  be  repre- 
sentative of  the  wet  and  dry  conditions  at  the  site.  The  laboratory 
tests  (uniaxial  and  triaxial)  and  data  reduction  procedures  are  docu- 
mented in  detail  in  references  4-8.  The  results  of  the  laboratory  tests 
are  briefly  discussed  below. 

87.  The  uniaxial  strain  tests  were  conducted  by  rapidly  loading  a 
laterally  constrained  soil  sample  and  measuring  its  axial  strain-time 
response.  The  results  from  a typical  test  are  shown  in  fig.  17  and  con- 
sist of  stress-time,  strain-time,  and  stress-strain  relations. 

88.  The  triaxial  tests  were  conducted  by  rapidly  loading  a con- 
fined (confining  pressures  from  0 to  approximately  68,950  kPa,  or 
10,000  psi)  soil  sample  to  failure  (i.e.  failure  in  15-40  msec).  The 
results  from  a typical  test,  shown  in  figs.  18-20,  consist  of  (a)  rela- 
tions of  deviator  stress  (i.e.  vertical  stress-radial  stress)  versus 
axial  strain  (fig.  18 ) , (b)  relations  of  mean  normal  stress  to  maximum 
stress  difference  (fig.  19),  and  (c)  relations  of  shear  to  normal  stress 
(fig.  20).  The  input  parameters  for  the  penetration  equations  were  de- 
rived from  the  laboratory  test  results  as  described  below. 

89.  Initial  Young's  modulus  of  elasticity,  E . This  parameter 
was  obtained  from  the  initial  slope  of  the  deviator  stress-axial  strain 
relation  as  shown  in  the  example  in  fig.  l8.  Mathematically,  E equals 

u.(o  O )/u.£  < lil  ‘tlic  0Xcunp?i_6  "fcliS  a. VC I'cLgG  dope  )Tor*  "tVO  CGll 

samples  is  given. 

90.  Compressibility  parameter,  p^/p  • The  compressibility  pa- 
rameter was  determined  from  the  equation  of  continuity  as  follows : 


= exp(e  ) 
P v 


(1) 


where 


p = current  wet  density  of  the  substrate  material 
P 

p = initial  wet  density  of  the  substrate  material 

e = vertical  or  axial  strain 
v 
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Fig.  3,9.  Dynamic  yield  surface  from  triaxial  compression  test 
airdrop  area  1,  test  pit  5 
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Fig.  20.  Mohr's  diagram,  airdrop  area  1,  test  pit  6,  dry  season 


The  experimental  values  of  e were  detennined  from  the  uniaxial 
strain-stress  strain  data  as  shown  in  the  example  in  fig.  IT. 

91.  Ultimate  shear  strength.  The  ultimate  shear  strength  (o^) 
was  computed  by  the  equation  . y 

(\  - »r) 


a = 
u 


max 


V3" 


(2) 


where 


= ultimate  shear  strength  of  soil 

a = axial  stress  used  in  the  triaxial  compression  test 
a 

o = radial  stress  or  confining  pressure  used  in  tri- 
axial compression  test 


(°a  - °r) 


= stress  difference  at  failure 


max 


The  maximum  stress  difference  (°a  ~ °rj  used  in  eq  2 was  determined 
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stress  difference  of  the  soil  is  assumed  to  be  the  maximum  soil  stress 
that  would  be  experienced  by  the  substrate  material  and  the  ADSID/S  dur- 
ing a penetration  event. 

92.  Drag  coefficient.  The  drag  coefficient  was  computed  using 
eq  C2  (Appendix  C),  the  ADSID/S  dimension,  and  the  computed  coefficient 
of  friction  (y)  between  the  projectile  and  the  soil.  The  coefficient  y 
was  obtained  as  follows : The  dynamic  angle  of  internal  friction  of  the 

soil  <f  was  determined  from  the  Mohr  diagram  (i.e.  shear  stress  versus 
normal  stress)  shown  in  fig.  20.  The  relation  shown  in  fig.  21  was  then 
entered  using  the  relative  amounts  of  sand  and  clay  in  the  soil  sample 
as  portrayed,  by  the  grain-size  curve  for  the  site  and  a value  of  6/ <f> 
was  determined  for  the  dynamic  case , where  6 is  the  angle  of  friction 
between  the  projectile  and  the  soil.  Since  the  angle  of  internal  fric- 
tion <}>  is  known  (fig.  20),  6 can  be  computed  and  the  tangent  of  6 

is  the  coefficient  of  friction  y . The  values  of  y and  the  ADSID/S 
dimension  in  conjunction  with  eq  C2  (Appendix  C)  were  then  used  to  com- 
pute the  drag  coefficient  for  each  airdrop  area.  The  values  of  y 

and  the  data  used  to  obtain  them  are  listed  in  table  7. 
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Fig.  21.  Skir  friction  betweer  smooth  steel  and  soil 
(static  and  dynamic) 


PART  III:  ANALYSIS  OF  DATA 


93.  To  reach  the  objectives  identified  in  paragraph  37,  the  anal- 
ysis of  the  data  was  directed  toward  answering  the  following  specific 
questions , 

a,.  Wnat  are  the  seismic  response  characteristics  of  the  SID 
test  areas  at  Fort  Huachuca?  How  do  these  characteris- 
tics vary  in  time  and  space? 

b_.  What  environmental  factors  should  be  measured  in  support 
' of  SID  testing,  i.e.  both  sensor  detection  performance 
and  implantation  performance,  to  provide  a basis  for  the 
interpretation  of  the  results?  What  are  the  temporal  and 
spatial  variations  in  terrain  factors  at  Fort  Huachuca? 

c . Do  the  airdrop  test  sites  at  Fort  Huachuca  provide  test 
environments  that  will  allow  evaluation  of  ADSID/S  im- 
plantation over  a wide  range  of  soil  conditions?  What 
are  the  impact  stresses  imposed  on  ADSID/3's  at  the  Fort 
Huachuca  airdrop  test  areas? 

d.  What  analytical  schemes  can  he  used  to  extrapolate  SID 
test  results  (both  sensor  detection  performance  and  sen- 
sor implantation  performance ) from  one  set  of  terrain 
(site)  conditions  to  another? 

94.  This  part  of  the  report  is  directed  toward  the  analysis  of 
the  data  on  a sixe-to-site  basis;  whereas  Part  IV  addresses  the  mapping, 

"tile  SpfcLU  i 5l1_  V6ti‘i 8/bi  Oil  OT  tliO  i"  O 4-0  V culb  Vi  T OniuOIl'u  ? 3013331?.  C TO  — 

sponse,  and  SID  performance  (detecmion  distance)  data. 

Seismic  Response  Data 


Selection  of  seismic 
response  descriptors 

95.  As  previously  stated,  the  seismic  response  is  the  ground 
motion  resulting  from  the  propagation  of  seismic  waves  generated  by  a 
target  or  some  background  disturbance , The  seismic  waves  generated  by 
various  targets  (drop  hammer,  walking  man,  moving  vehicle,  etc.)  used 
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in  this  study  were  measured  with  scientific  geophones  and  recorded  in 
analog  format.  Since  the  analog  seismic  signals  proved  unwieldy  for 
site-to-site  comparison  and  mapping,  descriptors  of  the  signal  meaning- 
ful to  SID  detection  performance  testing  ve.re  selected  for  analysis. 
These  descriptors  have  been  used  in  nearly  all  the  seismic  response 
studies  conducted  at  WES  including  those  conducted  for  TECOM,  i.e. 
studies  conducted  for  the  Tropic  Test  Center  and  the  Airborne  Communica- 
tion and  Electronics  Board  in  the  Panama  Canal  Zone"^  and  at  Fort 

14  pc 

Bragg,  North  Carolina.  5 ' The  descriptors  selected  for  analysis  are 
PPV,  PSPV,  and  f • all  are  defined  in  paragraph  43. 

9 6.  As  used  herein,  PPV  is  the  maximum  (peak.)  positive  particle 
velocity  (or  amplitude)  that  occurred  in  the  seismic  wave  for  each  re- 
spective distance.  PSPV  (peal  summed  particle  velocity)  can  best  be 
described  by  the  techniques  vised  to  obtain  it.  A fast  Fourier  computer 
program  was  used  to  convert  the  digital  data  for  each  record  from  the 
time  domain  to  the  frequency  domain,  resulting  in  a curve  of  particle 
velocity  versus  frequency  (f).  A curve  of  summed  particle  velocity 
(SPV)  versus  f was  obtained  by  applying  a simulated  T-Hz  filter  band 
to  the  fast  Fourier  output,  data.  The  simulated  filter  used  the  equation 


SPV.  = 7 

l 


/x.  _ + X.  , + X.  + X.  _ + X.x0 
1-2  l-l  l l+l  i+2 


where  is  the  particle  velocity  at  a frequency  i , to  average  the 

particle  velocity  amplitudes  occurring  within  the  band  and  to  determine 
the  total  response  that  occurred  within  the  band.  This  process  for  ob- 
taining PSPV  is  illustrated  in  fig.  22.  The  PSPV  is  the  maximum  value 
(amplitude)  of  SPV  that  occurs  on  the  SPV-f  curve.  The  frequency  used 
was  that  at  which  PSPV  occurred;  this  is  the  frequency  at  which  most  of 
the  energy  from  the  seismic  wave  is  propagated  at  a given  distance  from 
the  signal  source. 

97.  In  addition  to  PPV,  PSPV,  and  f two  other  seismic  descrip- 
tors, nns  velocity  (amplitude)  and  span  were  used.  They  are  also  de- 
fined in  paragraph  43. 


summed  Particle  velocity,  cm/SEC  x to' 3 particle  velocity,  cm/sfc  x ic3  Particle  velocity,  cm/SEC  x io‘‘ 
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C.  SUMMED  PARTICLE  VELOCITY  VS  FREQUENCY  DATA 


Fig.  22.  Process  for  obtaining  summed  particle 
velocity-frequency  data  from  analog  data 


Reproducibility  of  seismic  descriptors 

98,  A prerequisite  to  a meaningful  comparison  of  the  seismic 
characteristics  of  a specific  site  at  different  times,  such  as  in  vet 
and  dry  seasons,  or  of  the  seismic  responses  of  different  sites  is  the 
inherent  variability  in  the  signal  generator,  measurement  and  reduction 
equipment,  and  procedures  used  to  obtain  the  seismic  descriptors  (i.e. 
PPV,  PSPV,  and  f ).  The  controlled  source,  or  drop  hammer,  has  been 
used  consistently  as  the  signal  generator  for  all  TECOM  seismic  sensor 
studies.  The  variability  of  the  equipment  and  procedures  has  been 
studied  for  selected  sites  in  the  Panama  Canal  Zone  and  at  Fort  Bragg 
(paragraph  95). 

99-  The  results  of  the  studies  in  Panama  show  that  values  of  the 
three  seismic  response  parameters  (i.e.  PPV,  PSPV,  and  f ),  measured 
a.t  all  distances  tested  (10,  20,  30,  1+0,  and  60  m) , could  be  reproduced 
85  percent  of  the  time  within  jd.O  percent  of  their  mean  values.  The  re- 
sults of  the  drop-hammer  reproducibility  tests  at  Fort  Bragg  have  been 
summarized  and  are  included  in  fig.  23.  These  results,  especially  the 
PPV  and  PSPV  values,  show  a decrease  in  variation  with  distance.  To  il- 
lustrate the  percentage  of  variation  in  these  data,  the  same  three  sites 
(sites  30l,  1+92,  and  1+93)  were  selected  for  determining  the  standard 
deviation/mean  (SD/iu)  values.  These  SD/'m  percentages  are  tabulated 


below. 
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100.  The  valuations  in  the  Fort  Lragg  data  are  slightly  greater 
than  those  in  the  Panama  data.  In  summary,  it  appears  the  data  obtained 
with  the  drop  hammer  are  reproducible  within  about  15  percent  of  the 
mean,  most  of  the  time j therefore,  the  drop-hammer  data  should  provide  a 
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Fig.  23.  Comparison  of  minimum  values,  maximum  values,  arithmetic  means, 
and  standard  deviations  for  the  general  seismic  descriptors  (PPV,  PSPV, 
tmd  f) , Port  Bragg  repeatability  tests  (wet-season  conditions) 
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useful  basis  for  comparing  the  seismic  response  of  one  site  with  that  of 
another  and  also  for  comparing  the  seismic  responses  at  the  same  site 
under  varying  terrain  conditions.  Since  the  variability  in  the  seismic 
response  data  is  an  important  consideration  in  SID  testing,  at  least 
three  repetitions  of  each  seismic  measurement  should  be  made  in  future 


seismic  sensor-  studies , 


Reduction  of  dry-season 
seismic  response  data 


101.  Field-data-eollection  programs  were  conducted  in  August  1971 
(wet  season)  and  March  1972  (dry  season).  Seismic  l’esponse  data  result- 
ing from  drop-hammer,  roan-walking , and  vehicle  targets  were  obtained 
during  both  periods.  The  data  collected  in  the  wet  season  are  believed 


to  be  in  error  as  a result  of  an  undetected  instrumentation  calibration 


problem  and.  are  not  suitable  for  analysis;  however,  the  wet-season  drop- 
hammer  data  were  reduced  and  are  presented  in  table  8.  The  values  of 
PPV  and  PSPV  are  consistently  lower  than  expected  and  no  further  analy- 
sis was  performed. 

102.  From  the  dry-season  seismic  response  data,  typical  records 


for  each  site  were  selected  as  follows; 


Drop  hammer  and  man  walking.  Records  at  distances  of 


5,  10,  20,  30,  40 , and  60  m from  the  energy  source. 

b.  M151  vehicle.  Records  at  distances  of  50,  100,  150, 
200,  250,  300,  and  400  m from  the  energy  source. 

c.  M577  vehicle.  Records  at  distances  of  100,  200,  300, 
400 , and  500  m from  the  energy  source . 

d.  Ambient  and  induced  noise.  One  sample  record 
(0.682  sec). 

103-  The  dr.v-season  field-measured  seismic  response  (analog  sig- 
nature) data  obtained  from  the  man-walking,  drop-hammer,  vehicle,  and 
background-  (ambient  and  induced)  noise  tests  were  converted  to  digital 
form  by  computer  at  1500  points  per  second.  The  digitized  records  were 
then  processed  through  a fast  Fourier  transform  program  to  convert  the 


digital  data  from  the  time  domain  to  the  frequency  domain.  PPV,  PSPV, 


f , and  rms  velocity  were  then  computed  and  tabulated  for  the 


drop-hammer,  man-walking,  and.  vehicle  energy  sources  (tables  8-10,  re- 
spectively). Plots  of  SPV-f  relations  were  derived  for  selected  records 
from  the  vehicle  and  background-noise  tests.  These  data  and  their  anal- 
yses are  discussed  in  the  following  paragraphs. 

Site-to-site  variation  in 
dry-season  seismic  response  data 

104.  Graphs  portraying  the  seismic  response  parameters  PPV,  PSPV , 
and  f as  a function  of  distance  from  the  energy  source  are  presented 
in  figs.  24-26,  respectively.  The  relations,  are  discussed  below. 

105.  PPV-distance  data.  The  PPV  data  (fig.  24),  in  general,  show 
a fairly  good  correspondence  among  sites  within  the  hand-emplacement 
area  (i.e.  sites  335-34.2,  348,  and  349).  An  exception  to  this  trend  is 
the  curve  for  site  342,  which  shows  a very  small  peak  amplitude  of 

0.01  x 10  cm/sec  at  a distance  of  60  m.  The  curves  for  sites  within 
the  airdrop  areas  (i.e.  sites  34?  and  4lO-4l3)  show  considerable  data 
variation.  The  curve  for  Willeox  Playa  site  (site  413)  is  a great 
deal  different  from  those  for  sites  within  the  four  other  airdrop 
areas . 

10b.  PSPV-distance  data.  The  PSPV  data  (fig.  25)  vary  consider- 
ably for  distances  less  than  approximately  4o  m,  and  somewhat  less  for 
distances  of  40-60  m.  Data  for  the  10  sites  within  the  hand-emplacement 
area  and  airdrop  sites  4lO-4l2  show  a fairly  good  correspondence.  The 
Willeox  Playa  site  (site  4] 3)  shows  considerable  variation  from  the 
other  sites  and  also  shows  very  little  attenuation  in  the  seismic  wave 
between  5 and  60  m. 

107.  f-distance  data.  The  sum-frequency  data  (fig.  26)  portrayed 
for  the  sites  within  the  hand-emplacement  area  vary  considerably  (20- 
50  Hz)  at  close  distances,  but  are  in  relatively  good  agreement  for 
distances  of  30  m and  greater.  Sum-frequency  data  for  airdrop  sites 
4l0-4l2  are  not  too  much  different  from  the  data  for  the  sites 
within  the  hand-emplacement  area;  however,  there  is  more  variation  for 
certain  frequencies  (i.e.  30  and  40  Hz).  The  Willeox  Playa  site  again 
is  quite  different  from  the  other  sites  and  shows  a low  and  constant 
predominant  sum  frequency  of  7-9  Hz,  In  addition  to  the  variation  of 
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!ii.  Site-to-site  variation  in  PPV  at  various  distances  from 
the  drop-hammer  energy  source,  dry-season  conditions 


sum-frequency  data  among  the  sites,  there  is  also  considerable  variation 
in  sum-frequency  data  within  each  site,  with  the  exception  of  site  413- 

Variation  in  dry- 
season  ambient-noise  data 

108.  The  seismic  response  resulting  from  ambient  (background) 
noise  was  recorded  during  normal  working  hours  at  each  site  as  discussed 
in  paragraph  74 . One  sample  record  (0.682  sec)  each  was  selected  for 
sites  335-337,  339,  34l,  3^2,  and  3^9  to  study  the  range  of  varia- 
tion in  ambient  noise  found  in  the  hand-emplacement  test  area.  These 
records  were  used  to  generate  SPV-f  relations  (fig.  27)  to  depict  the 
amplitude  and  frequency  content  of  the  ambient  seismic  noise.  The  ambi- 
ent noise  present  within  the  hand-emplacement  area  varied  between 

-3  -3 

0.008  x 10  and  0.09  x 10  cm/sec.  in  SPV  amplitude  with  the  maximum 
amplitudes  generally  falling  within  a low-frequency  (1-20  Hz)  bandwidth. 
Although  the  background  noise  was  measured  in  diverse  locations  in  the 
hand— emplacement  are  a under  various  wind  ccnu.iuj.wns  turn  au  uiiierent 
times  of  the  day,  the  SPV-f  data  are  quite  similar  in  that  they  show 
relatively  low  background  noise  levels.  Wind  conditions  observed  at 
site  341  (see  fig.  27e ) were  the  highest  measured  (4. 1-6. 6 m/sec)  during 
the  background-noise  tests,  and  the  ambient-noise  level  at  all  frequen- 
cies is  higher  for  that  site  than  for  the  other  six  sites.  Even  so,  the 
noise  level  is  well  below  that  which  would  cause  difficulty  in  testing 
sensors  with  threshold  levels  and  a frequency  response  similar  to  those 
of  the  Phase  III  MINISID.  In  summary,  it  appears  that  the  hand- 
emplacement  area  should  provide  a seismically  quiet  area  for  SID  test- 
ing, insofar  as  ambient  noise  is  concerned. 

Variation  in  dry- 
season  induced-noise  data 

109.  SID  testing  will  always  be  affected  by  seismic  noises  in- 
duced by  activities  extraneous  to  the  test.  During  the  field  program, 
seismic  responses  caused  by  aircraft  (a  Huey  helicopter  and  a single- 
and  a double-engine  light  aircraft ) and  by  a freight  train  traveling 
across  Willcox  Playa  were  measured.  The  SPV-f  relations  derived  from  a 
selected  number  of  measurements  are  presented  in  fig.  28.  The  data 
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Fig.  27.  Dry-season  ambient -noise  data 
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displayed  in  fig.  28  were  obtained  on  an  opportunity  basis;  therefore, 
altitudes  and  distances  from  the  energy  source  were  not  measured.  As 
expected,  the  results  show  that  the  induced  noises  increase  the  general 
level  of  background  noise  and,  in  some  cases,  increase  it  markedly  as  a 
function  of  frequency.  Fcr  example,  the  curve  for  the  helicopter 
(fig,  28a)  shows  several  peeks  in  amplitude  occurring  at  frequency  in- 
tervals of  about  12  Hz.  The  Piper  Cub  (fig.  28b)  did  not  appear  to 
materially  affect  the  noise  level  at  frequencies  below  about  TO  Ha  (com- 
pare figs.  28a  and  28b),  but  increased  it  sharply  at  88  Hz.  The  noise 
level  of  the  Twin  Cessna  (fig.  28c)  showed  a near  constant  amplitude 
(0.03  * 10  ' cm/sec)  for  frequencies  between  0 and  100  Hz,  with  the  ex- 
ception of  the  small  peak,  at  approximately  88  Hz.  At  about  115  Hz  an 
abrupt  change  in  SPV  amplitude  occurred.  The  train  traveling  across 
Willcox  Playa  (fig.  28d)  caused  a large  increase  in  seismic  noise 
(0.3  x 10  cm/sec)  from  0 to  lb  Hz,  and  considerable  variation  in  noise 
levels  for  frequencies  between  l6  and  110  Kz. 

110.  In  summary,  induced  background  noise  will  most  likely  occur 
during  SID  testing  at  Fort  Huachuca.  While  the  induced  noise  data  pre- 
sented herein  are  representative  of  only  selected  cases  and  do  not  nec- 
essarily reflect  the  maximum  amplitude  and  predominant  frequencies  that 
might  be  encountered  at  Fort  Huachuca,  the  amplitudes  of  the  induced 
noises  will  be  above  those  of  the  ambient  noise,  and  the  frequencies  and 
amplitudes  of  induced  noise  will  span  a wide  range  of  values.  Since  the 
ambient  and  induced  background  noises  increase  or,  in  some  cases,  de- 
crease the  sensitivity  -'f  certain  types  of  SID's,  it  is  recommended  that 
the  background  noises  be  measured  during  the  times  of  SID  testing. 
Seismic  response  from  vehicles 

111.  Seismic  response  relations  in  terms  of  SPV  versus  f were 
obtained  for  the  two  vehicles  (M151  and  M377 ) that  were  tested  at  the 
two  cross-country  sites  (3^  and  3h5 ) in  the  hand-emplacement  area.  The 
vehicles  were  run  at  two  speeds,  8 and  32  km/hr.  These  data  are  pre- 
sented in  fig.  29.  The  M577  traveling  at  32  km./hr  generated  the  highest 
SPV  values  for  the  frequencies  shown.  Also,  except  for  some  frequencies 
from  70  to  78  Hz  for  site  3^5  (fig.  29b),  the  M151  traveling  at  8 km/hr 


generated  the  least  amount  of  seismic  response.  It  is  interesting  to 
note  that  for  both  sites,  there  is  a good  deal  of  similarity  between  the 
response  generated  by  the  M577  traveling  at  8 km/hr  and  the  M151  travel- 
ing at  32  km/hr.  The  genera],  similarity  in  the  SPV-f  data  for  sites  34k 
and  345  for  each  of  the  two  vehicles  suggests  that  the  terrain  condi- 
tions at  both  sites  are  quite  similar.  Some  of  the  more  subtle  differ- 
ences noted  among  the  3PV-f  curves  for  the  two  sites  cannot  be  fully  ex- 
plained, e.g.  the  lowest  SPV  value  for  the  MI51  traveling  at  8 km/hr 
occurs  at  48  Hz  at  site  344,  and  at  32  Hz  at  site  34 1 . However,  differ- 
ences in  the  roughness  of  the  surface  over  which  the  vehicles  traveled 
could  have  accounted  for  differences  of  this  kind  and  magnitude. 

112.  The  results  illustrate  clearly  that  the  seismic  responses 
emanating  from  vehicles  are  a function  of  both  vehicle  type  and  speed. 

In  addition,  it  is  well  known  that  the  type  of  terrain,  i.e.  both  sur- 
face and  substrate  conditions,  also  affects  the  seismic  responses  of 
vehicles,  although  in  this  case  the  two  sites  were  quite  similar.  The 
implication  of  this  is  that  if  SID  results  from  vehicle  tests  are  to  be 
extrapolated  to  conditions  other  than  those  occurring  during  the  tests, 
the  effects  of  vehicle  speed  and  terrain  must  be  quantitatively  known. 

Experimental  MINISID  and  ADSID/S  Test  Data 

MINISID  data. 

113.  Detection  distances  obtained  at  low-,  medium-,  and  high-gai.n 
settings  of  a Pnase  III  MINISID  are  presented  in  table  11.  Fig.  30 
shows  graphically  the  we l-  and  dry-season  MINISID  detection  values  at 
medium  gain . A point  especially  worthy  of  mention  is  that  two  different 
men  made  the  walking  tests  in  the  vet  and  dry  seasons;  therefore,  the 
data  for  these  tests  for  the  two  seasons  may  or  may  not  be  directly 
comparable . 

11 4.  The  detection  distances  presented  in  fig.  JO  chow  a differ- 
ence of  1-3  in  for  five  sites  (337-339,  342,  and  34b) > a difference 

of  5-7  m for  two  sites  (34l  and  3,J7 ) , and  a diffr- rerice  of  18-20  m for 
two  sites  (386  and  348).  For  the  nine  cases,  there  seems  to  be  no 


particular  trend  toward  either  season's  allowing  a consistently  greater 
or  lesser  detect!-  ige. 

ADSID/G  data 

115-  The  results  of  the  experimental  ADSID/S  tests  conducted  by 
CSEWS  are  presented  in  table  12-  These  tests  were  conducted  during 
August  1972  and  should  be  representative  of  wet-season  conditions  at  the 
CSEWS  test  site. 

116 . The  results,  in  general,  show  that  most  penetration  lengths 
(i.e.  depth  of  penetration  measured  along  the  longitudinal  axis  of  the 
sensor)  ranged  between  20  and  35  cm,  although  penetrations  as  deep  as 
50  cm  and  as  shallow  as  15  cm  were  obtained.  The  emplacement  configura- 
tions of  two  of  the  air-implanted  ADSID/S's  at  the  CSEWS  test  area  are 
shown  in  fig.  31. 


a.  Impact  angle  50  deg  (off 
verti cal ) 

Fig.  31 . ADSlD/S  implanted  in  CSEWS  airdrop  test  area 

Environmental  Data 


Ilf.  The  environmental  data  obtained  at  Fort  Huachuca.  test  sites 
were  those  directly  applicable  to  the  prediction  of  SID  detection  per- 
formance and  sensor-implantation  performance  as  discussed  previously. 

The  data,  along  with  the  theoretical  models  (developed  nr  in  development 
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at  WES),  provide  the  input  and.  analytical  tools  needs  for  interpreta- 
tion and  extrapolation  of  SID  test  results  from  site  to  site  and  from 
region  to  region.  The  environmental  data  collected  for  input  to  SID 
detection  performance  and  air- implantation  performance  models  are 
discussed  below. 


Data  for  SID  detection  per- 
formance  predictions  and  evaluations 

118.  The  environmental  data  measured  in  support  of  SID  detection 

performance  testing  and  prediction  were  obtained  during  August  1971  (wet 
season)  and  March  1972  (dry  season).  So  that  these  data  could  be  input 
directly  into  the  SID  detection  performance  prediction  model,  they  were 
.reduced  and  summarized  by  soil  layers:  surface,  refraction  layer  1,  and 

refraction  layer  2. 

119.  Surface  layer.  The  surface,  or  uppermost,  layer  of  soil  is 

the  thin  zone  of  material  that  comes  in  contact,  with  the  seismic  energy 
source  (drop  hammer,  footstep,  vehicle  track  or  wheel,  and  so  on)  end 
the  SID  geophore . Although  the  phenomena  dealing  with  the  interaction 
of  the  energy  source  and  the  substrate  have  not  been  adequately  quanti- 
fied, the  surface  terrain  factors  considered  important  in  the  transfer 
of  seismic  energy  from  the  source  to  the  ground  and  from  the  ground  to 
the  SJD  geophone  are:  \ 

a.  Soil  type  (iJSCS). 

■J.b.  Soil  particle  size,  percent  finer  by  weight. 

c_.  Moisture  content , percent-  • 

•? 

d_.  Soil,  vet  density,  g/cnt. 

e_.  Cone  penetration  resistance,  CPR  (a  measure  of  soil 
strength,  kPa. 

f_.  Height  (cm)  and  density  of  vegetation  ground,  cover. 

g_.  Geometry  of  ground  surface  l roughness  ) . 

120.  The  data  for  the  factors  listed  above  with  the  exception  of 
surface  geometry  data  are  given  in  tables  2-6.  Of  the  factors  .listed, 
only  soil  type,  soi.c  p article  size,  and  surface  geometry  can  be  consid- 
ered static  or  constant  with  time.  Moisture  content,  wet  density,  CPR, 
and  height  and  density  of  vegetation  change  as  a function  of  site  and 


weather  conditions.  From  table  2,  it  can  be  seen  that  the  majority  of 
surface  soils  a.t  the  various  sites,  in  general,  exhibit  the  same  charac- 
teristic gradation  and  most  are  classified  SC.  Exceptions  to  this  trend 
are  sites  339,  390 , and  3^+9,  where  the  surface  soils,  are  classified  £M, 
and  sites  9 10  and  4ll,  which  contain  53.5  and  69.4  percent  silt  and 
clay,  respectively,  and  are  classified  CL.  Particle-sine  .distributions 
for  surface  soils  at  Fort  Huachuca  sites  generally  are  quite  similar; 
however,  there  were  considerably  more  small . surface  rocks , at  sites  338, 
4l0,  and  4ll  than  at  the  other  sites.  The  surface  rocks  contained  in 
the  substrate  material  at  these  sites  xrere  not  included  in  the  samples 
used  to  determine  the  particle-siae  distributions  of  the  materials. 

(See  ground  photographs  in  figs.  B2b,  B7b,  and  E8a,  respectively.  Appen- 
dix B.  ) 

121.  Site-to-site  and  season-to-season  comparisons  of  surface 

moisture  content  and  CPR  are  shown  in  figs.  32  and  33,  respectively. 
Similar  comparisons  for  wet  density  are  not  presented  herein  since  wet 

density  is  directly- related  to  moisture  content,  i.e.  p --  (l  + W.K^  . 

where  p = soil  wet  density,  W - soil  moisture  content  (percent),  ancl 

= soil  dry  density.  The  dry-season  site-to-sioe  moisture  values 
within  the  hand-emplacement  area  ranged  from  a low  yf  09  percent  (site 

39l)  to  a high  of  5 percent,  (site  338).  Tne  next  highest  moisture  con- 

tents were  2.8  percent  (site  392)  and  2.2  percent  (site  336);  for  J of 
the  10  sites,  the  dry-season  moisture  content  was  less  than  2 percent. 

122.  As  expected,  the  moisture  values  within  the  hand -emplacement 
area  were  higher  during  th  ,•  we  i,  set  son,  ranging  from  5.2  percent  (site 
336)  to  11. T percent  (site  338).  In  general,  the  surface  moisture  data 
(fig.  32)  show  that  there  yas  more  variation  from  season-tc-season  than 
from  site-to~site . 

123.  Tne  dry-season  surface  CPR  values  (fig.  33)  ranged  from  100 
to  1990  kPa,  with  considerable  variation  among  the  sites.  The  season- 
to-seasen  variation  was  large  for  sites  335,  336,  and  3'i0 ; moderate  for 
sites  338,  392,  and  399;  and  insignificant  for  sites  337,  339,  391,  and 
398.  In  some  cases  the  wet-season  data  exhibited  greater  CPR  values 
than  the  dry-season  data  (e.g.  sites  335,  330,  339,  392,  and  399),  which 
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seems  to  indicate  that  the  addition  of  moisture  created  some  cohesion  in 
the  surface  soil. 

12k.  The  vegetation  ground  cover  (see  table  6)  occupying  the 
sites  within  the  hand-emplacement  area  was  quite  uniform  in  terms  of 
density,  height,  and  percent  of  ground  covered.  Only  a few  isolated 
areas  (near  the  arroyos ) were  without  some  type  of  ground  cover  (sites 
335 » 337,  and  33 9)-  As  a result  of  this  homogeneity,  the  relative  ef- 
fects of  ground  cover  on  SID  testing  were  approximately  the  same  for 
all  sites  within  the  hand-emplacement  area. 

125.  The  geometry  of  the  ground  surface  had  very  little  effect 

on  SID  detection  performance  because  of  the  minor  irregularities  in  sur- 
face configurations  of  the  test  sites  (see  paragraph  83). 

126.  Subsurface  refraction  layers  I and  2.  Refraction  layers  1 
and  2 are  the  first  two  subsurface  layers  containing  different  primary 
(or  compression)  wave  velocities.  These  two  layers  are  determined  by 
techniques  of  refraction  seismology  and  often,  but  do  not  necessarily, 
correspond  to  the  apparent  stratification  in  the  soil  prol'ixe , as 
stated  in  paragraphs  l6-l8. 

127.  Values  of  four  terrain  factors  pertaining  to  refraction 
layers  1 and  2 are  required  as  input  to  the  WES  Theoretical  Seismic 
Propagation  Model, ^ The  four  factors  are  listed  below;  their  wet-  and 
dry-seasons  values  are  contained  in  table  2 . 

a_.  Compression  wave  velocity  of  layer  1 (Vp  ) and  layer  2 
(Vpp),  m/sec. 

b_.  Shear  wave  velocity  of  layer  1 (V  ) and  layer  2 (V  ), 
m/sec . 

£•  Thickness  of  layer  1 (T  ) and  layer  2 (T  ) , cm. 

d_.  Wet  density  of  layer  1 (pp)  and  layer  2 (p  ),  g/cm  . 

128.  Site-to-site  and  season-to-season  comparisons  of  compression 
wave  velocities  in  refraction  layers  1 and  2 are  presented  in  fig.  34. 
For  layer  1,  the  dry-season  velocities  ranged  from  250  to  500  m/sec, 
whereas  the  wet-season  velocities  ranged  from  230  to  455  m/sec,  indi- 
cating that  the  range  of  variation  does  not  change  appreciably  from  the 
dry  to  the  wet  season.  Fig.  34  also  shows  that  the  site-to-site 
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variation  for  layer  1 was  greater  than  the  senson-to-season  variation. 
The  same  trends  hold  for  the  layer  2 compression  wave  velocities,  i.e., 
the  site-to-site  variation  was  greater  than  the  season-to-season.  For 
all  sites  for  which  data  are  presented,  the  compression  wave  velocity  of 
the  second  layer  was  greater  than  the  compression  wave  velocity  of  the 
first  layer.  The  maximum  season-to-season  variation  was  75  m/sec  (site 
338)  for  layer  1 and  98  m/sec  (site  34l)  for  layer  2. 

129.  The  vet-  and  dry-season  values  of  shear  wave  velocities  in 
layers  1 and  2 are  portrayed  in  fig.  35.  From  these  data  it  can  be  seen 
that  the  shear  wave  velocity  in  the  first  layer  was  consistently  less 
than  that  in  the  second  layer.  For  layer  1 the  dry-season  site-to-site 
shear  wave  velocities  ranged  from  129  to  207  m/sec  (sites  336  and  338, 
respectively),  whereas  the  wet-season  shear  wave  velocities  ranged  from 
137  to  182  m/sec  (sites  339  and  348,  respectively).  For  layer  2 the 
dry-season  site-to-site  variation  ranged  from  256  to  371  m/sec  (sites 
337  and  349,  respectively)  and  the  vet-season  variation  ranged  from  245 
to  ^4u  m/sec  (sices  3 3c  and  342 , respectively ) . The  maximum  measured 
season-to-season  variation  (1*0  m/sec)  occuri'ed  at  site  336,  and  the  min- 
imum (12  m/sec)  occurred  at  site  335-  In  general,  the  season-to-season 
shear  wave  velocity  data  are  greater  in  layer  2 than  in  layer  1. 

130.  The  layer  thicknesses  of  refraction  layers  1 and  2 during 
the  wet  and  dry  seasons  are  compared  in  fig.  36.  The  data  for  layer  1 
show  a site-to-site  thickness  range  of  50--140  cm  for  dry-season  condi- 
tions and  a site-to-site  thickness  range  of  55-150  cm  for  wet-season 
conditions.  These  extremes  occurred  at  sites  340  and  339,  respectively. 
The  maximum  season-to-season  thickness  variation  for  layer  1 was  25  cm 
at  site  349 , and  the  minimum  thickness  variation  was  5 cm  at  sites  337 
and  340.  The  site-to-site  thickness  variation  for  layer  2 was  much 
greater  than  for  layer  1 with  values  for  layer  2 ranging  from  135  to 
>1575  cm  and  110  to  >1555  cm  for  the  dry  and  wet  seasons,  respectively. 
These  extremes  occurred  at  sites  336  and  338,  respectively.  Sites  338, 
339,  34l,  and  342  had  very  thick  second  layers  (i.e.  >1500  cm)  in  both 
seasons.  (Values  shown  for  these  sites  were  limited  by  the  length  of 
the  refraction  line;  therefore,  they  were  not  used  to  determine  maximum 
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and  minimum  thicknesses.)  The  maximum  season-to-season  thickness  varia- 
tion in  layer  2 was  100  cm  at  sites  348  and  349.  The  minimum  season-to- 
season  layer  thickness  variation  (15  cm)  also  occurred  in  layer  2 at 
site  340. 

131.  A review  of  the  wet  densities  in  table  2 shows  very  little 
site-to-site  variation  in  the  first  and  second  layers,  i.e.  values 
ranged  from  about  1.7  to  2.0  g/em^.  As  stated  previously,  the  wet  den- 
sity varies  as  a function  of  moisture  content,  and  the  maximum  seasonal 
variation  in  measured  moisture  content  was  4.8  percent  (site  349, 
layer  l).  The  rest  of  the  moisture  content  values  varied  less  than 

3 per.ent.  Therefore,  the  seasonal  variation  in  wet  density  for  the 
majority  of  sites  was  quite  small. 

132.  Summary . The  preceding  paragraphs  have  discussed  the  tempo- 
ral and  spatial  variation  of  the  terrain  factor  either  known  (modeled) 
or  assumed  to  have  some  effect  on  SID  detection  performance.  The  rela- 
tive effects  of  terrain  factors  (and  changes  in  them)  cannot  be  speci- 
fied in  simple  terms,  because  the  propagating  medium-seismic  response 
interactions  are  exceedingly  complex,  and  a practical  and  general,  theo- 
retical solution  requires  extensive  computer  modeling.  A general  solu- 
tion to  this  problem  is  being  sought  at  WES,  and  a first-generation 
model  has  been  developed  (see  reference  17).  It  is  recommended  that 
this  model  be  used  to  perform  sensitivity  analyses  on  the  Fort  Huachuca 
data  to  establish  the  relative  effects  of  the  various  terrain  factors  or. 
seismic  wave  generation  and  propagation.  Once  this  is  accomplished, 
only  those  factors  causing  significant  changes  in  SID  detection  ranges 
would  require  measurement  in  support  of  SID  testing. 

Data  for  sensor  air-implantation  per- 
fomiance  predictions  and  evaluations 

133.  Data  derived  for  use  as  input  to  equations  (Appendix  C)  that 

predict  implantation  performance  of  air-delivered  SID's  axe:  the  ini- 

tial Young's  modulus  of  elasticity,  ultimate  shear  strength,  compressi- 
bility parameter,  soil  mass  density,  and  drag  coei i icient . The  methods 
for  deriving  these  factors  from  the  uniaxial  strain  and.  triaxial  compres- 
sion tests  are  presented  in  paragraphs  07-92.  The  laboratory -obtained 
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data  on  these  factors  for  airdrop  areas  1,  4,  and  5 are  listed  in 
table  13-  Airdrop  areas  2,  3,  and  6 (see  fig.  4 for  location  of  areas  2 
and  3 and  fig.  5 for  the  location  of  area  6)  were  also  selected  for 
stuoy  of  penetration  characteristics  of  air-delivered  SID's,  but  time 
and  funding  constraints  prohibited  any  work  on  these  areas , Airdrop 
area  2 had  been  selected  to  provide  on  area  that  would  be  impenetrable 
because  of  the  numerous  rock  outcrops  and  shallow  soil.  Airdrop  area  6 
had  been  selected  to  present  a target  of  soft  material  to  study  the  ef- 
fects of  deep  implantation  of  air-delivered  SID's;  and  area  3 had  been 
selected  as  a target  of  intermediate  difficulty,  along  with  aii’drop 
areas  1,  4,  and  5.  The  following  paragraphs  briefly  discuss  the  input 
data  obtained  for  airdrop  areas  1,  4,  and  5 for  use  with  equations  (Ap- 
pendix C)  that  predict  penetration  depths  of  air-implanted  SID's, 

134.  Average  initial  Young's  modulus.  The  average  initial 

Young's  modulus  E was  obtained  by  taking  an  average  E value  from  the 

deviator  stress-axial  strain  relations,  such  as  those  illustrated  by  the 

example  in  fig.  10.  The  E values  for  airdrop  areas  1,  4,  and  > are 

given  in  table  13-  The  dry-season  E values  varied  from  560  to 
2 

1400  kg/cm  with  test  pit  5 in  airdrop  area  1 naving  the  lowest  value 
and  test  pit  6 in  airdrop  area  1 and  test  pit  2 in  airdrop  area  5 having 
the  highest  values.  The  wet-season  E values  are  considerably  lower, 
especially  for  airdrop  areas  1 and  5- 

135,  Ultimate  shear  strength.  The  ultimate  shear  strength 

was  obtained  from  the  relations  -of  maximum  stress  difference  (o  - a ) 

a r 

versus  mean  normal  stress  (o  + 2o  )/3  presented  in  fig.  37 • The  value 

cl  T 

or  ou  was  determined  by  the  procedure  given  in  paragraph  yi.  for  tne 
dynamic  yield  envelopes  (fig-  37)  corresponding  to  wet-season  conditions 
(i.e.  high  moisture  values)  for  all  airdrop  areas  and  dry-season  condi- 
tions (i.e.  low  moisture  values)  for  airdrop  area  5,  pit  2,  the  ultimate 
shear  strengths  were  determined  from  the  maximum  obtainable  values  of 
o_  — 0 _ portrayed  in  fig.  37?  utilizing  equation  2.  For  the  dynamic 
yield  envelopes  corresponding  to  dry-season  conditions  for  airdrop 
area  1,  pit  6,  airdrop  area  4,  pits  3 and  4,  and  airdrop  area  b,  pit  3 , 
a different  procedure  was  used.  The  dynamic  yield  envelopes  portrayed 
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by  these  airdrop  areas  did  not  reach  a maximum  value  of 


(i  .e . 

the  slopes  of  the  curves  did  not  approach  zero)  for  the  range  of  mean 
normal  stresses  considered  in  the  triaxial  test.  Since  additional  high 
pressure  (i.e.  mean  normal  stress)  stress-strain  data  would  have  been 
required  to  obtain  a maximum  value  of  o - o for  these  soil  materials, 

£L  37 

and  since  the  system  of  equations  used  to  predict  the  depth  of  penetra- 
tion of  the  ADSID/S  was  determined  to  be  very  insensitive  to  ultimate 

2 

shear  strengths  larger  than  approximately  85  kg/cm  for  the  velocity 
range  considered  in  this  study,  it  was  decided  that  oo  - a,,  values 

a * 

corresponding  to  a mean  normal  stress  of  about  300  kg/cm  would  be  suf- 
ficient for  use  in  computing  the  ultimate  shear  strength  of  each  of 
these  airdrop  areas  (i.e.  airdrop  area  1,  pit  6,  airdrop  area  4,  pil  s 3 
and  4,  and  airdrop  area  5,  pic  1?  during  dry-season  conditions).  The 
o values  used  to  determine  penetration  depths  for  the  various  airdrop 
areas  are  given  in  table  13. 

136.  The  largest  value  of  was  obtained  at  airdrop  area  4 and 

the  lowest  aL  airdrop  area  1 during  both  the  wet  and  dry  seasons.  In- 
crease in  soil  moisture  reduced  the  value  of  significantly.  For 

the  three  airdrop  areas  for  which  ultimate  shear  strength  data  were  ob- 
tained, the  wet-seuson  to  dry-season  variations  were  greater  than  the 
variations  among  airdrop  areas  1,  4,  and  5* 

137.  Compressibility  parameter ■ The  compressibility  parameter 

p / p was  computed  using  the  uniaxial  strain  data  summarized  in  fig.  38 
3? 


and  the  following  field-measured  values  for 
cedure  is  given  in  paragraph  90. 


The  computation  pro- 


initial  Soil  Density 


(p) , s/cm3 


Drop 

Test 

Dry 

Wet 

Area 

Pit 

Season 

Se as on 

1 

5 

1.99 

2.14 

6 

1.73 

1.84 

4 

3 

1 .61 

1-73 

4 

1.63 

1.7b 

5 

1.60 

1.70 

2 

1.70 

1.81 
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The  computed  values  in  table  13  for  the  dry  season  vary  from  1.18  to 
1.33  and  follow  the  same  ranking  trend  as  the  ultimate  shear  strength 
data,  in  that  the  lowest  values  are  for  airdrop  areas  1 and  ‘3  and  the 
highest  values  are  for  airdrop  area  4.  The  season-to-season  variation 
in  the  compressibility  parameter  data  is  smallest  for  airdrop  area.  5 
and  greatest  for  airdrop  area  1. 

133.  Drag  coefficient.  The  drag  coefficient  computed  as 

discussed  in  paragraph  92  resulted  in  a value  of  0.2  for  all  sites  and 
conditions,  as  shown  in  table  13. 

139.  Supplementary  information.  In  addition  to  the  information 
discussed  in  the  previous  paragraphs,  soil  profiles  for  airdrop  areas  1 
and  4-6  were  prepared  (fig.  39)-  Also  some  additional  descriptive  in- 
formation taken  in  support  of  the  seismic  response  tests  is  presented 
in  table  2 for  airdrop  areas  1 and  3-6. 
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F rediction  and  Evaluation  of  MIMISID  Detection  Performance 
SID  detection  performance  model 

l40.  A theoretical  simulation  of  the  phenomena  involved  in  sensor 
detection  performance  as  depicted  in  fig.  1 is  a complex  mathematical 
problem,  and  a practical  solution  will  require  extensive  use  of  the  com- 
puter. As  stated  in  paragraph  39,  the  detection  problem  has  been  di- 
vided into  four  component  parts  as  follows : 

el.  The  development  of  a forcing  function  that  simulates 

transfer  of  energy  from  the  seismic  signal  source  to  the 
substrate  material. 

b_.  The  development  of  a seismic  wave  propagation  model  to 
simulate  the  propagation  of  microseismic  waves  within 
the  substrate  materials. 

c_.  The  development  of  functions  that  simulate  the  transfer 
of  energy  from  rhe  vibrating  medium  (substrate)  to  the 
sensor . 

d_.  The  development  of  functions  that  simulate  the  logic 
response  of  the  sensor  to  the  vibrating  medium. 
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141.  Significant  progi'ess  has  been  made  to  date  on  each  component 
part  of  the  overall  sensor  detection  model.  It  is  felt  that  once  two 
serious  omissions  are  incorporated  into  parts  a_  and  b_  above,  the  perfor- 
mance (i.e.  detection  range)  of  a SID  system  for  both  man-walking  and 
moving- vehicle  targets  can  be  predicted  realistically  and  quantitatively 
by  utilizing  the  characteristics  i f the  energy  source  and  the  input  ter- 
rain surface  and  subsurface  factor  data.  The  omissions  in  parts  a.  and 
b_,  respectively,  are  : 

a.  Component  part  a.  A simulation  of  the  stress  applied  to 
the  ground  from  the  source  (i.e.  a realistic  force-time 
history  of  the  source  moving  over  the  ground)  has  not 
been  modeled. 

b_.  Component  part  d . Macrogeometry  coefficients  represent- 
ing local  site  conditions  have  not  been  available. 

142.  Functions  describing  the  force-time  histories  of  the  energy 
sources  are  presently  being  developed  by  WES  and  will,  hopefully,  ac- 
count for  the  transfer  of  energy  from  specified  sources  (i.e.  drop  ham- 
mer, footstep,  vehicle,  etc-.)  into  the  propagating  substrate  material. 

In  support  of  b_  above,  a literature  review  has  been  completed,  and  con- 
siderable information  on  the  effects  of  macrogeometry  on  Rayleigh  wave 
propagation  has  been  compiled.  However,  these  data  have  not  been  veri- 
fied, and  at  the  time  of  this  report,  they  had  not  been  incorporated 
into  the  overall  detection  prediction  scheme. 

SID  detection  distance 
predictions  and  evaluations 

1 4 3 • Forcing  function  and  seismic  wave  propagation  submodel  pre- 
dictions , To  demonstrate  the  applicability  of  the  forcing  function  and 
seismic  wave  submodels  (i.e.  comp>onent  parts  a_  and  b_  of  the  overall  de- 
tection model  presented  in  paragraph  l40),  predictions  of  the  analog 
signal  as  a function  of  distance  from  the  source  were  derived  for  the 
drop-hammer  energy  source.  These  predictions  are  discussed  briefly  in 
the  following  paragraphs. 

144 . Laboratory  measurements  indicate  that  the  restoring  force 
associated  with  the  drop  hammer-soil  system  is  a nonlinear  function  of 
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the  displacement  of  the  drop-hammer  weight.  The  time-dependent  forces 
exerted  on  the  ground  we*  calculated  by  using  the  Rungo-Kutta  method  to 
numerically  solve  the  differential  equation  of  motion  of  the  drop-hammer 
weight  subject  to  a nonlinear  restoring  force.  The  calculations  indi- 
cate that  the  forcing  function  can  be  approximated  by  a cosine  squared 
pulse  whose  amplitude  and  pulse  width  are  controlled  by  site  conditions. 
To  illustrate  the  effects  of  using  a cosine  squared  pulse  to  simulate 
the  force-time  history  of  the  drop-hammer  energy  source,  it  was  experi- 
mentally determined  (see  reference  16)  that  a pulse  width  of  0.016  sec 
and  a peak  amplitude  of  71.7  * 10^  dj nes  were  appropriate  for  the  condi- 
tions for  a selected  site  (site  3*4l ) at  the  time  of  test.  (The  cosine 
squared  pulse  obtained  for  site  5**1  is  approximately  equivalent  to  the 
stress  produced  by  a footstep  of  a man  weighing  73  kg  and  wearing  size  9 
boots.)  This  information  was  then  used  in  conjunction  with  the  field- 
measured  properties  listed  below  to  obtain  estimates  of  the  span,  peak 
amplitude,  and  rms  amplitude  within  the  predicted  analog  signal. 


j a.  P-wave  velocity  of  first  refraction  layer  = 365  m/sec. 

I b.  S-wave  velocity  of  first  refraction  layer  = 20*4  m/sec. 

3 

c_.  Wet  density  of  first  refraction  layer  - 1.75  g/cm  . 
d_.  Thickness  of  first  refraction  layer  = 1.15  m. 
c_.  ?-wave  velocity  of  infinite  half-space  (.i.e.  no  second 
layer  was  encountered)  = 8l8  m/sec. 

f_.  S-wave  velocity  of  infinite  half-space  = 333  m/sec. 

3 

Wet  density  of  infinite  half-space  = 1.9  g/cm  . 
i h_.  Soil  damping  coefficient  c — 0,03  . 

l*+5.  The  resulting  predictions  of  span,  peak,  amplitude,  and  rms 
amplitude  are  plotted  along  with  the  measured  values  in  fig,  *j0.  As  can 
be  seen,  the  measured  and  predicted  signals  compare  well  in  amplitude , 
frequency  content,  and  duration.  The  frequency  content  of  the  cosine 
squared  pulse  appears  to  simulate  the  force-time  history  of  the  drop- 


hammer  energy  source , as  indicated  by  the  close  agreement  between  the 
predicted  and  measured  rms  amplitudes. 

l*+6.  The  results  presented  (fig.  *40 ) are  for  one  site  only,  but 
similar  results  have  been  derived  for  other  sites  using  similarly 
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derived,  forcing  functions.  In  addition,  some  predictions  using  forcing 
functions  derived  for  the  M151  and  M113  'vehicles  have  been  made  and  com- 
pared with  measured  values.  From  the  results  obtained  to  date  at  WES, 
it  seems  reasonable  to  assume  that  current  research  will  result  in  use- 
ful and  workable  solutions  to  the  problems  identified  as  component 
parts  a and  b_  in  paragraph  l40. 

147-  SID  log: c modeling.  To  study  the  problems  associated  with 
simulation  of  the  transfer  of  energy  from  the  vibrating  medium  to  the 
SID  and  the  logic  response  of  the  SID  (component  parts  c_  and  A_  in  para- 
graph l40),  use  was  made  of  the  Phase  111  MINISID  and  analog  man -walking 
test  results.  To  derive  a medium-gain  detection  distance,  the  field- 
measured  analog  recordings  of  the  seismic  response  from  a man-valking 
test  were  analyzed  by  an  analog  computer  programmed  to  duplicate  (model) 
the  functioning  logic  of  the  MINISID.  (A  description  of  the  analog  com- 
puter model  is  given  in  Appendix  D. ) 

148.  Measured  MINISID  detection  distances,  the  detection  dis- 
tances predicted  by  the  analog  computer  model,  and  the  differences  be- 
tween the  two  are  presented  in  table  14.  A plus  Sign  (+)  indicates  that 
the  model-predicted  distance  was  greater  than  the  measured  distance,  and 
a minus  sign  (-)  indicates  that  it  was  less.  Also  presented  are  the 
mean  difference,  standard  deviation,  and  the  number  of  sites  underpre- 
dicted and  ovei-predicted  for  each  geographical  area  (i.e.  Fort  Bragg, 
wet  and  dry  seasons;  Fort  Huachuca,  dry  season;  and  several  sites  in 
West  Germany). 

149-  From  table  l4 , it  f'an  be  seen  that  the  measured  and  pre- 
dicted values  from  the  Fort  Huachuca  test  sites  were  extremely  close, 
with  the  mean  and  standard  deviation  of  the  differences  being  3.1  and 
i+.U  m,  respectively.  Also,  little  bias  toward  overprediction  or  under- 
prediction was  observed,  since  eight  sites  were  overpredicted  and  seven 
were  underpredicted.  The  data  from  West  Germany  show  a strong  bias  for 
overprediction  and  a relatively  poor  comparison,  i.e.  mean  and  standard 
deviation  of  12  and  7-5,  respectively.  The  Fort  Bragg  results  were  in- 
termediate between  the  Fort  Huachuca  and  the  West  Germany  data. 

150.  The  measured  and  predicted  (modeled)  detection  distances 
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(table  14 ) for  sites  at  each  of  the  three  geographical  areas  are  quite 
close,  especially  when  the  following  items  are  considered.  The  pre- 
dicted distances  were  derived  from  one  test  (or,  in  some  cases,  two), 
whereas  the  MIN1S1D  detection  distances  were  averages  of  from  5 to 
15  tests.  Also,  the  MlNlSlD's  used  were  production  line  items  and  were 
not  the  same  in  each  area;  the  relative  performances  of  the  devices  are 
not  Known.  Further  the  background  noise  from  one  test  to  another  could 
not  he  controlled,  and  as  the  wind  or  cultural  activity  changed,  the 
background  noise  changed  from  test  to  test,  and  this  induced  significant 
effects  on  the  test  results  (see  paragraph  110 ). 

151.  Summary . The  implication  of  these  tentative  analog  computer 
results  is  that  component  parts  c_  and  q_  (paragraph  lho)  can  be  solved 
analytically.  The  transfer  functions  for  any  particular  sensor  can  be 
readily  derived  and  programmed  on  an  analog  or  digital  computer,  thus 
giving  a means  for  predicting  SID  performance  if  the  analog  wave  train 
is  available.  The  effect  of  geophone  coupling  was  not  adequately  tested 
in  these  experiments  because  good  coupling  can  usually  be  accomplished 
with  the  external  geophone  of  the  M1NIS1D.  Further  study  is  needed  to 
determine  coupling  efficiency  as  a function  of  sensor  mass  and  geometry. 

Prediction  and  Evaluation  of  Penetration  Depths  and  Impact 
Deceleration  Values  for  the  ADSID/S 

152.  To  derive  an  estimate  of  the  range  of  air-implantation  dif- 
ficulty and  the  amount  of  deceleration  that  a typical  air-delivered  SID 
(Phase  III  ADSID/S,  MA3T ) would  experience  at  airdrop  areas  1,  4 , and  1 , 
a theoretical  study  was  performed  using  the  system  cf  equations  in  Ap- 
pendix C.  The  depth  of  penetration  for  vertical  (normal)  impact  and  for 
an  angle  of  U5  deg  (oblique  impact)  and  peak  deceleration  values  for 
normal-impact  conditions  were  predicted  and  are  discussed  in  the  follow- 
ing paragraphs.  For  both  the  normal  and  oblique  impact  configurations, 
the  penetration  calculations  were  terminated  when  the  fin  assembly  (see 
fig.  13)  of  the  ADSID/S  sensor  came  in  contact  with  the  surface  of  the 
soil.  This  was  done  since  no  information  was  available  on  the  behavior 
of  the  fin  assembly  and  the  tailplate  to  determine  whether  the  assembly 
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would  break,  bend,  or  shear  off  when  it  came  in  contact  with  the  sub- 
strate material  during  the  penetration  event. 

Prediction  of  depth  of  penetration 

l‘>3.  Normal -impact  predictions.  Theoretical  predictions  of  depth 
of  penetration  and  deceleration  were  obtained  for  a normal-impact  con- 
figuration of  the  ADSID/S  for  each  of  the  six  sites  (i.e.  six  soil  pits) 
within  airdrop  areas  1,  h,  and  5.  Two  moisture  contents  were  used  that 
were  assumed  to  be  representative  of  the  wet  and  dry  seasons,  respec- 
tively. The  predicted  depths  are  presented  as  a function  of  impact  ve- 
locity in  fig.  1*1.  These  relations  reflect  the  terrain  variations  en- 
countered within  each  of  the  areas.  The  location  of  the  test  pits 
represented  by  these  predictions  is  shown  in  fig.  16. 

15l*.  Oblique-impact  predictions.  Theoretical  predictions  were 
also  made  for  oblique  (off-normal)  impact  configurations  for  each  of  the 
six  sites  within  airdrop  areas  1,  U,  and  5,  again  with  two  moisture  con- 
tents used.  A theoretical  solution  for  predicting  the  penetration  of  a 
projectile  striking  the  soil  at  an  oblique  angle  has  not  been  adequately 
formulated,  however,  experience  at  ml 3 has  shown  that  for  shallow  pene- 
trations a reasonabxe  approximation  of  the  depth  of  penetration  for 
oblique  impacts  can  be  obtained  by  resolving  the  velocity  vector  (i.e. 
velocity  along  the  trajectory  of  the  projectile)  into  its  horizontal  and 
vertical  components  and  then  determining  the  vertical  penetration  (i.e, 
the  vertical  distance  from  the  surface  of  the  ground  to  the  tip  of  the 
nose  of  the  projectile)  corresponding  to  the  vertical  component  of  the 
velocity  vector.  In  order  to  obtain  the  penetration  path  length,  the 
vertical,  penetration  depth  must  be  divided  by  cos  0 , where  0 is  the 
angle  of  impact  off  vertical.  Mathematically,  the  vertical  (normal)  im- 
pact velocity  is  expressed  as 


V = V cos  0 
v a 


= normal  component  of  impact  velocity  vector 


where 


V = velocity  vector  along  longitudinal  axis  of  the  ADS  uD/S 
(i.e.  along  the  projectile  trajectory) 

6 = angJe  of  impact  (off  vertical) 

155.  This  technique  was  used  to  compute  the  d ;pth-of-penetration 
predictions  for  an  impact  angle  of  45  deg  (fig.  42).  The  impact  ve- 
locity portrayed  in  fig.  42  is  the  velocity  vector  along  the  45-deg 
trajectory. 

156.  Prom  the  data  (fig.  4l)  relating  depth  of  penetration  to  im- 
pact velocity  for  vertically  impacting  ADSID/S's,  it  is  noted  that  there 
is  considerable  variation  in  the  penetration  depths  within  and  between 
the  drop  areas,  and  that  the  penetration  depths  are  very  dependent  upon 
the  amount  of  soil  moisture.  To  illustrate  the  effects  of  moisture  on 
the  prediction  depths,  let  it  be  assumed  that  20  cm  cf  penetration  are 
required  to  ensure  that  the  ADSID/S  is  emplaced  sufficiently  deep  to  op- 
erate as  designed.  From  fig.  4l,  the  impact  velocities  required  to  ob- 
tain a 20 -cm  depth  of  penetration  for  each  airdrop  area  are  as  follows  : 


Impact  Velocity 
(m/sec)  for  20-cm 
Penetration 

Area 

Pit 

Wet 

Dry 

1 

5 

25 

38 

6 

50 

>100 

4 

3 

78 

>100 

4 

89 

>100 

5 

1 

68 

>100 

2 

48 

>100 

157-  The  wet-season  impact  velocities  for  a depth  of  penetration 
of  20  cm  vary  from  25  m/sec  (area  1,  pit  5 ) to  89  m/sec  (area  4,  pit  4); 
whereas  the  dry-season  velocities  vary  from  38  m/sec  (area  1,  pit  5)  to 
>100  m/sec  (all  other  five  pits).  The  data  tabulated  above  suggest  that 
during  the  wet  season,  area  1 offers  the  least  penetration  resistance; 
whereas  during  the  dry  season  only  that  portion  of  area  1 represented  by 
pit  5 offers  the  least  resistance.  During  the  dry  season,  area  1,  pit  6, 
and  areas  4 and  5 require  impact  velocities  greater  than  100  m/sec  to 
obtain  a 20-cm  penetration.  The  impact  velocity  normally  used  (i.e.  the 
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anticipated  velocity)  for  air-implantation  of  the  ADSID/S  is  approxi- 
mately 50  m/sec.  The  predicted  depth-of-penetration  curves  show  that  at 
an  anticipated  optimum  impact  velocity  of  50  m/sec,  a depth  of  penetra- 
tion of  20  m would  occur  only  within  area  1,  pit  5,  in  loth  vet  and  dry 

seasons,  and  within  area  1,  pit  6,  and  area  5,  pit  2,  during  the  wet 
season. 

158.  In  general,  the  ADSID/S  impacts  the  ground  at  some  angle  off 
the  vertical.  Therefore,  the  results  shown  in  fig.  4l  show  the  highest 

value  of  penetration  that  could  be  expected  at  each  drop  area.  Pig.  42 

shows  similar  results  for  impact  angles  at  45  deg  off  the  vertical. 
Normally,  the  impact  angles  are  less  than  45  deg  off  the  vertical  (see 
table  12  for  experimental  data);  therefore,  these  curves  should  repre- 
sent the  minimum  amount  of  penetration  expected.  As  can  be  seen  from 
the  predicted  curves,  the  impact  velocity  to  allow  a penetration  of 
20  cm  would  be  gx'eater  than  100  m./sec  for  all  sites  (except  area  1, 
pit  5)  during  the  dry  season  and  less  than  100  m/sec  for  area  1 and 
pit  2 of  area  5 during  the  wet  season.  If  the  impact  velocity  is  as- 

{ sumed  to  be  50  m/sec,  a 20-cm  penetration  path  length  would  occur  only 

. at  airdrop  area  .1,  pit  5. 

! Prediction  of  maximum  deceleration 

j 159-  The  maximum  rigid-body  decelerations  calculated  for  impact 

j velocities  of  50  m/sec  for  airdrop  areas  1,  4,  and  5 are  listed  in 

table  15.  These  calculations  are  for  normal  impacts  of  the  ADSID/S- 
The  specified  maximum  allowable  deceleration  for  a workable  Phase  III 
ADIGD/S,  MA37,  is  approximately  9*85  x 1C^  m/sec^.  All  computed  deeel- 

j 

; eration  values  (see  table  15 ) were  over  the  maximum  allowable  decelera- 

tion value  for  a workable  (i.e.  operational)  ADSID/S,  except  those  for 
area  1,  pit  5,  during  the  wet  and  dry  seasons,  and  area  1,  pit  6,  and 

area  5,  pit  2,  during  the  wet  season. 

Evaluation  of  theoretical 
1 depth-of-penetration  pre- 

dietjons  for  oblique  impacts 

l60.  Since,  as  stated  in  paragraph  154,  an  adequate  theory  for 
predicting  penetration  of  obliquely  impacting  projectiles  is  not 
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available,  it  was  desired  that  the  theoretical  depth-of-penetration  pre- 
dictions made  herein  be  experimentally  verified.  The  ADSID/S  tests  con- 
ducted by  CSEWS  at  Fort  Huaehuea  provided  a limited  amount  of  data  that 
could  be  used  to  estimate  how  well  the  theoretically  predicted  values, 
for  obliquely  impacting  projectiles,  matched  the  measured  ones.  To  de- 
termine the  terrain  factor  data  in  the  CSEWS  airdrop  area,  i.e.  E , 


, Pp/P  > V+  » and  Cp  , laboratory  tests  were  conducted  on  one  bulk 


surface  soil  sample  obtained  from  the  area.  Because  of  the  small  amount 
of  soil  available  on  the  CSEWS  area,  laboratory  tests  (uniaxial  and  tri- 
axial)  were  conducted  at  only  one  moisture  content,  i.e.  10  percent. 
These  results,  along  with  the  results  from  two  standard  laboratory 
tests,  grain-size  distribution  and  Atterberg  limits,  were  compared  with 
corresponding  data  from  all  the  sites  in  drop  areas  1,  4,  and  5-  It  was 
determined  that  the  terrain  conditions  at  the  CSEWS  airdrop  area  were 
similar  to  terrain  conditions  at  area  1,  pit  5,  during  the  wet  season 
(i.e.  l4  percent  moisture).  This  can  be  seen  in  fig.  1+3,  where  gro.in- 
size  uistributions  for  the  two  sites  correspond  quite  closely,  and  in 
fig.  41+ , where  the  plasticity  index  and  liquid  limit  data  correspond 
quite  closely.  More  significant  is  the  similarity  of  the  dynamic  yield 
surfaces  from  the  triaxial  compression  tests  for  the  two  sites  as  shown 
in  fig.  45, 

l6l . Based  on  the  similaritv  between  the  test  site  conditions  and 


the  results  of  the  triaxial  compression  tests,  it  was  assumed  that  the 
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area  1,  pit  5 (fig-  4l),  could  be  extrapolated  to  the  CSEWS  area,  once 
they  were  corrected  for  moisture  conditions  at  the  time  of  test.  From 
fig.  1+5  it  can  be  seen  that  the  dynamic  yield  surfaces  derived  from  the 
triaxial  tests  were  nearly  identical  for  the  soils  at  the  CSEWS  airdrop 
area  at  10  percent  moisture  content  and  at  airdrop  area.  1,  pit  5 ? at 
10-ll+  percent  moisture  content.  For  this  reason  it  was  assumed  that 
penetration-impact  velocity  relations  generated  at  l4  percent  moisture 
content  for  area  1,  pit  5,  would  correspond  to  the  penetration-impact 
velocity  relations  at  the  CSEWS  airdrop  area  at  10  percent  moisture 


content.  The  relation  for  the  CSEWS  area  is  shown  as  a solid  line  in 
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(pit  5)  and  CSEWS  airdrop  area 
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Fig.  46.  Prediction  of  normal-impact  penetration  curves  for  various 
moisture  contents , CSEWS  airdrop  area 

fig.  46.  Table  l6  shows  that  the  experimental  ADSID/S  tests  (where  im- 
pact velocities  were  measured)  were  conducted  in  soils  whose  moisture 
values  ranged  from  5-5  to  8.6  percent.  To  provide  a basis  for  estimat- 
ing the  depth  of  penetration  at  these  moisture  contents , the 
penetration-impact  velocity  data  shown  as  the  solid  line  in  fig.  46  were 
extrapolated  to  cover  moisture  contents  at  the  time  of  the  experimental 
ADSID/S  tests.  These  extrapolated  penetration  impact  relations  are 
shown  as  dashed  lines  in  fig.  46.  The  extrapolation  was  linear  based  on 
the  spread  of  the  10  and  14  percent  moisture  curves  shown  in  fig.  4l  for 
area  1,  pit  5- 

162.  Near-surface  impact  velocities  were  measured  for  11  tests 
(i.e.  tests  5,  9,  11,  12,  34,  36,  75,  76,  and  78-80)  at  the  CSEWS  air- 
drop area.  These  tests  were  used  as  a basis  for  comparing  the  predicted 
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penetration  values.  The  test  data  obtained  at  the  CSEWS  area  include. 

(a)  impact  velocity  vector  along  the  trajectory  of  the  ADSiD/S,  (b) 
ADSID/S  impact  angle,  and  (c)  surface  soil  moisture  content  and  strength. 
The  predicted  depth  of  penetration  derived  from  the  graphs  shown  in 
fig.  46  were  compared  with  the  field-measured  depth  of  penetration  val- 
ues (table  16).  Because  all  the  test  sensors  (i.e.  ADSID/S)  impacted 
the  ground  at  an  angle  off  the  vertical,  the  predicted  length  of  pene- 
tration (as  opposed  to  vertical  depth  of  penetration)  is  used  so  that 
the  predicted  and  measured  values  will  be  more  directly  comparable. 

163.  The  method  used  to  predict  the  depth  of  penetration  was 
identical  with  that  described  in  paragraph  154  and.  involved  determining 
the  vertical  impact  velocity  (see  bottom  of  table  16 ) from  the  impact 
velocity  vector.  Next,  this  velocity  was  used  to  enter  the  curves  in 
fig.  46  at  the  appropriate  field  moisture  content  to  determine  the  cor- 
responding depth  of  penetration  from  vhi'.n  the  length  of  penetration  was 
calculated.  This  method  of  obtaining  a prediction  of  the  length  of 
penetration  does  not  adequately  simulate  the  penetration  phenomena  tor 
oblique  impacts , but  experience  has  shown  that  it  yields  relatively  good 
approximation  of  penetration  lengths , especially  for  shallow 
penetrations . 

l6U.  The  results  of  this  study  appear  to  confirm  this  also  be- 
cause it  can  be  seen  that  the  deviations  from  measured  penetration 

lengths  range  from  an  underprediction  of  12  cm  (test  75,  table  16)  to 
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analytical  procedures  used  to  derive  the  penetration  estimates  are  con- 
sidered, the  predicted  results  are  surprisingly  compatible  with  the 
measured  values.  Although  this  comparison  cannot  be  interpreted  as  a 
verification  of  the  theoretical  depth-of-penetration  predictions 
(figs.  4l  and  42)  it  does  increase  confidence  in  the  predicted  penetra- 
tion depths  of  the  air-implanted  SID's. 

Summary 

.165.  The  airdrop  areas  (areas  1-6  and  the  CSEWS  area)  provided 
terrains  covering  a wide  range  of  penetration  characteristics.  This 
variation  in  terrain  conditions  allows  the  nesting  oi  air-implanted 


105 


I 


SID's  under  varying  stress  conditions  (i.e,  high  to  low  rigid-body  de- 
celeration values)  and  under  a variety  of  soil  strength  conditions. 


i 
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Extrapolation  of  SID  Performance  Vaises 


166.  The  WES  concept  for  the  extrapolation  of  SID  test  results  is 
based  on  the  assumption  that  adequate  simulation  models  describing  the 
phenomena  under  test  can  be  formulated.  In  the  case  of  SID  detection 
performance,  this  involves  simulating  the  generation,  propagation,  and 
coupling  of  the  seismic  energy  from  the  source  to  the  sensor.  In  the 
case  of  air-implantation  of  sensors,  it  involves  simulating  the  penetra- 
tion event  of  an  air-delivered  SID  and  the  coupling  of  the  SID  to  the 
soil.  The  extrapolation  concept  also  requires  that  the  terrain  inputs 
to  the  simulation  models  be  constitutive  arid  measurable  properties  of 
the  terrain,  such  that  they  can  be  classified  and  portrayed  on  an  areal 
basis.  Based  on  the  data  presented  in  this  part  of  the  report,  it  is 
felt  that  adequate  simulation  models  required  to  implement  the  concept 
can  he  derived.  For  a discussion  of  the  mapping  of  terrain  factors  on 
an  areal  basis  see  Part  IV  of  this  report. 


PART  IV:  MAPPING  OF  TERRAIN  AND  SEISMIC  RESPONSE 

FACTORS  AND  SID  PERFORMANCE  VALUES 


167-  The  theoretical  detection  prediction  scheme  and/or  concepts 
presented  in  this  report  provide  an  analytical  framework  for  obtaining 
SID  performance  predictions  for  man-walking  targets  in  various  world  ter- 
rains. A point  worthy  of  reemphasizing  is  that  the  prediction  scheme  re- 
quires specific  terrain  and  seismic  input  data  for  execution.  In  pursuit 
of  a general  SID  detection  prediction  capability,  ais  illustrated  in 
fig.  2,  this  portion  of  the  report  pertains  to  the  development  of  proced- 
ures for  preparing  terrain  factor  maps  that  will  delineate  uhe  geograph- 
ical distribution  of  the  individual  factors  describing  the  relative  seis- 
mic response  of  the  area  and  those  factors  required  or  hypothesized  as 
being  required  as  input  to  the  SID  detection  performance  model  (i.e. 
forcing  function  and  seismic  wave  propagation  submodels). 

168.  The  terrain  and  seismic  factor  mapping  discussed  herein  per- 

tains to  the  Fort  lluachuca  hand-emplacement  test  area  and  consists  of 
preparing  four  types  of  maps:  (a)  a terrain  factor  complex  map  for  sub- 

surface layers  1 and  2 to  provide  inputs  to  the  theoretical  seismic  wave 
propagation  model,  (b)  surface  condition  maps  to  be  used  witn  the  forcing 
function  model  under  development,  (c)  seismic  response  factor  complex 
maps  to  show  the  relative  seismic  response  of  the  hand-emplacement  area, 
and  (d)  SID  performance  maps.  The  terrain  and  seismic  response  factor 
complex  maps  and  the  SID  performance  maps  are  based  on  data  obtained  by 
WES  field  teams  during  August  1971  and  March  1972  and  are  assumed  to  he 
representative  of  the  wet-  and  dry-season  conditions  at  the  test  sites. 

169.  The  terrain  and  seismic  factors  and  the  SID  performance  val- 
ues (i.e.  detection  distances)  mapped  for  the  hand-emplacement  area  at 
Fort  Huachuca  provide  a direct  basis  for  comparing  (i.e.  correlating)  the 
terrain  and/or  seismic  factor  conditions  in  other  geographic  regions,  in- 
cluding the  test  areas  at  Fort  Bragg'*'1*’1'  and  in  the  Panama  Canal 
Zone.10,1^  The  following  paragraphs  describe  the  detailed  procedures  used 
in  preparing  the  terrain  and  seismic  factor  maps  and  the  SID  performance 
map  for  the  Fort  Huachuca  hand-emplacement  area. 
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Terrain  Factor  Mapping  for  Subsurface  Soil  Layers 


170.  The  four  terrain  factors  that  have  been  mapped  are : 

l 

a_.  Thickness  of  layer  1 (T  ) and  layer  2 (T0). 

b_.  Compression  wave  velocity  for  lsyer  1 (Vp^)  and  layer  2 

(VP2 ^ * 

c.  Shear  wave  velocity  for  3.ayer  1 (V  ) and  layer  2 (V  ) . 
d_.  Wet  density  for  layer  1 (p  ) and  layer  2 (p^). 

Terrain  factor  mapping  classes 

171.  Since  classification  is  a prerequisite  to  any  mapping  pro- 
cess , and  the  adequacy  of  uhe  proposed  map  depends  upon  its  class  inter- 
vals , considerable  effort  was  devoted  to  determining  appropriate  terrain 
mapping  classes.  Consideration  also  was  given  to  selecting  classes  that 
could  be  used  in  comparing  terrain  conditions  e.t  Fort  Ruachuea  with 
other  TECOM  SID  test  areas.  The  class  intervals  or  class  ranges  of  each 
of  the  terrain  factors  were  chosen  to  minimize  the  effect  of  stratifica- 
tion. This  was  accomplished  by  selecting  class  ranges  than  would  mini- 
mize the  error  when  midpoint  values  on  which  to  make  predictions  of  SID 
detection  performance  are  used;  however,  class  intervals  too  small  to  be 
mapped  adequately  were  avoided. 

172.  Two  sets  of  terrain  factor  classes  were  established  prior  to 
mapping,  one  set  for  layer  3.  and  the  othet'  for  layer  2.  Before  any  map- 
ping was  accomplished,  the  midpoint  values  of  the  tentatively  selected 
class  ranges  were  input  to  the  theoretical  seismic  wave  propagation 
model  to  detennine  the  adequacy  of  the  selected  ranges . If  the  PPV  and 
rms  amplitudes  of  the  predicted  signal  obtained  using  the  theoretical, 
model  compared  closely  with  the  PPV  and  rms  amplitudes  resulting  from 
the  field-measured  drop-hammer  signal,  the  mapping  class  ranges  were 
considered  adequate.  If  the  comparison  showed  differences  greater  than 
15  p>ereent  (for  distances  greater  than  20  m),  smaller  mapping  class 
ranges  were  selected,  and  the  above-mentioned  process  was  repeated  until 
adequate  classes  were  established  for  all  terrain  factors  (Vp,Vg,T.  and 
p).  The  four  terrain  factors  and  their  selected  mapping  classes  are 
listed  in  table  17. 
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Mapping  of  terrain 
factors  for  layers  1 and  2 


173-  The  terrain  factors  for  layers  1 and  2 of  the  test  area  were 
mapped  by  using  the  environmental  data  collected  during  the  field  pro- 
gram, air-photo  interpretation  techniques,  a topographic  map  cf  the  area, 
and  the  interpreter's  personal  knowledge  of  the  mapped  area, 

IT1*  • The  air-photo  interpretation  was  accomplished  through  a 
series  of  steps,  the  first  being  to  obtain  the  necessary  overlapping 
(stereo)  photos  of  the  hand-emplacement  area.  This  was  followed  by  a 
stereoscopic  examination  of  the  overlapping  prints,  whereby  various  pho- 
tographic terrain  patterns  were  identified,  according  to  tone  and  tex- 
ture, together  with  their  topographic  relations . For  the  photographic 
patterns  within  which  ground  truth  (i.e.  site)  data  were  available,  the 
class  range  as  portrayed  by  the  data  for  that  particular  site  was  used. 

In  those  patterns  without  ground  truth  data,  the  class  ranges  for  eacn 
terrain  factor  were  assigned  by  extrapolation  from  the  points  of  known 
ground  truth  data,  through  associations  of  similar  patterns  and  through 
the  interpreter's  knowledge  of  the  area.  After  all  the  identified  pat- 
terns ha^been  outlined  on  the  air  photos,  a map  unit  symbol  represent- 
ing a factor  value  class  was  assigned  to  each  respective  pattern.  In 
effect,  the  resulting  map  portrayed  the  four  factor  value  classes 
(table  17)  for  one  terrain  factor,  such  as  wet  density  (p).  Similar 
maps  were  then  prepared  for  V , V0  , and  T , respectively. 

X O 

175.  After  the  four  terrain  factox-s  (V  V„,T,  and  p)  for  each 

l iD 

of  the  two  layers  had  been  mapped , they  were  combined  into  a factor 
complex  map  (fig.  47)  by  superimposing  each  of  the  individual  maps  on  a 
common  base  map.  To  simplify  the  identification  of  each  delineated  area 
on  the  factor  complex  map,  the  areas  were  identified  by  numbers  from  1 
to  40.  These  numbers  represent  discrete  combinations  of  classes  of 
each  of  the  four  terrain  factors  for  the  two  soil  layers  in  the  hand- 
emplacement  area,  and  ai’e  contained  in  the  factor  complex  map  legend 
shewn  in  table  18. 

Evaluation  of  mapping  classes 

176.  To  determine  the  adequacy  of  the  selected  factor  map  class 
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ranges,  a check  was  made  hy  cumparing  the  theoretical  predicted  results 
(using  20-rn  measured  seismic  (analog)  signals  and  the  prediction  models 
described  in  Appendix  D)  for  sites  335-342,  348,  and  349  with  the  pre- 
dicted values  obtained  using  midpoint  values  of  the  class  ranges  from 
map  units  1,  3,  10,  45,  51,  64,  89,  96,  103,  and  115,  respectively.  Tne 
factor  values  used  in  the  model  to  obtain  the  predictions  for  the  site 
and  mapped  conditions  are  given  in  table  19.  The  results  of  this  com- 
parison are  shown  in  table  SO.  Excellent  correlation  of  results  was  ob- 
tained with  the  20-m  reference  signals,  the  largest  difference  being  2 m 
observed  at  site  337,  and  map  unit  103.  In  general,  these  results  indi- 
cate that  the  selection  of  the  factor  mapping  classes  was  a good  compro- 
mise between  prediction  accuracy  and  mapping  practicality. 

Mapping  of  Surface  Terrain  Factors 

177«  The  relative  effect,  of  surface  terrain  factors  on  the  trans- 
fer of  seismic  energy  from  the  source  to  the  substrate  has  not  been  quan- 
titatively determined.  Nevertheless,  surface  moisture  content,  soil 
type,  and  surface  CPR  almost  surely  control  the  transfer  of  seismic 
energy  for  man-walking  energy  sources.  For  this  reason,  these  terrain 
factors  will  almost  certainly  be  required  as  inputs  to  the  forcing  func- 
tion model  under  development.  These  three  surface  terrain  factors  have 
been  mapped  for  the  hand-emplacement  area,  and  are  discussed  in  the  fol- 
lowing paragraphs . 

Surface  soil  mo i sture  content 

178.  Six  factor  classes  of  sciu  moisture  content  for  the  dry  and 
wet  seasons  were  established  for  the  mapping  of  the  hand-emplacement 
area,  as  follows: 

Surface  Soil  Moisture 
Factor  Factor  Value  Range 

Class  percent 

1~3 
4-8 
9-12 


1 

2 

3 


(Continued) 


Factor 

Class 


Surface  Soil  Moisture 
Factor  Value  Range 
percent 

U 13-16 

5 17-20 

6 21-1*0 

Soil  moisture  classes  l*-6  did  not  occur  in  the  hand-emplacement  area 
during  either  the  wet  or  dry  season,  but  were  established  to  provide 
the  class  ranges  necessary  for  comparison  of  the  Fort  Huachuca  SID 
hand-emplacement  area  with  other  TECOM  SID  test  areas.  The  wet-  and 
dry-season  surface  soil  moisture  content  factor  map  for  the  SID  hand- 
emplacement  area  is  shown  in  fig.  1*8.  The  units  on  this  map  represent 
factor  classes  of  each  moisture  factor  class  (1-3)  range  for  the  dry 
and.  wet  seasons. 


Surface  soil  type 
(or  classification) 

179.  Five  factor  classes  of  surface  soil  type  (i.e.  soil  classi- 
fication by  USCS)  were  established  as  follows: 


Surface  Soil  Type 
Factor  Factor  Value  Range 

Class  USCS  Classification 


1 

2 

3 

1* 

5 


SP 

SC 

0L 

SM 

GM 


The  SP  material  consisted  of  a relatively  clean  sand  with  little  or  no 
fines,  the  SC  material  consisted  of  a sand  witn  an  appreciable  amount  of 
fines  (i.e.  a sand-clay  mixture),  the  0L  soil  was  primarily  an  organic 
silt,  the  SM  was  a sand-silt  mixture  with  low  plasticity,  and  the  GM  was 
a gravel-sand-silt  mixture,  (Soil  type  class  3 (Oh)  did  not  occur  in 
the  hand-emplacement  area. ) The  soil  type  map  for  the  hand-emplacement 
area  is  presented  as  fig.  1*9- 

Surface  cone  pene- 
tration  resistance  (CPR) 

180.  Six  factor  classes  of  surface  CPR  were  also  established  for 
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the  mapping  of  the  hand-emplacement  area  as  follows : 


Surface  Cone  Penetration 
Factor  Resistance  (CPR)  Factor 

Class  Value  Range , kPa 


1 

2 

3 

4 

5 

6 


0-200 

200-400 

400-600 

600-900 

900-1200 

1200-1800 


The  mapping  units  portrayed  hy  the  various  patches  on  the  surface  CPR 
map  (fig.  50)  include  the  factor  classes  for  both  the  wet-  and  dry- 
season  conditions , 


Mapping  of  Dry-Season  General  Seismic  Response  Factors 

181.  The  three  dry-season  general  descriptors  (i.e.  PPV,  PSrv, 
and  f as  functions  of  distance)  a.s  determined  by  the  seismic  signal 
generated  by  the  drop-hammer  (controlled)  energy  source  were  used  for 
mapping  the  seismic  response  of  the  hand-emplacement  area.  (Wet-season 
data  could  not  be  used,  as  discussed  in  paragraph  101.  ) Prior  to  map- 
ping, best-fit  curves  representing  the  field-measured  drop-hammer  PPV, 
PSl'Y.  and  f data  (table  8)  were  generated.  The  PPV,  PSPV , and  f 
data  for  selected  distances  were  then  obtained  from  the  best-fit  curves 
(table  8)  and  these  data  were  used  as  a basis  for  constructing  the  seis- 
mic response  complex  map  (fig.  5l).  The  factor  classes  established  for 
the  mapping  of  the  three  descriptors  were  designed  to  be  responsive  to 
the  operational  logic  of  Phase  III  SID’s  (see  fig.  A3,  Appendix  A)  and 
the  distribution  of  the  factor  classes  that  occurred  within  the  hand- 
emplacement  area. 

Seism.'  c response 
factor  mapping  classes 

182.  PPV  and  PSPV  mapping  classes.  The  factor  mapping  classes 
(table  2l)  were  selected  such  that  classes  1,  2,  and  3 would  generally 
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correspond  to  the  Phase  III  SID  threshold  values  on  high,  medium,  and 
low  gain,  respectively.  The  threshold  values  at  30  Hz  are  approximately 
0.08  x 10  ^ , 0.18  x io  ^ , and  0.50  x 10  ^ cm/sec,  respectively  (see 
fig.  A3,  Appendix  A).  These  values  are  relatively  conservative,  since 
they  offer  a reasonable  compromise  between  a frequency  bandwidth  at 
9-45  Hz  and  peak  values  at  approximately  l4  Hz.  Classes  4-8  for  PPV  and 
classes  4-10  for  PSPV  were  selected  to  encompass  the  drop-hammer  field- 
measured  values  (table  21). 

183.  f mapping  classes . The  mapping  classes  for  f were  se- 
lected such  that  they  would  encompass  the  frequency  response  of  the 
Phase  III  SID  and  also  the  field-measured  values.  These  factor  class 
ranges  also  are  listed  in  table  21. 

Mapping  of  seismic  response  factors 

184.  The  seismic  response  factor  complex  map  (fig.  51 ) was  gener- 
ated by  first  preparing  individual  factor  maps  for  PPV,  PSPV,  and  f 
for  distances  of  5,  10,  20,  30,  40,  and  60  m,  respectively.  This  was 
accomplished  by  examining  the  drop-hammer  best-fit  data  (table  8)  for 
PPV,  PSPV,  and  f at  the  various  distances  and  selecting  a factor  class 
for  each  corresponding  distance  (i.e.  a factor  class  for  5,  10,  20,  30, 
4o,  and  60  m,  respectively).  The  factor  class  for  each  distance  was 
then  noted  on  a base  map.  At  this  stage,  boundaries  were  drawn  delin- 
eating zones  or  patches  in  which  all  the  class  numbers  were  the  same. 

In  effect,  the  map  differentiated  the  area  in  terms  of  PPV,  PSPV,  and 
f for  a selected  distance,  such  as  20  m. 

185.  All  the  individual  PPV,  PSPV,  and  f maps  for  each  of  the 
selected  distances  (i.e.  5,  10,  20,  30,  40 , and  60  m)  were  consolidated 
into  a single  map  (i.e.  a factor  complex  map).  All  factor  class  ranges 
identified  on  the  map  in  the  order  of  PPV,  PSPV,  and  f factors  were 
tabulated,  and  each  array  (i.e.  unique  combinations  of  values  of  PPV, 
PSPV,  and  f)  was  assigned  a map  unit  number.  The  resulting  tabulation 
(table  21)  is  the  legend  for  the  seismic  response  factor  complex  map 
(fig.  51);  a unit  on  this  map  represents  an  array  of  class  ranges  as 
identified  in  the  legend,  and  delineates  areas  exhibiting  similar 
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combinations  of  seismic  conditions,  at  least  insofar  as  PPV,  PSPV , and 
f are  concerned. 


186.  A MINISID  III  performance  (detection  distance)  map  of  the 
hand-emplacement  area  (fig.  52)  was  constructed  using  the  experimental 
detection  range  data  (table  9)  obtained  during  wet-  and  dry-season  con- 
ditions at  the  ten  sites  (sites  335-342,  348,  and  349),  air-photo  inter- 
pretation techniques , terrain  data  obtained  within  the  area,  and  the 
interpreter's  personal  knowledge  of  the  hand-emplacement  area.  The  map- 
ping classes  were  selected  to  minimize  the  difference  between  the  mapped 
values  and  those  portrayed  by  the  experimental  detection  range  data. 

The  legend  for  the  MINISID  III  performance  map  is  given  in  table  22. 

The  legend  for  the  SID  performance  map  (tab.le  22)  does  not  differentiate 
between  wet-  and  dry-season  conditions.  As  can  be  seen  in  fig.  30  (a 
bargraph  showing  a comparison  of  the  wet-  and  dry-season  detection  dis- 
tances at  medium  gain),  there  appears  to  he  no  trend  toward  any  one  sea- 
son giving  a consistently  greater  or  lesser  detection  range.  For  ex- 
ample, the  detection  distance  data  presented  for  the  wet  and  dry  seasons 
in  fig.  30  show  a difference  of  1-3  m for  five  sites  (337-339> 

3^2,  and  349),  a difference  of  5-7  m for  two  sites  (34l  and  347),  and  a 
difference  of  18  and  20  m for  two  sites  (336  and  348,  respectively).  A 
more  distinct  trend  may  have  been  observed  if  the  same  man  had  been  used 
for  the  walk  tests  during  both  seasons.  As  this  was  not  the  case,  the 
experimental  MINISID  III  data  for  the  two  seasons'  conditions  cannot  be 
considered  absolutely  compatible.  The  characteristics  of  the  two  per- 
sonnel used  for  these  tests  were:  The  wet-season  MINISID  III  data  are 
based  on  a man  weighing  79-2  kg,  wearing  boots  size  11,  and  walking  an 
average  pace  of  1 m/step;  whereas,  the  dry-season  MINISID  III  data  are 
based  on  a man  weighing  65.8  kg,  wearing  boots  size  9,  and  walking  an 
average  pace  of  1 m/step. 

187.  To  illustrate  the  effects  of  personnel  characteristics  on 
the  MINISID  III  detection  ranges , the  resxilts  of  a few  selected  tests 
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NOTE*.  MAP  BASED  ON  MAN  WEIGHING  66-79  KG,  WEARING  BOOT 
SIZE  9-11,  WALKING  AN  AVERAGE  STRIDE  OF  1 M/STEP. 

FOR  LEGEND  SEE  TABLE  22. 


Fig.  5 2.  SID  detection  performance  map,  hand-emplacement  area, 
a larger  print  of  this  figure,  see  map  U inside  hack  cover) 
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PART  V:  SUMMARY  OF  RESULTS,  CONCLUSIONS,  AND  RECOMMENDATIONS 


Results  and  Conclusions 


188.  As  stated  in  paragraph  37,  the  overall  objective  of  this 
seismic  and  environmental  study  was  to  improve  procedures  for  conducting 
engineering  tests  to  evaluate  various  SID  systems.  The  following  dis- 
cussion will  show  how  well  the  specific  objectives  were  achieved. 

Significant  terrain  and  seismic 
response  factors  and  their  variation 

189.  The  significant  terrain  factors  required  to  determine  sensor 
detection  performance  are : 

a_.  Surface  terrain  factors  (paragraph  119): 

(1)  Soil  type 

(2)  Soil  particle  size 

(3)  Soil  moisture  content 

(4)  Soil  wet  density 

(5)  Soil  strength  (cone  penetration  resistance) 

(6)  Height  and  density  of  vegetation  ground  cover 

(7)  Geometry  of  ground  surface  (roughness) 
b_.  Subsurface  terrain  factors  (paragraph  127): 

(1)  Compression  wave  velocity  (each  layer) 

(2)  Shear  wave  velocity  (each  layer) 

(3)  Thickness  of  e iC'n  layer 

(4)  Soil  wet  density  of  each  layer 

190.  The  relative  effects  of  terrain  factors  (and  changes  in 
them)  cannot  be  specified  in  simple  terms.  However,  it  is  felt  that  a 
sensitivity  analysis  performed  on  the  Fort  Huachuca' data  by  a first- 
generation  model  developed  at  WES  would  indicate  those  factors  that 
cause  significant  changes  in  SID  detection  ranges.  As  a result,  only 
those  factors  indicated  by  the  analysis  would  have  to  be  measured  in 
future  studies  (paragraph  132). 

191.  To  improve  the  overall  procedures  for  prediction,  extrapola- 
tion, and  comparison  of  SID  performances  (detection  and  implantation). 
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meteorological  ar.d  vegetation  data  should  also  be  obtained  within  the 
general  area  of  SID  testing  as  follows : 

a.  Meteorology  factors  (paragraph  76): 

(1)  Wind  speed 

(2)  Wind  direction 

(3)  Precipitation  (rainfall  and  snowfall) 

b . Vegetation  factors  (description  of  plants  within  a 

radius  of  100  in  of  the  SID  test  site)  (paragraph  7^): 

(1)  Average  plant  base  diameter  and  number  of  stems  per 
plant 

(2)  Average  plant  height 

(3)  Average  dominant  plant  spacing 

(4)  Plant  crown  characteristics  (average  diameter  and 
depth ) 

192.  The  seismic  response  parameters  (i.e.  PPV , PSPV , and  f as 
a function  of  distance)  resulting  from  the  drop-hammer  energy  source 
(which  is  assumed  to  be  comparable  to  thul  ix-sulting  from  a footstep 
source)  indicate  that  the  hand-emplacement  area  provides  terrain  condi- 
tions that  are  somewhat  ideal  for  performing  engineering  tests  of  hand- 
emplaced  SID's  (paragraphs  105-107).  The  seismic  responses  as  measured 
at  the  10  test  sites  (see  Appendix  B)  show  some  site-to-site  variation, 
and  this  variation  should  be  considered  in  SID  engineering  testing;  how- 
ever, it  is  believed  that  this  variation  can  be  accounted  for  analyti- 
cally and  therefore  will  not  necessarily  result  in  any  erratic  or  incon- 
sistent results  in  engineering  tests  performed  on  the  hand-emplaced 
SID's.  The  variations  in  seismic  response  parameters  portrayed  by  air- 
drop areas  1 and  4-6,  o.s  measured  in  the  dry  season,  indicate,  in 
general,  that  the  detection  performance  of  air-implanted  SID's  can  be 
evaluated  for  a large  range  of  terrain  conditions  (paragraphs  105-107). 
The  variation  in  these  conditions  allows  the  testing  of  air-implanted 
SID's  under  varying  stress  conditions  (i.e.  high  to  low  rigid-body  de- 
celeration values ) and  under  a variety  of  soil  strength  conditions 
(paragraph  165 ) . 

193.  The  seismic  responses  resulting  from  induced  background 
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noises  created  by  low-flying  fixed-wing  and  rotary-wing  aircraft  are  of 
sufficient  frequency  and  duration  within  the  East  Bange  to  have  signifi- 
cant effects  upon  certain  types  of  SID  testing;  therefore , measurements 
of  the  induced-noise  level  must  he  made  during  SID  testing  (para- 
graph 110).  The  seismic  response  resulting  from  wind  action  on  plant 
structures  within  the  hand-emplacement  area  and  the  airdrop  areas  did 
not  produce  seismic  waves  of  sufficient  magnitude  that  were  considered 
significant  to  the  engineering  testing  of  most  SID  systems  employing  a 
Phase  III  type  or  similar  logic  (paragraph  108). 

194.  Based  on  the  surface  and  subsurface  terrain  data  (end  , in 
some  cases,  seismic  data)  obtained  during  wet-  and  dry-season  conditions 
at  Fort  Huachuca,  it  was  determined  that  the  selected  SID  test  sites 
portrayed  significant  variation  from  site  to  site,  but,  in  general , very 
little  variation  between  wet-  and  dry-season  conditions  (paragraphs  121- 
125,  128-131,  13^-138  and  figs.  32-38). 

Data  collection  procedures 

195.  The  procedures  used  by  WES  for  acquisition  and  recording  of 
the  significant  terrain  and  seismic  response  data  listed  in  para- 
graphs l89  and  191  are  being  documented  in  a separate  loose-leaf  report 
and  will  be  provided  to  Fort  Huachuca  upon  completion. 

Factor  maps  of  hand-emplacement  area 

196.  The  terrain  factor  complex  maps  (figs.  ^7-50)  prepared  for 
the  hand-emplacement  area  provide  surface  aad  subsurface  factor  data  on 
an  areal  basis  for  use  witn  the  forcing  function  and  seismic  wave  propa- 
gation submodels  to  estimate  SID  performance  (l.e.  detection  distances; 
for  man-walking  energy  sources  (paragraph  167 ) • The-  seismic  response 
factor  complex  map  (fig.  51 ) prepared  for  the  hand-emplacement  area  pro- 
vides seismic  data  on  an  areal  basis  on  PPY  , PSPV,  a id  f as  a function 
of  distance,  and  will  be  useful  in  selecting  representative  areas  (i.e, 
areas  that  will  provide  the  least  amount  of  terrain  variation)  for  con- 
ducting DID  engineering  tests. 

Extrapolating  3ID  test  results 

197-  Based  on  the  WES  concept  (fi.g.  2)  for  the  extrapolation  of 
SID  test  results,  adequate  simulation  models  describing  the  phenomena 


under  test  are  being  formulated,  (paragraphs  151  and  1 66).  In  the  case 
of  SID  detection  performance,  this  involves  simulating  the  generation, 
propagation,  and  coupling  of  the  seismic  energy  from  the  source  to  the 
sensor.  In  the  case  of  implantation  of  air-delivered  SID's,  this  in- 
volves simulating  the  terrain-penetration  event  of  an  air-delivered  SID 
and  the  coupling  of  the  SID  to  the  soil.  The  extrapolation  concept  also 
requires  that  the  terrain  inputs  to  the  simulation  models  be  constitu- 
tive and  measurable  properties  of  the  terrain  such  that  they  can  be 
classified  and  portrayed  on  an  areal  basis. 

Depth  of  penetration  and  decel- 
eration  of  air-implanted  sensors 

1.98.  The  predicted  (and  unvalidated)  depth-of-per.etration  results 
for  airdrop  areas  1,  and  5 indicate  that  the  AD3T.D/S  air-delivered 
sensor  will  implant  satisfactorily  (i.e.  it  will  obtain  a penetration 
depth  sufficient  for  detections  of  a seismic  disturbance)  within  airdrop 
area  1 (test  pit  5)  during  both  the  wet  and  dry  seasons;  and  within  air- 
drop area  1 (test  pit  6)  and  airdrop  area  5 (test  pit  2),  it  will  im- 
plant satisfactorily  only  during  the  wet  season  (paragraph  157)-  These 
predictions  are  based  on  the  assumption  that  the  anticipated  impact  ve- 
locity of  the  ADSID/S  will  be  approximately  50  m/sec.  For  airdrop 
area  1 (test  pit  6)  and  airdrop  areas  and  5 during  the  dry  season,  the 
penetration  results  indicate  that  the  ADSID/S  at  50  in/sec  impact  veloc- 
ity will  penetrate  the  soil  only  a few  centimeters  and  then  probably 
fall  down  (paragraph  157)-  All  computed  deceleration  values  were  over 
the  maximum  allowable  for  a workable  ADSID/S,  except  those  for  area  1 
(pit  5)  during  the  wet  and  dry  seasons,  area  1 (pit  6)  during  the  wet 
season,  ajid  area  5 (pit  2)  during  the  wet  season  (paragraph  159)- 

199-  The  very  limited  experimental  penetration  data  obtained  at 
the  CSEWS  airdrop  area  compared  reasonably  well  with  theoretical  predic- 
tions of  penetrations  (paragraph  164 ) . The  comparison  of  these  data 
also  shows  that  the  depth  of  penetration  of  a typical  air-implanted  SID 
(i.e.  ADSID/S)  is  very  dependent  upon  terrain  moisture  conditions  (at 
the  time  of  test)  (paragraph  l6l). 

200.  The  input  parameter  data  required  by  the  equations 
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(presented  in  Appendix  C)  used  to  provide  theoretical  estimu  es  for  the 
depth  of  penetration  and  impact  deceleration  of  air-implanteu  JlD's  are 
obtained  from  special  laboratory  tests  on  undisturbed  soil  samples.  The 
inputs  to  the  equations  are : 

a_.  Initial  Young's  modulus  of  elasticity,  E . 

b_.  Ultimate  shear  strength,  . 

c_.  Compressibility  parameter,  p^/p  • 

<1.  Soil  wet  density,  p , 

e_.  Coefficient  of  friction  between  the  projectile  and  the 
soil,  p . 

f_.  Drag  coefficient,  C^  . 

Recommendations 


201.  Since  the  Fort  Huachuca  hand-emplacement  area  exhibited 
rather  small  differences  in  the  four  subsurface  terrain  factors  (i.e. 
VpjVgjT,  ancl  p)  between  the  wet-  and  dry-season  conditions,  it  is 
recommended  that  measurements  of  these  factors  be  made  only  for  new 
sites  established  for  engineering  testing  of  a SID  system.  It  is  fur- 
ther recommended  that  soil  moisture  content  and  cone  peneiration  resis- 
tance (soil  strength)  be  measured  at  all  sites  (including  the  ones 
sampled  by  WES)  that  are  to  be  used  for  engineering  tests.  To  obtain 
moisture  and  strength  data  on-site  in  the  proper  form  for  use,  a calcium 
carbide  gas  pressure  moisture  tester  (i.e.  Speedy  Moisture  Tester)  and 
the  WES  cone  penetrometer  are  recommended.  (The  procedures  for  using 
this  equipment  are  being  documented  as  noted  in  paragraph  195  and  will 
be  furnished  to  Fort  Huachuca  upon  completion.) 

202.  WES  has  a mathematical  computer  model  (see  reference  3)  that 
provides  estimates  of  moisture  content  and  strength  of  the  0-  to  15-cm 
and  15-  to  30-cm  soil  layers.  It  is  recommended  that  this  model  be  uti- 
lized to  estimate  daily  moisture  contents  and  strengths  for  specific 
areas  to  be  used  for  wet-season  testing  at  Fort  Ruacnuca,  since  local 
rainfall  records  axe  available  for  the  hand-emplacement  and  the  airdrop 
areas,  and  these  data  are  required  as  the  basic  input  to  the  prediction 
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model.  It  is  also  recommended  that  measurements  of  wind  speed,  wind  di- 
rection, and  precipitation  (rainfall  and  snowfall)  be  obtained  during 
SID  testing  for  use  in  interpreting  test  results  (paragraph  77). 

203.  The  input  parameter  data  for  use  with  the  sensor  implanta- 
tion equations,  which  provide  penetration  depths  and  impact  stresses  of 
air-delivered  SID's,  are  presently  obtained  from  the  results  of  special 
uniaxial  strain  and  triaxial  compression  tests  performed  in  the  labo- 
ratory on  undisturbed  soil  samples.  This  procedure  for  obtaining  t'ne 
necessary  parameter  data  is  very  time-consuming  and  too  costly  for  prac- 
tical. application,  such  as  for  use  in  day-tc-day  sensor  testing.  There- 
fore , it  is  recommended  that  research  be  continued  to  determine  better 
and  faster  methods  for  obtaining  input  data  required  by  the  sensor  im- 
plantation equations , as  well  as  by  all  other  models  being  developed  at 
WES  to  estimate  and  extrapolate  SID  performances  for  world  terrain 
conditions . 

204.  It  is  recommended  that  a WES  type  drop-hammer  apparatus  be 
used  by  Fort  Huacnuca  personnel  to  calibrate  a test  site  for  SID  engi- 
neer ing  evaluation.  It  is  also  recommended  that  at  least  three  succes- 
sive drops  of  the  hammer  be  made  at  each  data  collection  site,  and  the 
average  of  the  resulting  seismic  response  data  be  used  for  comparing  SID 
performance  at  one  site  with  that  of  another  site  (paragraph  100). 

205.  Since  the  theoretical  procedures  briefly  presented  in  this 
report  appear  to  offer  a rational  basis  for  interpreting  and  extrapola- 
ting SID  test  results,  it  is  recommended  that,  work  be  continued  on  the 
further  development  of  the  WES  theoretical  seismic  wave  propagation 
model  to  allow  for  the  following  considerations: 

a.  Forcing  function  to  account  for  the  transfer  of  seismic 
energy  from  the  source  to  the  propagating  medium. 

b.  Viscous  damping  coefficients  that  better  approximate 
site  conditions . 

c_.  Forcing  function  that  more  adequately  describes  near- 
ficld  wave  propagation. 

d_.  Coefficients  that  describe  the  geometry  of  the  ground 
surface . 
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Specia.1  studies  to  determine  the  effects  of  sensor  mass 
•sensor  attitude  or  configuration,  and  terrain  coupling 
on  SID  performance.  ••  ? ' . . v is'.  ' . 
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Table  3 

Summarv  of  Moisture  Data 


Moisture  Content,  ,o  1 

Site 

No. 

Depth 

cm 

Wet* 

Season 

Dry** 

Season 

Difference! 

335 

Surface 

5.3 

1.7 

+3.6 

50 

13.2 

10.5 

+2.7 

100 

7.6 

7.4 

+0.2 

336 

Surface 

5.2 

2.2 

+3-0 

50 

6.4 

6.0 

+0.4 

100 

10.8 

9-3 

+1.5 

337 

Surface 

6.5 

1.5 

+ 5.0 

50 

7.7 

7.2 

+0.5 

100 

5.5 

4.2 

+1-3 

33o 

Surface 

11.7 

5.0 

r — 
( 

50 

10.1 

7.1 

+3.o 

100 

9.3 

6.1 

+3.2 

339 

Surface 

5.5 

1.4 

+4.1 

50 

9-1 

7.4 

+1.7 

100 

7.2 

6.1 

+1.1 

340 

Surface 

8.1 

1.4 

+6.7 

50 

3.0 

4.5 

-1.5 

100 

5.6 

6.5 

-0.9 

341 

Surface 

6.4 

0.9 

+5.5 

I 

50 

5.9 

9-7 

-3.8 

100 

7.9 

8.3 

-0.4 

342 

Surface 

5.3 

2.8 

+2.5 

50 

7.0 

5.0 

+2.0 

100 

8.2 

8.3 

-0.1 

* Data  obtained  March  1972. 

**  Data  obtained  August  1971. 

1 Difference  obtained  by  subtracting  dry-season  data  from  wet-season  data, 
+t  Area  inundated  during  wet  season  ( August- October ) . 


Table  U 


Summary  of  Cone-1 enetrat ion  Resi  stance  (Soil  Strength)  D>:ta 


Season 

Cono-l’enclrotion  Resistance,  KT:a,  by 

Various  Pc 

rUiS> 

Site  No. 

Surface 

9-0 

7.5 

10.0 

312.5 

335 

Wet  season 

aO 

2 300 

3590 

3410 

4000 

Dry  season 

100 

3 W0 

U090 

5170+ 

Difference 

+400 

+1090 

- 500 

-1--30 

33b 

Wet  season 

8e0 

24/0 

3220 

2790 

2530 

2480 

Dry  season 

240 

3710 

*>17CH 

Di  fforence 

+L?0 

-1?40 

-1950 

33  < 

Wet  season 

430 

1500 

2300 

2480 

2340 

Dry  season 

540 

3580 

4710 

'•170+ 

Differeuco 

-no 

-2080 

-?4l0 

-20»9O 

330 

Wet  season 

520 

1/.90 

1030 

1800 

1920 

22u0 

Dry  season 

81 3 

?8oo 

5170+ 

Difference 

-290 

-1110 

-3540 

339 

Wet  season 

620 

?4oo 

3>.70 

3720 

'2940 

2300 

Dry  Reason 

790 

4700 

5170+ 

Difference 

+ 30 

-2300 

-1500 

340 

Wot  Season 

83O 

21  /0 

3230 

3210 

2140 

30u0 

Dry  season 

1990 

4560 

51  70+ 

Difference 

-11 80 

-23/0 

-1940 

34 1 

W**t  reason 

4yu 

Irsu 

?7'Aj 

«OnA 

er\J  JO 

-.f  UA 

Dry  season 

'•■10 

4240 

4000 

3940 

31 70+ 

Difference 

-120 

-2520 

-2310 

-2170 

-2340 

342 

Wet  season 

1000 

3000 

4440 

4480 

4490 

4390 

Dry  season 

■•',0 

9000 

51"0+ 

Difference 

+ 350 

-1400 

- 730 

348 

Wet  seas -mi 

oro 

2110 

2940 

3250 

30S0 

2380 

Dry  season 

4250 

51,0+ 

Di  lTcrci.cc 

- 20 

-2140 

-2230 

349 

Wet  season 

770 

2?>.0 

3030 

3010 

2610 

2940 

Dry  season 

410 

3080 

5120 

51  .’0+ 

Difference 

4 300 

-14?0 

-POQO 

-21i,0 

347 

Wet  season 

370 

7.0 

1140 

1450 

1940 

2450 

(Airdrop  1) 

Dry  season 

370 

2050 

4910 

'•I  70 1 

Difference 

+ 1*0 

-I89O 

-3770 

-3720 

410 

Wet  season* 

_ _ 



„ _ 

_ _ 

.. 

(Airdrop  *0 

Dry  season 

30 

370 

1200 

1800 

2010 

2010 

Difference 

— 

— 

„ 

— 

UlJ 

Wet  season' 







„ _ 



(Airdrop  V) 

Dry  season 

320 

1?90 

1900 

1930 

2350 

2*. >40 

Difference 

— 

— 

— 

— 

— 

412 

Wet  season' 

_ „ 



_ _ 

_ _ 



( Airlroi  3) 

Dry  season 

80 

4‘  0 

1230 

2010 

2400 

2810 

Difference 

— 

— 

— 

— 

-- 

-- 

4.13 

Wot  reason* 

_ __ 

.. 

(Airdrop  6) 

Dry  season 

540 

K00 

I9I1O 

1410 

1410 

13'  0 

Pi f Terence 

-- 

_ „ 

Dote:  Wft-ceuson  data  obtain*. 

jd  during;  August  19 '1; 

dry- season  d 

2la  cjl.triiiicd  'luring  Mnrcli  19 

Maximum 

reading  of  cone 

penetrometer  is 

values  prose 

i.tod  are  averages  ol 

17  indi 

penetration  readings. 

Wo  data  obtained  during  wet  season. 


Table  b 

Soil  Data  for  oT»eclai  CGIfrfS  Airdrop  Area* 


lrop 

Moisture 

Cone  Penetration  j 

Airdrop 

Moisture 

Cone  Penetration 

t .No. 

Content, 

Besi stance,  kVa** 

Test  Mo. 

Content,  '/• 

Resistance,  KPa 

3 

6.o 

3630 

2 

-- 

-- 

3 

_« 



it 

8.8 

4400 

5 

6.3 

4930 

C 

7.6 

4340 

7 

5.0 

3930 

8 

6.1 

5140 

9 

6.6 

4920 

10 

6.1 

-- 

11 

6.4 

5:30 

12 

8.6, 

37  JO 

13 

4.7 

4220 

lit 

8.7 

2960 

39 

9.0 

1950 

36 

9.1 

4690 

17 

8.6 

tuYO 

18 

9-9 

4030 

19 

8.3 

1.770 

20 

h.2 

5140 

21 

h.?_ 

5150 

22 

h.C 

5170 

23 

9.3 

5170 

24 

6.5 

4620 

29 

— 

-- 

n 

c.yj 

20i0 

27 

7.1 

34.  V 

28 

10.2 

1670 

29 

11.2 

1030 

30 

10.7 

51701 

31 

-- 

5120 

33 

6.7 

46/V 

8.2 

5050 

34 

5-7 

5170* 

36 

8.0 

5170t 

37 

-- 

5030 

38 

4.2 

5170+ 

39 

7.9 

3310 

ho 

4.3 

5170) 

i . 

*+j 

b 'h 

c 1 1 

J ± 1 

h2 

5.6 

403  0 

43 

3.5 

4420 

44 

9.0 

4 1]  0 

49 

6.7 

5170+ 

A 6 

Y.o 

4360 

47 

5.8 

31704 

48 

9.6 

2630 

4? 

8.6 

51704 

90 

5.1 

5110 

91 

3.8 

51704 

92 

5.0 

4 630 

93 

3-7 

31704 

94  • 

4.6 

5170+ 

65 

8.5 

5 ICO 

96 

4.8 

4910 

57 

5.2 

51704 

98 

7.1 

4 3 00 

99 

6.v 

31  JOh 

6o 

*4 , 6 

51i0 

* Soil  type:  a gravel- sand- silty  mixture-  classified  an  SC  according  to  Unified  Soil  Classi  ficatior 

System;  surface  wet  density  1.7  to  1.9  g/i’m0'  (based  on  10  samples). 

* * A i'lu-3  sign  (O  denotes  maximum  penetrometer  reading. 
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Table  7 


Summary  of  Soils  Data,  Airdrop  Areas 


Air 

tirop 

Aren 

Test 

Pit 

Season 

Relative 
Amount,  % 
Sand  Clay 

0,  deg 

6,  u°g 

Tan  6 

1 

Dry 

53 

47 

0.465 

0 

0 

0 

1 

5 

Wet 

53 

47 

0.465 

0 

0 

0 

1 

6 

Dry 

68 

32 

0.475 

18. 5 

5.94 

0.103 

1 

6 

Wet 

68 

32 

0.475 

0 

0 

0 

4 

3 

Dry 

56 

44 

0.468 

0 

0 

0 

4 

3 

Wet 

56 

44 

0.468 

0 

0 

0 

4 

4 

Dry 

28 

72 

0.450 

8.5 

3.83 

0.067 

4 

4 

Wet 

28 

72 

0,!j?0 

2.0 

0.90 

0.016 

5 

1 

Dry 

l4 

86 

0.443 

8.0 

3.5!+ 

0.061 

5 

-1 

Wet 

l4 

86 

0.443 

0 

0 

0 

5 

2 

Dry 

4o 

60 

0.458 

*“*  r. 

c. . J 

1.15 

0.020 

5 

CL 

Wet 

4o 

60 

0.458 

0 

0 

0 

Note:  0 = angle  of  internal  friction,  deg. 

6 = angle  of  friction  between  projectile  and  soil,  de§*. 
tan  6 = coefficient  of  friction 
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Wet- season  data  obtained  August  197'.:  dry-season  date  obtained  March  1972. 
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Response  Measurements 
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Table  10 


OuiTwiary  of  Vehicle  Seismic  Response  Test  Data,  Hand-Emplacement  Area 


Cross- 

Country 

Si te  Vehicle 


Distance 
f rom 
Energy 
Source 


Peak  Particle 
Velocity 

10"  ^ cni/b?c 


;.  on  so  Measurements  (Dry  Season  ) 
Peak  hummed 


1 article  Velocity 
10”3  cm/sec 


Frequency 7 Hy 


Root  Mean 
Square 
Velocity 
10" 3 cm/sec 


m 

Measured 

Curve 

Measured 

Curve 

Measured. 

Curve 

Field  Men.8 

50 

0.l‘i 

0.15 

0.11 

O.lh 

13 

13-3 

0.04 

100 

o.ofc 

0.10 

O.C6 

0.0/ 

15 

13.3 

0.0? 

150 

e.09 

<j.  07 

0.07 

o.ou 

13 

13.3 

0,0? 

2CO 

O.Gfl 

0.06 

0.07  • 

0.03 

10 

13.3 

0.0? 

?50 

O.OU 

0.05 

0.01 

o.or 

16 

13-3 

0.01 

300 

O.OU 

O.f'lt 

0.01 

0.01 

13 

13.3 

0.01 

Uoo 

•- 

0.03 

— 

0.01 

» 

-- 

V 

500 

— 

0.0? 

— 

0.00 

— 

— 

50 

0.51 

0.5U 

O.UG 

0.53 

18 

18. 5 

0.18 

100 

0.?8 

0.31 

0.2G 

0.30 

21 

17.1 

0.10 

150 

-- 

0.?1 

-- 

o.po 

*- 

-- 



POO 

O.23 

o.lll 

0.31 

O.lH 

0 

H*.  3 

0.09 

250 

— 

0.10 

— 

0.10 

— 

— 

300 

O.05 

0.08 

O.OU 

0.08 

1? 

11.5 

0.0? 

1*00 

O.05 

O.C* 

0.05 

O.O1! 

10 

8.7 

0.0? 

500 

— 

0.03 

— 

0.03 

~ 

— 

— 

50 

.... 

1.10 



1.69 

__ 

100 

0.63 

0.68 

1*35 

o.y 

15 

15. It 

0.38 

150 

— 

0.1*8 

— 

0.58 

-- 



200 

0.33 

0.36 

0.U0 

0.39 

15 

15.3 

0.1? 

?50 

— 

o.pO 

0.57 

..  - 

— 

300 

0.21 

o.?3 

0.13 

0.19 

13 

13.1 

0.05 

400 

0.22 

0.15 

0 . 10 

1? 

11.0 

0.04 

500 

0.12 

0.10 

0.00 

U.05 

10 

10.7 

0.03 

50 

2.95 

__ 

It.  97 

- - 

_ „ 

100 

l.*»3 

I.r-G 

1 .26 

P.37 

15 

12.1 

0.50 

150 

-- 

1.08 

— 

1.31 

— 

— 

200 

1.00 

0.74 

1.3U 

0.77 

15 

14.2 

0.4o 

250 

— 

0.53 

-- 

0.I1O 

-- 

__ 

-- 

300 

0.3? 

0.39 

0.29 

o,?9 

1? 

1C. 3 

0.12 

1)00 

0.?4 

0.21 

0.27 

o.U 

1? 

18.  U 

0.09 

500 

0.11 

0.12 

0.0? 

O.05 

28 

20.5 

0.05 

50 

0.22 

0.18 

0.20 

0.16 

21 

17.  4 

0.07 

100 

0.09 

0.12 

0.05 

0.09 

13 

3-5.1 

0.02 

lvO 

O.OS 

0.08 

0.0C 

0.06 

13 

14.8 

0.0? 

200 

0.07 

0.0G 

O.OG 

O.Oll 

1? 

13* r- 

0.0? 

250 

0.05 

0.05 

0.02 

0.03 

15 

32.2 

0.01 

300 

-- 

0.04 

— 

0.0? 

— 

— 

1*00 

— 

r\  ij  j 

-- 

0.01 

— 

500 

— 

0.0? 

— 

0.01 

-- 

— 

— 

50 



o.r.i 

o.C) 

__ 

100 

0. 20 

0.35 

0.37 

0.35 

13 

lU,s 

0.10 

3 50 

— 

0.?3 

— 

0.P3 

-- 

— 

200 

0. -J9 

0.1-5 

0.29 

o.ir, 

15 

13-0 

0.09 

250 

-- 

o.l? 

— 

0.1? 

— 

-- 

300 

C.07 

0.09 

0.05 

0.09 

1? 

n .5 

0.02 

uoo 

0.05 

0.05 

0.0? 

0.03 

9 

10.0 

0.01 

500 

— 

0.03 

-- 

0.03 

— 

— 

50 

... 

1.2U 



p.llt 



100 

1.24 

0.7? 

1.8? 

1,13 

33 

14.8 

0.0? 

150 

-- 

0.U8 

-- 

0.71 

— 

— 

— 

200 

0.25 

0.35 

0.37 

O.'iY 

15 

13. 0 

0.12 

250 

-- 

o.?r, 

— 

0.3? 

— 

— 

— 

300 

0.16 

0.19 

0.17 

O.P? 

1? 

n.? 

0 . 07 

Uoo 

0.08 

0,11 

0.07 

0.11 

9 

9* 11 

0.0? 

500 

0.07 

c.07 

0.07 

0.06. 

7 

7.6 

0.03 

50 

1 .7*' 

2.35 

1.58 

P.68 

25 

2? » 5 

O.yB 

100 

1.69 

1.3f- 

2.UU 

1-59 

1? 

IG.3 

0.8*  • 

150 

0.90 

0.90 

O.QO 

1.09 

).? 

13.5 

0.3^ 

200 

0.91 

0.64 

1.19 

0.79 

9 

8.0 

0.4  c 

250 

— 

0.1*7 

-- 

O.59 

-- 

*- 

— 

300 

0.30 

0-35 

o.i*-., 

O.I15 

G 

5.1 

0.1  *; 

1*00 

0.17 

0.20 

O.30 

0 . ?7 

7 

7.0 

o/°. 

300 

0.13 

0.12 

0.1U 

0.17 

7 

— 

O.O*) 

¥ 

VI 


Hoi  id - 
**npl  uce- 
acnt 


335 


336 


3>«? 


Table  H 

Smeary  of  MIMTnfP  Data  for  Man-Walk lng  and  Vehicle  Tests 


337 


333 


339 


3**0 


3tl 


3^6 


3**9 


Wo. 

of 

Test 

\a3j  (Inin* 
Dry 
Season 

Med i un 
V.'ot 
Season 

Gain 

Dry 

Season 

High  Gain' 
Dry 
Season 

Man-Walking 

Data 

1 

9.0*  * 

28.0” 

C-3.0" 

? 

JO. 0 

26.5 

63. 0 

3 

10.0 

— 

?3.0 

68, n 

k 

10.0 

— 

27.5 

65.0 

5 

-- 

-- 

27.0 

67.5 

Avg 

: j 

— 

?7.b 

65.3 

i 

11.0-' 

by.  Of 

P8.0-* 

00.0“ 

? 

1P.0 

1*7.0 

3A.0 

(.(  .0 

3 

12,0 

58.0 

3U.0 

CC-.  0 

h 

13.0 

55.0 

35-0 

63.0 

6 

— 

52.0 

35 .0 

— 

Avg 

\?.o 

51.8 

3b.? 

65.0 

1 

3.0“ 

25.  Of 

v3.0” 

2 

5.0 

21.0 

23.0 

52.0 

3 

5.0 

23.0 

T3.0 

52.0 

U 

6.0 

r-r.o 

25. C 

55-0 

Avg 

I-.8 

23.0 

83.8 

53.8 

1 

13.0* « 

PP.<H 

30.0“ 

1*9.0-- 

? 

lo.o 

30.0 

33.0 

50.0 

3 

lf'.o 

3»’.o 

28.0 

c0.0 

U 

ir.o 

31.-0 

25.0 

51-0 

5 

— 

35.0 

-- 

-- 

Avg 

12.3 

31.'. 

29.0 

so.o 

1 

— 

lil.Of 

21.5“ 

to. O' • 

2 

26.0 

20.0 

**1.5 

3 

.. 

2?.0 

22. 

t2.5 

ii 

-- 

2^.0 

23.0 

*n.o 

5 

— . 

25.0 

20.0 

41.0 

Avg 

— 

23.0 

21.li 

»*1.6 

1 

13. 0“ 

__ 

38-0.. 

89.0- ■ 

j 

0.0 

-- 

I10.O 

8*1.0 

3 

11.0 

-- 

38.0 

89.0 

l 

9.0 

-- 

38.0 

75.0 

5 

-- 

-- 

37.0 

89.0 

Avp 

10.0 

-- 

38.? 

86.0 

1 

e.o" 

23.01 

33.0" 

7f  .0*  • 

p 

8.0 

2£.0 

35.0 

71.0 

3 

8.0 

?6.0 

35.0 

bO.o 

t 

-- 

29.0 

29.0 

80. 0 

5 

-- 

-- 

33-0 

99.0 

Avp 

8.0 

27.2 

33.0 

75.2 

1 

i*.o-« 

10.  Ck 

20.0” 

311.0" 

? 

3.0 

?L.o 

20.0 

iaO.0 

* 

1.0 

??.o 

20.0 

38.0 

u 

5.o 

?5.o 

19.0 

3U.0 

5 

3.0 

O’J  A 

Avg 

3.2 

22.  A 

19.8 

36.2 

1 

r«.o*" 

35.ot 

55. 0” 

70.0“ 

p 

l’t.0 

33.0 

5b.  0 

Bl.o 

3 

17.0 

■53.0 

56.0 

79.0 

It 

15.0 

35.0 

56. 0 

76.0 

6 

J8.0 

-- 

69.0 

Avp 

v,,ft 

3b  .5 

55.3 

00. 0 

1 

7.0“ 

27.  f» 

P7.0* • 

? 

( .9 

30.0 

pr, . r.. 

5 6.0 

3 

5.0 

30.0 

27.0 

55.0 

J, 

-- 

27.0 

27.0 

55.0 

•> 

-- 

30.0 

-- 

-- 

Avg 

(.,? 

28.9 

20. *■ 

5I1.R 

j Airdrop  3 


MINISJT)  Deception  Distance,  ^ 


Airdrop  1 3'*  7 


Airdrop  U 


irdrop  3 


Airdrop  0 


CC 

cc 


Wet -Sen sun  data  obtained  in  August  1971,  dry-nenaon  dot  a obtained  in  March  lr>7PJ 
Ho  test  conducted  In  Ur:  vr* 1 r.cocon . 

P?t,rMw  distance  for  e man  w 1 gh ’ ng  6.5 kg  (vearine  boots*  si?'-  9'  and  walking 


0 avciag**  of  one  Ftnp  per  peter. 
Detection  distan-c  for  a man  weighing  79.?  kg  (wearing  toot*,  size  11'  and  walking  an  overn re  of  on**  ntep  p«w  meter. 
Detection  distance  for  o mnn  weighing  70.7  kg  (w-oaring  boots,  r.lr.e  7-1/?)  uni  walking  an  average  of  one  ytej  par  0,83  m. 


K1  77 

vqhi cle , 

speed 

3?  k*n/hr. 

M'.7V 

VPll|e]». 

speed 

0 kjn/hr . 

Ml'  J 

\«’hi 

speed 

1?  fcn/h.  . 

Ml  51 

vehicle , 

speed 

8 Jon/ hi*. 

WO. 

Lou-  Gain 

Medium 

Gain 

High  Gain 

of 

Dry 

V.’ct 

f)rv 

Dry 

Site 

Test 

Season 

Season 

Season 

S'*  a k on 

Man- 

Walking  Data  (Continued ) 

l 

y. 

3'*  7 

1 

10.0” 

29.f>t 

23.0” 

56.0” 

2 

9.0 

Pt,0 

60.0 

3 

30.0 

?3.o 

60.0 

4 

12.0 

-- 

22.0 

59.0 

5 

13.0 

-- 

-- 

-- 

; 1 

Avg 

30.6 

2$  A - 

23-0 

50.8 

i*10 

1 

36.0“ 

30.0” 

6s.  0" 

2 

13.0 

-- 

30.0 

Co.o 

3 

33-0 

— 

35.0 

Cl.o 

j 

t 

36.0 

-- 

31-0 

61.0 

1 

5 

12.0 

-- 

-- 

-- 

A vc 

lb.0 

-- 

31*5 

Cl.O 

l 

kU 

1 

17-0“ 

25.0” 

c*»  ,0*  * 

-4 

2 

17.0 

-- 

25.0 

*15 .0 

3 

17.0 

-- 

25.0 

145.^ 

l 

?0.0 

-- 

?8.0 

hr> , 0 

5 

-- 

-- 

27.0 

1*9.0 

Avg 

37.9 

-- 

?6.0 

U8.8 

112 

] 

33-0“ 

3.3.0“ 

50.0“ 

2 

1U.0 

-- 

bO.  0 

51.0 

3 

33.0 

-- 

3i.  0 

5U.0 

11.0 

-- 

37.0 

9b. 0 

5 

-- 

-- 

31.9 

51.0 

Avg 

3?.B 

-- 

35 .0 

52.0 

'*33 

J 

1*0.  Otl 

-- 

71.0tt 

101, Oft 

, 

2 

br.o 

-- 

6.9.0 

98.0 

3 

I15.0 

-- 

71.0 

105.0 

4 

45.0 

-- 

69. 0 

lt>4  .'J 

5 

hi*. 0 

-- 

-- 

105.0 

• 'J 

Avg 

1*3-2 

”* 

70.0 

102,6. 

Vehi cle 

DaL.a 

1 

200.0 

350.0 

310.0 

1*00.0 

! 

2 

-- 

310.0 

-- 

*450,0 

IWg 

200.0 

?30.0 

310.0 

i*S9.o 

3Wi*f 

1 

390.0 

P98.0 

p8o.o 

1*79.0 

2 

260.0 

300.0 

3?5.0 

500.0 

Avg 

225  0 

299.0 

302.5 

'*87.9 

3*4l|f 

1 

60.0 

120.0 

175.0 

.'6r*  ,f) 

2 

60.0 

315.0 

r.0.0 

275.0 

Avy 

co.o 

117.5 

162.5 

250,0 

3’*bj5 

1 

10.0 

75.0 

90.0 

12-3.0 

. 

O 

1 ' . . n 

A',,o 

75  .0 

1 75.0 

; 

Avg 

12.5 

80.0 

82.5 

130.0 

•* 

3'*5S 

.1 

60.0 

— 

210.0 

250.0 

3'i95S 

1 

15.0 

**5.0 

125-0 

■- 

2 

-- 

__ 

__ 

"1 

Avg 

15-0 

-- 

1*5-0 

2 26.0 

\ 

31.9, 

1 

225.0 

-- 

300.0 

3' 1.5  X* 

■i 

2 

2?5.q 

^00.0 

1*11.0 

Avg 

225.0 

-- 

30c.  (» 

1*03.0 

31.5*1 

1 

-738.0 

-> 

333-0 

ll57.0 

2 

222.0 

-- 

798.0 

nCO.O 

■A 

Avg 

2?,0 . 0 

"" 

295-5 

bpS.  5 

(* 

Length  of  IVne- 
traticm  Klon.’ 
Tert  Sensor  Impact 
No.  Trajectory,  c;r. 


Sensor 

lnpuet  Airy  raff  Wind  Speed, 
An^lc  Xltitudc  nr/r.ce,  and 

*3cg  lyVec  m Direction 


Note:  For  soil  rnositure  and  slrcr'Pth  daea  corresponding  to  each  test  number,  sec  tabl*  t . 

* Direction  of  aircraft,  w ,-nt  tu  east. 

**  Near-surfaoc  impact  velocity  measurements  obtained  using  high~si«ed  ground  photography  (aec  ‘-able  1C}. 


Penetration  Predictions  of  ADSID/S  Sensor 
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Comparison  of  Measured  and  Wed  lo  ted  Detect Distances  (Medium  9gin) 
Using  Analog  Comi  uter  Model?  and JKo n-WnlKlng  Dai. a 


Kite  Nn. 

T)i»toci 
Measured  by 
Mirar.nj 

..ion  Distance,  m 
Predicted  Difference 

Remarks 

Vet-Season  Conditions 

(Aug  l‘>7? ) . IV  v 'v/aag 4 C.* 

370 

64 

to  -4 

371 

105 

97  -8 

High  background  noise  daring  MINTS ID  tests 

37? 

93 

87  -6 

373 

',.6 

74  *18 

High  background  noise  during  analog  man-walking  tests 

y/k 

83 

8 1 *1 

374 

43 

4i  .2 

376 

61 

73  -8 

Wind  justing  during  MUITOTb  teats 

377 

73 

63  -10 

37  3 

37 

3f  -1 

371 

43. 

55  *14 

y 

380 

m0 

27  -13 

Vina  g’rsJ.ing  during  MINI  STD  tests 

3«1 

29 

?s  *1 

38? 

85 

77  -8 

303 

19 

24  .5 

1*9-3 

81 

81  0 

493 

36 

64  *28 

Dry- Season  Conditions 

(Mar  l‘/72).  Port  Bragg,  N.  C.** 

373 

4C 

40  -6 

374 

uu 

54  .m 

37'. 

28 

31  *3 

37f. 

to 

73  *4 

377 

64 

6?  -2 

373 

26 

65  *39 

. 379 

24 

29  -*5 

3S0 

2? 

23  *1 

3S1 

25 

29  *4 

- . 

3B2 

?6 

52  *26 

3«3 

12 

IP  c 

Dry -Season  Condition'- 

(Mur  1-/71),  lort  Huachuca,  Ariz.t 

33'. 

21 

?4  -3 

337. 

34 

?9  -s 

337 

24 

20  -4 

33‘i 

?9 

27  -2 

339 

21 

??  -a 

340 

38 

34  -4 

34 1 

33 

35  *1 

34s 

2C 

2?  -*2 

347 

23 

24  •] 

3U8 

55 

5?  -3 

31,0 

?C 

?8  *2 

4.10 

32 

33  -1 

U11 

2C 

39  *13 

High  background  noise  during  analog  man-walking  tests 
(backhoe  machine  digging  hole) 

4is 

35 

38  *3 

413 

70 . 

7?  1? 

V/et- S'- a son  < Tend  i t i on 5 (dan  3 973  ^ , Oe  rmanyt  t 

4}b 

0 

5 ‘5 

hit 

26 

4c  *3.4 

High  background  noise  during  analog  nan-walking  tests 
(factory,  road  traffic ^ 

hyr 

17 

31  *14 

1*10 

14 

23  -9 

U2(- 

1 

7 *6 

High  wind  noise  during  analog  man-walking  tests 

431 

4o 

53  *14 

4?? 

0 

Jl  4 ll 

4 4 5 

18 

30  -*1? 

High  background  noise  during  analog  man-wnjking  tests 

L4« 

l1* 

33  -19 

High  background  noise  during  analog  man-walking  tests 

4u 

2? 

51  *29 

High  background  noise  during  analog  man-walking  tests 

4 63 

31 

39  *8 

hrs 

3 

5 *? 

Moan  difference,  7.7;  std  d*'v,  *11.0;  number  of  sites  overprodictod , 6;  number  of  sites  underpredicted , 8, 

My  an  difference,  7»1»  std  dev.  *13.?;  number  of  sites  ov^rpredicled , B;  number  of  sites  und^rpredl  cted  , 2. 

Mean  difference,  3-1  m;  ntd  dev,  . 4 m;  number  of  sites  overprcdlcted,  8;  number  of  sites  unde  rprcdic  ted , 7. 

Mean  difference,  12. Oj  std  dev,  7,5;  number  of  sites  overpredicted,  1?;  nun.bei  of  sites  unde rprcdi cted,  0. 


Table  15 


Prediction  of  Maximum  Deceleration  of  ADSID/S  at 
50-m/sec  Impact  Velocity 


Maximum  Deceleration  of 


ADSID/S  at  50-m/sec  0 
Impact  Velocity,  m/sec 


Airdrop 

Test 

Wet 

Dry 

Area 

Pit 

Season 

Season 

- 1 

5 

1.927  x 103 

3.769  X 103 

6 

6.420  x 103 

2.730  x 104 

4 

3 

1.519  x io4 

2.555  X 104 

4 

4 

1.771  X 10 

ll 

2.687  X 10 

c; 

s 

1 

4 

1.175  x 10 

4 

2.355  x 10 

2 

5.968  X 103 

h 

2.408  X 10 

Note:  Maximum  allowable  deceleration  for  ADSID/S  i: 

approximately  9*85  X ic3  m/sec^  . 


Table  16 

Comparison  of  Measured  and  Predicted  Penetration  Depths  of  Air- 


Delivered  (ADSIP/S)  Sensors at  CSEWS  Airdrop  Test  Area 


Impact 

Length 

of  Penetration,  cm 

Impact 

Ang.i  e 

Difference , 

Velocity 

(Off 

Moisture 

Measured 

Test 

Vector 

Vertical ) 

Content 

from 

No.* 

m/sec 

deg 

oj 

Measured 

Predicted  Predicted 

47-795 

37 

6.3 

25.4 

24.4 

-1.0 

46.032 

39 

5.6 

24.1 

22.1 

-2.0 

51.227 

35 

6.4 

28.0 

27.4 

-0.5 

44.603 

28 

8.6 

27.9 

30.0 

+2.1 

52.962 

52 

8.2 

22.9 

30.9 

+8.1 

51.297 

60 

8.0 

26.7 

24.1 

-2.6 

45 . 129 

59 

7.6 

30.0 

18.5 

■12.0 

50.224 

62 

6.4 

24.1 

14.7 

-9.4 

49.104 

55 

7.9 

22.9 

25.9 

+3.0 

49.409 

55 

5-5 

22.9 

17.5 

-5.4 

47.809 

46 

6.5 

27.9 

20.3 

-7.6 

V 

Va 

x 

V 

V = V cos  9 
v a 

v = velocity  in  the  vertical 
direction 

\! 

k 

pl  = ycos 

P ” 1<J  11^  oil 

e 

of  penetration 

V\ 

Pp  - depth  ■ 

of  penetration 

; WiW. 

\ DEPTH  OF 

\ PENETRATION,  Pp 

i n ur  \ 

PENETRATION,  'V  ✓ 

P.  (CM) 

See  table  2 for  site  descriptions. 


Table  17 


Factor  Classes  and  Factor  Value  Ranges  for  Terrain  Factor  Complex 
Map,  Hand- Emplacement  Area,  Wet  and  Dry  Seasons 


Factor 

Class 

Layer 

Thickness  T 
cm 

Compression 
Wave  Velocity  Vp 
m/sec 

Shear  Wave 
Velocity  Vg 
m/sec 

Density  p 
g/cm3 

Layer  1 

1 

40-65 

220-260 

75-115 

1.00-1.45 

2 

66-95 

261-290 

116-135 

1.46-1.60 

3 

96-125 

291-320 

136-155 

1.61-1.80 

4 

126-165 

321-360 

156-175 

1.81-2.10 

5 

166-200 

361-400 

176-190 

2.11-2.30 

6 

201-240 

401-460 

191-220 

7 

241-280 

461-525 

8 

281-320 

9 

321-600 

Layer  2 

1 

100-150 

470-570 

235-285 

1.61-1.80 

2 

151-240 

571-64S 

286-335 

1.81-2.10 

3 

241-330 

646-750 

336-385 

2.11-2.30 

4 

331-480 

751-950 

5 

481-600 

951-1100 

6 

601-800 

1101-1300 

7 

801-1100 

1301-1700 

8 

1101-2000 

1701-2100 
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Table  20 


Comparison  of  Theoretical  Predictions  Based  on 
Site  and  Factor  Complex  Data 


Predictions  Based  or. 


Site 

No. 

Ma.p 

Unit 

Site 

Data  and  20-m 
Deference  Signal 

Factor  Complex 
Data  and  20-m 
Deference  Signal 

Difference,  m 

335 

12 

27 

27 

0 

336 

7 

25 

24 

+1 

337 

25 

17 

15 

+2 

33Q 

13 

23 

28 

0 

339 

30 

17 

16 

+1 

3^0 

3 

25 

25 

0 

341 

22 

31 

31 

0 

3'+2 

21 

19 

19 

0 

348 

9 

27 

28 

-1 

349 

1 

22 

21 

+1 

Table  21 


Legend  for  Dry-Season  Drop-Hammer  Seismic  Response  Map, 
Hand-Emplacement  Area 


Drop-Hammer  Seismic  Response  Factor  Class 


PFV  PSPV  f 

Map  5 1 5 20  50  CO  oO  5 10  20  30  Co  So  5 10  ~20  30  Co  £o 

Unit  m n.  m m m m m m m m ;n  m m m m m m m 


1 

2 

3 

4 

5 

6 

7 

8 


4 

4 

4 

4 

4 

4 

4 

4 


3 

4 
4 
4 
4 
4 


2 

o 
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2 

3 

3 

3 

3 
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2 

2 

2 

2 

2 

2 

3 


1 

1 

2 

1 

1 

2 

2 

o 


1 
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1 

1 

1 

1 

1 

2 
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3 

2 

2 

3 

3 

4 
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o 

2 

2 

2 

2 

3 


2 

1 

2 

2 

2 

2 

2 


1 
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1 

] 

1 

1 

2 

2 
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1 

1 

1 
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1 

2 


1 

1 

1 

1 

1 

1 

1 

2 


3 

5 

4 

4 

5 

4 
9 

5 


3 
5 

4 
4 

9 

4 

7 

5 


2 

4 

4 

4 


4 


3 
5 

4 


2 

3 

3 

3 

3 

3 

3 

4 


2 

3 

3 

2 

3 

3 

2 

3 


2 

1 

2 

1 

1 


Sei  ami  c- 
Response 
Factor 
C]  ass 


Seismic  Response  Factor  Value  Range 


PFV 
10  ^ cm/sec 


P3FV 

_ O 

10  3 cm/sec 


f 

Hz 


1 

0. 

,0 

to 

0.10 

0. 

.0 

to 

0.10 

11 

to 

16 

2 

0. 

.11 

to 

0.30 

0. 

.11 

to 

0.30 

17 

to 

22 

3 

0. 

• 31 

to 

o 

o 

o 

0. 

• 31 

tp 

0.60 

23 

to 

29 

4 

0. 

.61 

to 

1.50 

0. 

.6] 

to 

1.00 

30 

to 

36 

5 

1, 

.51 

to 

3-50 

1 , 

.01 

tc 

1.50 

37 

to 

43 
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3- 

■ 51 

to 

7-00 

• 51 

to 

2.50 

44 

to 

50 
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7. 

.01 

to 

13.00 

2 

■ 51 

to 

3-50 

51 

to 

56 

8 

13 

.01 

to 

28.00 

3 

.51 

to 

6.oo 

57 

to 

62 
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6, 
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to 
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to 

68 

10 

13. 
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Table  22 

Legend  for  SID  Performance  Map,  Hand-Emplacement  Area 
Wet,  and  Dry  Seasons 


Detection 

Di stance 

Class 

Map 

Unit 

Low  Medium 

Gain  Gain 

High 

Gain 

Class 

Detection  Distance 
Class  Range,  m 

1 

1 

3 

5 

1 

1 to  4 

2 

1 

3 

6 

2 

5 to  10 

3 

1 

k 

5 

3 

11  to  20 

4 

1 

4 

6 

4 

21  to  30 

5 

1 

5 

8 

9 

31  to  40 

6 

ti. 

3 

5 

6 

4l  to  50 

7 

2 

4 

5 

7 

51  to  60 

8 

2 

4 

6 

8 

6l  to  80 

9 

2 

4 

8 

q 

61  to  110 

10 

2 

5 

6 

10 

111  to  140 

ll 
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8 

12 

2 
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9 

13 
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i4 
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15 
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16 
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APPENDIX  A;  MINISID  AND  ADSID/S  DEVICES 


1.  The  Phase  III  MINISID  is  a hand-emplaced,  seismic  intrusion 
detection  device,  whereas  the  Phase  III  ADSID/S  (MA37 ) can  either  he 
hand-emplaced  or  emplaced  by  dropping  from  a fixed-wing  or  rotary-wing 
aircraft.  The  major  components  of  the  signal  processing  circuitry  of 
the  MINISID  ar 1 ADSID/S  are  the  geophone  (MINISID  contains  an  external 
cable -attached  geophone  and  the  ADSID/S  contains  an  internal  geophone), 
band-pass  amplifier,  and  an  integration  circuit  (i.e.  a detection  and 
alarm  circuit).  Both  sensing  devices  can  be  disabled  by  the  destruction 
of  internal  components.  Such  disabling  of  the  MINISID  can  be  caused  by 
continuous  transmission,  tampering,  battery  run-down,  end  of  programmed 
life,  and  improper  recovery  techniques.  The  disabling  criteria  for  the 
ADSID/S  are: 

a.  Tilting  (i.e.  improper  implantation)  or  handling.  The 
unit  is  disabled  if  the  implantation  angle  is  greater 
than  60-65  deg  from  the  vertical. 

b.  End  of  programmed  life:.  Time  is  preselected  as  7,  17 , 

28,  ^5,  65,  or  90  days,  ana  an  OFF  position  can  be 
specified  for  its  operational  life. 

c_.  End  of  battery  life.  Approximately  150  days . 
d_.  Tampering.  The  unit  is  disabled  by  an  antitamper  circuit 
if  disassembly  is  attempted  after  impact. 

Since  the  ADSID/S  is  a ballistic  penetration  device,  its  impacting  envi- 
ronment is  critical  to  the  survival  of  its  internal  components.  The 
characteristics  for  survival  of  its  interna  components  of  the  aircraft- 
delivered  device  are  given  below: 

a_.  Shock . Half-size  wave  of  6-m/sec  duration  and  less  than 
1000-g's  p>eak  deceleration . 
b_.  Vibration^,  20  to  2000  Hz  with  20-g’s  peak. 
c_.  Pressure . 5 to  35  psi  or  approximately  1/3  to  2-1/3 

atmospheres . 

2.  The  operational  components  of  the  signal  processing  circuitry 
of  the  MINISID  and  ADSIC/S  devices  are  discussed  briefly  below. 


A1 


Geophone 


3-  The  Phase  III  SID  goophone  consists  of  a coil  of  conducting 
wire  and  a magnet,  one  rigidly  fixed  with  respect  to  the  earth  and  the 
other  suspended  from  a fixed  support  by  a spring.  Any  relative  motion 
between  the  coil  and  magnet  produces  a voltage  across  the  coil  termi- 
nals, which  is  proportional  to  the  velocity  of  the  motion.  The  geophone 

is  considered  to  move  as  .if  it  were  part  of  the  earth,  and  thus  produces 

electrical  signals  whenever  seismic  oscillations  are  detected.  The  geo- 
phor.e  sets  the  low-frequency  operating  point  of  the  system  and,  in  the 
i case  of  a Phase  III  seismic  sensor,  has  a resonant  frequency  of  approxi- 

| mately  l4  Hz  as  illustrated  in  fig.  Al.  The  Phase  III  geophone  is  also 


Fig.  Al.  Response  at  SID  geophone 

underdamped  and  has  a damping  factor  of  approximately  0.2  (i.e.  20  per- 
cent of  critical,  damping). 

Band-Pass  Amplifier 

4.  The  electrical  signal  from  the  geophone  is  of  a very  low  level 
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and  must  be  increased  in  the  band-pass  amplifier  before  further  process- 
ing can  take  place.  Most  SID  amplifiers  can  be  set  for  various  gains 
(low,  medium,  and  high),  which  ore  separated  by  8 db.  The  sensor  ampli- 
fier reduces  some  ui  the  background  noise  and  enhances  the  seismic  sig- 
nal from  the  seismic  sources  by  virtue  of  its  limited  bandwidth  of  from 
9 to  Hz,  as  shown  in  fig.  A2.  These  are  the  frequencies  at  which  the 
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Fig.  A2.  Response  of  SID  band-pass  amplifier 


gain  is  reduced  by  5 db  from  that  at  the  center  frequency.  The  particl 
velocity-frequency  logic  for  low,  medium,  and  high  gains  of  a Phase  III 
SID  is  shown  in  fig.  A3. 


Integration  Circuit 


5-  The  integration  circuit  is  designed  tc-  reduce  false  alarms 
caused  by  background  noise.  The  signal  is  rectified,  and  only  signals 
whose  amplitude  is  0.3  v or  larger  are  allowed  to  be  processed  further. 
This  corresponds  to  a particle  velocity  value  of  38.1  x 10  b cm/sec  or 
larger  at  high  gain  at  a frequency  of  iU  Hz,  as  indicated  by  fig.  A3. 


At  other  frequencies  or  other  gain  settings,  the  value  of  particle  ve- 
locity applied  to  the  geophone  would  he  correspondingly  higher.  The 
SID  alarm  system  operates  in  one  of  two  modes  as  follows:  (a)  the  alarm 

is  set  if  approximately  four  signals  are  detected  within  a 6-sec  period; 
or  (h)  the  alarm  is  set  if  a signal  longer  than  2 sec  is  detected  (sin- 
gle seismic  signals  less  than  2 sec  long  are  ignored;  multiple  seismic 
signals  less  than  2 sec  long  are  processed  in  the  same  manner  as 
footstep-type  signals).  The  maximum  alarm  rate  for  the  SID  is  one  alarm 
every  10  sec.  This  rate  conserves  battery  life  and  still  conveys  the 
needed  information  to  the  monitoring  point. 

6.  The  SID  alarms  are  recorded  on  an  annunciator  or  display  unit 
(portatale).  These  devices  allow  an  operator  to  monitor  a number  of 
SID's  at  the  same  time  and  thereby  assess  the  type  of  seismic  distur- 
bances in  an  area.  In  general,  the  alarm  rate  can  be  used  to  identify 
the  source  of  the  seismic  disturbance  and  the  direction  of  travel  of  the 
disturbance  if  proper  emplacement  is  made.  For  example,  consider  an 
array  of  SID's  emplaced  along  a trail.  If  the  SID's  alarm  in  sequence, 
the  operator  can  identify  the  direction  of  travel  as  being  in  the  direc- 
tion of  the  SID  alarms.  The  velocity  of  travel  can  be  determined  from 
the  rate  at  which  the  alarms  are  recorded.  If  all  the  SID's  alarmed 
simultaneously,  the  operator  can  probably  identify  the  source  as  an 
aircraft  (e.g.  helicopter)  or  exploding  shell.  Unfortunately , with  the 
SID's  presently  available,  interpretation  of  the  alarm  pattern  is  an  art 
rather  than  a science. 
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APPENDIX  B:  TERRAIN  SITE  DESCRIPTIONS  AND  SUPPLEMENTARY 

REFRACTION  SEISMIC  DATA 


General  Site  Descriptions 

1.  The  seismic  and  environmental  data  included  in  this  report 
vere  collected  by  WES  field  teams  at  12  sites  within  the  hand- 
emplacement  test  area  and  5 sites  within  the  6 airdrop  areas.  These 
sites  are  discussed  briefly  in  the  following  paragraphs. 

Hand-emplacement  area  (sites 
335-342,  344 , 345  , 348,  and  349) 

2.  Site  335  wa.s  situated  on  flat  terrain  and  within  an  area  where 
the  vegetation  consisted  of  acacia  and  mesquite  that  reached  heights  of 
1 to  2 m.  The  substrate  material  was  a loose,  silty,  coarse  sand,  clas- 
sified SC  by  USCS  in  the  surface  layer  (0-5  cm)  and  a reddish-brown 
sand-clay-silt  mixture  (GC)  in  the  5-  to  100-cm  layer.  A view  southwest 
along  the  110-m  line  is  shown  in  fig.  Bla. 

3.  Site  336  was  located  on  flat  terrain  with  no  surface  irregu- 
larities. The  vegetation  surrounding  it  was  composed  of  desert  grass 
(species  unknown)  20-30  cm  in  height  and  with  stem  diameters  of  2-3  mm; 
a few  very  widely  spaced  mesquite  bushes  1-2  m in  height  occurred  near 
sta  +80  (fig.  Bib).  The  substrate  material  in  the  surface  layer 

(0-7  cm)  was  a silty  sand  (SC)  and  in  the  subsurface  layer  (7-75  cm)  was 
a reddish-brown  sandy  clay  (CL). 

4.  Site  337 , situated  on  flat  terrain  and  within  approximately 
100  m of  a major  drainageway,  was  surrounded  by  closely  spaced  mesquite 
and  acacia  bushes  that  reached  heights  of  1 to  2 m (fig.  B2a).  The 
substrate  material  was  a coarse,  silty  sand  (SC)  in  the  surface  layer 
0-5  cm.)  and  a reddish-brown  clayey  sand  (SC)  in  the  5-  to  100-cm  layer. 

5.  Site  338  was  situated  on  terrain  that  had  a slightly  downward 
slope  of  approximately  3 percent  (fig.  B2b ) . The  vegetation  occupying 
the  site  was  predominantly  desert  grass  50-75  cm  in  height,  and  occurred 
in  clumps  spaced  approximately  50-75  cm  apart.  A few  widely  spaced 
acacia  bushes,  1 m in  height,  surrounded  the  site.  The  substrate 
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material  in  the  surface  layer  (0-5  cm)  was  a loose  silty-sand  and  gravel 
mixture  (SC).  Several  large  stones  (rocks)  up  to  5 cm  in  diameter  were 
embedded  in  the  ground  surface.  The  material  in  the  5~  to  00-cm  layer 
was  a hard,  compact,  reddish-brown  clayey  sand  with  smell  gravel  (SC). 

6.  Site  339  (fig.  B3a)  was  on  flat  terrain  and  within  200  m of  a 
major  drainageway.  No  surface  rocks  occurred  at  the  site.  The  vegeta- 
tion surrounding  the  110-m  line  consisted  of  widely  spaced  acacia  and 
mesquite  that  reached  heights  of  1-2  m.  The  substrate  material  was  a 
silty  sand  (SM)  in  the  0-  to  4-cm  layer  and  a clayey  sand  (SC)  in  the 
4-  to  200-em  layer. 

7 ■ Site  340 , located  in  grass-covered  terrain  and  within  the  up- 
permost part  of  a drainageway,  sloped  slightly  downward  at  about  2 per- 
cent in  the  direction  of  sta  +100  (fig.  B3b).  In  addition  to  the  grass 
vegetation,  some  rather  large  mesquite  trees  2 m in  height  surrounded 
the  site.  Substrate  material  was  a silty  sand  (SM)  in  the  0-  to  7-cm 
layer  and  a reddish-brown  clayey  sand  (SC)  in  the  7~  to  100-cm  layer. 

S.  Site  34 1 (fig.  B4a)  had  a relatively  smooth  ground  surface  ar.d 
was  approximately  30  percent  covered  with  a low-growing  (height  30  cm) 
desert  grass.  A few  widely  spaced  mesquite  trees  surrounded  the  site. 
The  substrate  material  was  a coarse  sand  (SC)  in  the  surface  layer 
(0-7  cm)  and  a reddish-brown  clayey  sand  with  some  gravel  (SC)  in  the 
subsurface  layer  (7-100  cm). 

9-  Site  342  was  also  situated  in  an  area  covered  with  grass  ap- 
proximately 60  cm  in  height.  The  ground  surface  along  the  110-m  line 
portrayed  no  surface  irregularities.  The  substrate  material  was  a 
coarse  sand  (SC)  in  the  0-  to  4-em  layer  and  a reddish-brown  clayey 
sand  (SC)  in  the  4-  to  100-cm  layer.  A view  of  site  342  from  sta  0 
is  shown  in  fig.  B4b. 

10.  Site  348  was  in  an  area  where  the  topography  was  very  flat  an  1 
the  vegetation  consisted  of  a dense  stand  of  grass  (fig.  B5a),  A few 
widely  spaced  mesquite  trees  1-2  m in  height  occurred  within  75  m of  the 
site.  The  substrate  material  was  a loose,  silty  sand  (SC)  in  the  0-  to 
7-cm  layer  and  a reddish-brown  clayey  sand  (SC)  down  to  a depth  of 
100  cm, 
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11.  Site  3^9  was  covered  (approximately  90  percent)  with  a very- 

dense  grass  (height  approximately  60  cm)  (fig.  B5b ) . The  ground  surface 
along  the  line  sloped  slightly  downward  (approximately  2 percent)  in  the 
direction  of  sta  +100.  and  there  were  a few  large  and  widely  spaced 
me s quite  trees  near  sta  +50.  The  substrate  material  was  a coarse  sand-  | 

silt  mixture  (SM)  in  the  0-  to  7-cm  layer  and  a reddish-brown  clayey 
sand  with  small  gravel  (SC)  in  the  7-  to  100-em  layer. 

12.  Cross-country  site  3M  covered  an  area  100C  m long  and  was 
positioned  from  about  i+00  m southwest  of  site  3*+2  to  250  m southwest  of 
site  339.  The  site  sloped  gently  downward  (approximately  5 percent)  in 
the  direction  of  sta  +100  and  contained  very  few  surface  irregularities. 

Grass  approximately  30-60  cm  in  height  occurred  along  the  entire  length 
of  the  line.  The  substrate  material  varied  from  a silty  sand  (SM)  to  a 
clayey  sand  (SC)  in  the  surface  layer  (0-7  cm).  The  material  below  the 


surface  layer  and  down  to  a depth  of  approximately  100  m was  predomi- 
nantly a clayey  sand  (SC).  A general  view  of  the  site  is  shown  in 
fig.  B6a. 

13 .  Cross-country  site  3>i5  extended  from  about  300  m north  of 
site  3^8  to  about  150  m west  of  site  336.  The  terrain,  in  general, 
sloped  slightly  downward  toward  sta  +1000.  The  desert  grass  that  oc- 
curred along  the  line  varied  in  height  from  about  20  to  75  cm.  A few 
scattered  areas  contained  no  vegetation.  The  substrate  material  in  the 
surface  layer  (0-7  cm)  was  primarily  a clayey  sand  (SC)  and  in  the  sub- 
surface layer  (7-100  cm)  ranged  from  a sandy  clay  (CI<)  to  clayey  sand 
(SC).  A general  view  of  the  cross-country  site  near  sta  +50  is  shown 
in  fig.  B6b. 

Airdrop  area  1 (site  3^7) 

lB.  Site  3Vf  was  situated  in  flat  terrain  and  within  an  area  con- 
taining a rather  close  spacing  of  large  (3  m in  height)  mesquite  trees 
(fig.  B7a).  The  substrate  material  was  a very  loose,  coarse , ^ilty  sand 
(SC)  in  the  0-  to  5-cm  layer  and  a red  sandy  clay  (CL)  in  the  5-  to 
100-cm  layer. 
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j a.  Cross-country  site  3H 

i 


b.  Cross-country  site  3^5 


Fig.  Eb.  Cross-country  sites  3V-+  and  3h5  , 
band-emplacement  area. 
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Airdrop  area  2 (rock 
outcrop  site;  unnumbered) 


15.  Exposed  bedrock  and  boulders  covered  the  ground  surface  at 
tills  site.  A few  small  areas  (i.e.  patches)  contained  a dark  brown, 
coarse,  gravelly  sand  material  (GW).  The  vegetation  occupying  the  site 
consisted  of  grass  and  scattered  clumps  of  acacia. 

Airdrop  area  3 (site  412) 

16.  Site  bl2  (fig.  BOb)  was  situated  on  flat  terrain  and  within 
a large  arroyo  (a  drainageway).  A wash  channel  2 m wide  and  5 m deep 
occurred  within  the  arroyo  and  paralleled  the  110-m  line.  The  vegeta- 
tion was  a dense  stand  of  reeds  (a  type  of  grass)  approximately  1- 

1.5  m in  height.  The  substrate  materials  consisted  of  loose  silty  sand 
(SC)  in  the  0-  to  20-cm  layer,  a coarse  sand  (SW)  in  the  20-  to  45-cm 
layer,  and  a silty  sand  (SM)  in  the  1+5-  to  100-cm  layer. 

Airdrop  area  1+  (si  te  4l0 ) 

IT-  Site  UlO  (figs.  BTb  arid  B9a)  was  situated  on  a slope  of  ap- 
proximately 2 percent.  The  vegetation  was  composed  primarily  of  creo- 
sote hushes  approximately  1-1.5  m in  height;  their  average  spacing 
varied  from  1.5  to  3*0  in.  The  sustre  -e  material  consisted  of  a gray, 
loose  silty  sand  with  pebbles  (0-  to  T-cm  layer),  and  a firm,  light- 
brown  clay  (CL)  with  scattered  caliche  down  to  a depth  of  100  cm. 
Airdrop  area  5 (site  til) 

18,  The  ground  surface  at  site  1+11  varied  from  flat  (from  sta  0 
to  sta  +50)  to  sloping  (sta  +50  to  sta  +100).  The  vegetation  was  com- 
posed of  small  woody  plants  (acacia)  that  reached  heights  of  1 to  2 m 
and  had  stem  diameters  of  1-2  cm.  The  acacia  vegetation  covered  ap- 
proximately 80  percent  of  the  ground  surface  and  made  walking  through 
the  area  very  difficult.  The  substrate  material  in  the  surface  layer 
(0-15  cm)  was  predominantly  a silty  sand  with  white  caliche,  and  con- 
tained some  rock  fragments  (figs.  B8a  and  B9b ) . The  material  in  the 
subsurface  layer  (15-100  cm)  was  a red  clay  (CH)  and  some  white  caliche 
Airdrop  area  6 (site  4l3) 

19-  Site  413  (fig.  B10)  was  situated  in  a rather  large,  flat 
lake  bed  (Willcox  Playa).  The  substrate  material  in  the  surface  layer 
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| Fig.  BIO,  Airdrop  area  6,  site  413 

! (0-20  cm)  was  a silty  clay  (CL)  and  in  the  subsurface  layer  (20-100  cm) 

j was  a dark  gray  clay  (CH). 

| CHEWS  airdrop  area 


20,  The  CSEWS  airdrop  area  was  situated  on  an  unimproved  air- 
strip (fig,  Bll ) . The  substrate  material  consisted  of  a sandy  silt  to 
clayey  sand  (SC)  with  various  amounts  of  gravel  and  small  pebbles, 

Refraction  Seismic  Data 


21.  Wet-  and  dry-season  refraction  seismic  data  for  sites  335- 
342,  3^5  (wet  season  only),  347-349,  and  4lQ-4l2  (dry  season  only)  with- 
in the  hand-emplacement  and  airdrop  test  areas  at  Fort  Kuachuca  are  in- 
cluded as  figs,  B12-B37.  These  data  are  portrayed  in  the  conventional 
time -versus -distance  format..  . 
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Fig.  Bl8.  Wet-season  refraction  seismic  data,  sice  3tl 


Fig.  B19.  Wet-season  refraction  seismic  data,  site  3^2 


, ih1  '.’>*•**  '.*■  ‘ 


6 


12 


. B30. 


Distance,  meters 

Dry-season  refraction  s< 
site  3ni 


li.l  ...  I. . .11  -.1  l . . J i . . I 


APPENDIX  C:  EQUATIONS  FOR  PREDICTING  PENETRATION  DEPTHS 

AND  DECELERATION  VALUES  FOR  AIR- IMPLANTED  SENSORS 


1.  The  Fundamental  concepts  of  the  mathematical  penetration  model 

used  to  predict  the  depth  of  penetration  of  the  ADSID/S  sensor  into  the 

Fort  Huuchucu  terrain  materials  are  documented  in  this  appendix.  The 

penetration  model  is  based  on  the  dynamic  cavity  expansion  theory  for  an 

elastic-plastic  locking  medium  formulated  by  Ross  and  Hanagud,"^  and  the 
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soil  penetration  theory  derived  by  Brooks  and  Reis.  ' The  dynamic  cav- 
ity expansion  theory  has  been  used  successfully  in  penetration  studies 
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of  ice,  frozen  ground,  and  various  types  of  soils.  The  soil  pene- 

tration theory  of  Brooks  and  Reis,  modified  to  account  for  the  compres- 
sibility of  the  substrate  materials  by  utilizing  the  expression  for  the 
dynamic  pressure  inside  a spherical  cavity  as  derived  by  Ross  and 
Hanagud,  was  used  in  this  study. 

2.  By  assuming  a simple  elastic-plastic  constitutive  model  for 
soil.  Brooks  and  Reis  derived  the  following  differential  equation  of 
motion  for  an  ogival  projectile  penetrating  into  a homogeneous  soil  me- 
dium in  the  vertical  Direction  z . 


2 

r.i  d z 


+ 


* I v(i)‘ 


= 0 


(Cl) 


m = muss  of  projectile 
A = cross-sectional  area  of  projectile 
d2z 

— — - = acceleration  of  prnjpt't.ile  at.  time  t 
dt 

= ultimate  strength  of  substrate  material  (.equivalent  to 
Von  Mj ses  yield  limit) 

E = initial  elastic  modulus  of  substrate  material  (equivalent  to 
Y oung ' s mo  dul us ) 

Cp  = drag  coefficient 

p = initial  mass  density  of  soil 

-rr  = velocity  of  projectile  at  time  t 


Cl 


3.  Using  the  Newtonian  form  of  aerodynamic:  drag.  Brooks  and  Reis 
derived  the  following  expression  for  the  drag  coefficient  in  terms 
of  the  geometry  of  the  frontal  faec  (fig.  Cl)  of  the  projectile  and  the 
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Fig.  Cl.  Projectile  parameters 
frictional  forces  acting  on  it. 


where 

R = radius  of  ogive 

P 

= the  caliber  ratio  of  the  ogive 
o 

r = radius  of  projectile 
a = ogive  angle,  see  fig.  Cl 

p = coefficient  of  friction  between  the  projectile  and  the  soil 


1.  The  differential  eq  Cl,  in  conjunction  with  the  expression  l’or 
the  drag  coefficients  C (eq  CP),  w;ls  used  to  predict  the  penetration 
of  the  ADGIP/S  into  earth  materials.  The  ultimate  strength  of  the  sub- 
strate material  , however,  was  modified  to  account  for  the  compres- 
sibility of  the  material.  The  following  initial  and  end  conditions  were 
used  to  integrate  eq  Cl: 


— = V . when 

at  1 


dz 

at 


= 0 


when 


(C'3) 


whore 

V = impact  velocity  of  projectile 

Z = maximum  depth  of  penetration 
max 

5-  To  account  for  the  compressibility  of  the  material,  the  ultj- 

0/3 

mate  strength  was  modified  using  the  expression  o [(P/3)(E/o  )]*"'  in 

t X \1 

eq  Ci,  (This  expression  corresponds  to  the  quasi-static  pressure  in  a 
spherical  failure  zone  at  the  interface  of  the  projectile  and  the  soil.) 
This  was  accomplished  by  equating  the  expression  above  to  a correspond- 
ing expression  derived  by  Ross  and  ilanagud  l’or  a compressible  elastic- 
plastic  material.  Assuming  an  elastic-plastic  material  of  the  looking 
type  characterized  by  the  pressure-volumetric  strain  relation  and  shear- 
ing stress-shearing  strain  relation  shown  in  fig.  C2 , Ross  and  Hanugud 
derived  the  following  expression  for  the  quasi-static  pressure  inside  a 
spherical  cavity. 


D = - E[1  - exp  (-33)]  - j YZn6  + |y  ^ Ktn 
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The  parameters  in  eq  aj’e  as  follows : 
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6=1-—  exp(-33) 
PP 


(C5) 


n=  I 


n=l  n 
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Pig.  C2.  Idealized  stress-strain  curves  for  a locking 
elastic-paastic  material 
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P ~ 2E  " 3 
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P = p exp(e^) 


(C8) 


where  (for  eq.  C*+-C8) 

B = Young's  modulus  of  elasticity,  corresponding  to  locked  elas- 
tic region  in  fig.  C2,  = 3G,  where  G is  the  shear  modulus 

Y = yield  strength  of  substrate  material,  shown  in  fig.  C2 

n - dummy  variable 

E = strain-hardening  modulus,  corresponding  to  locked  plastic 
X'  region  in  fig.  C2,  - 3G 

p = initial  density  of  substrate  material 

e.  = volumetric  strain  related  to  the  elastic  region  of  the 
pressure-volumetric  strain  curve  in  fig.  C2 

= volumetric  strain  related  to  plastic  region  of  the  pressure 
volumetric  strain  curve  in  fig.  C2 

6.  It  is  to  be  noted  that  due  to  the  locking  characteristic  of 
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the  volumet.ric-strnin-pressure  relation  in  fig.  C2 , the  built  modulus  of 

the  material  in  both  the  elastic  and  the  plastic  range  is  infinitely 

large,  and  the  shear  modulus  G is  equal  to  E/3  • If  the  parameters 

e.  and  e in  eq  C4-C8  are  taken  uo  be  zero,  an  incompressible  ma- 
i P 

terial  will  be  simulated.  To  simulate  a perfectly  elastic-plastic  ma- 
terial, the  strain-hardening  parameter  E+  should  be  equated  to  zero. 

7.  In  order  to  be  compatible  with  the  Brooks  and  Reis  assumption 
of  perfectly  elastic-plastic  material  in  utilizing  eq  C4,  the  strain- 
hardening parameter  E was  taken  to  be  zero.  The  system  of  eq  C4, 

C5,  C7 , and  08  was  then  studied,  parametrically  to  determine  the  effects 
of  compressibility  of  the  material,  on  the  pressure  inside  the  cav- 

ity. It  was  found  that  for  the  range  cf  variables  considered  in  this 
study,  the  pressure  Pg  calculated  for  an  incompressible  material  { i . e . . 
e = = 0)  is  approximately  Pp/P  times  the  corresponding  value  for 

the  compressible  material.,  i.e. 

P 

(P, . ) incompressible  « ^ (P  , ) compressible  (C9 ) 


Now  equating  the  two  expression  for  the  quasi-static  pressure  and 

p ° 

ou[(2/3)(E/ou)]  ''  and  solving  for  yields 


o 

u 


(CIO) 


In  view  of  eq  C 9 and  CIO,  the  modified  value  of  the  ultimate  strength  of 

the  material,  denoted  by  (a  ) compressible  , then  becomes 

Vi 


(o  ) compressible  = zr  (o  ) 

V-’i  u 

For  an  incompressible  material  p^/p  = 1 and  no  modification  of 
necessary. 


(Oil) 


u is 
u 


8.  Eq  Cl,  C2,  and  Cll  contain  eight  parameters  that  define  a 
terrain  impact  condition  for  a prediction:  m , A , E , °u  ’ ^p  ’ 
p , p , and  . The  parameters  E and  Ou  for  input  to  the  equa- 
tions were  determined  from  the  triaxial  compression  test  data  (fig.  18) 
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presented  in  Part  III  of  the  main  text.  The  parameter  p was  deter- 
mined from  the  uniaxial  strain  test  results  (fig.  IT),  also  presented  in 
Part  HI  of  the  main  text,  by  idealizing  them  with  the  locking  mod;  1 

shown  in  fig.  C2.  The  friction  coefficient  p was  determined  from  the 

9 

triaxial  test  results  using  the  procedure  given  by  Potyondy.  Th<  fron- 
tal face  (nose)  of  the  ADSID/s  was  approximated  as  a 4.35  caliber  ogive. 


APPENDIX  D:  COMPUTER  MODELS  FOR  PREDICTION  OF 

SID  DETECTION  PERFORMANCE 


Theoretical.  Detection  Performance  Computer  Model 

1.  The  overall  theoretical  detection  performance  computer  model 
used  in  this  study  is  composed  of  two  submodels,  the  seismic  wave  propa- 
gation model  and  the  sensor  logic  model.  The  detailed  steps  in  the 

formulation  of  this  model  will  not  be  presented  in  this  report,  since 

IT 

they  have  been  documented  in  two  other  reports.  ’ A brief  descrip- 
tion of  the  model,  however,  is  presented  in  the  following  paragraphs. 
Seismic  wave  propagation  model 

2.  A general  expression  has  been  developed  for  predicting  the 
characteristics  of  a seismic  signal  (i.e.  particle  velocity,  particle 
acceleration,  or  particle  displacement)  a.t  selected  ranges  from  a speci- 
fied energy  source  such  as  a walking  man,  moving  vehicle,  or  other  dis- 
turbances . The  general  equation  is 


At 

L ,p  m 


cos 


w 

V + * ♦ | (P  - i)  - ^ (Bm  - up  * eg 


(Dl) 


where 


A ( R ) 
li,p'  m 


particle  motion,  in  cm.  nt  a selected  ranere  R 
7 ' ~ w ] 

and  time  t 


L = vertical  component  of  three-  dimensional  seismic 
wave 


p = 0,  1,  2;  0 for  particle  displacement,  1 for 

particle  velocity,  and  2 for  particle  acceleration 

R = range  at  which  the  new  signal  is  desired.  The 
range  is  assumed  to  be  jn  the  far  field  (i.e. 
farther  than  2-1/2  wavelengths  from  the  source). 

t = time,  sec 


Dl 
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n 


<*> 

n 


e 

i 


frequency  numbers , A maximum  of  102U  frequency 
data  points  are  normally  used  to  describe  a 
particle  velocity  curve. 

circular  frequency,  radians/sec 

Fourier  transform  of  a measured  particle  motion, 
in  cm,  at  range  and  time  t 

magnitude  of  surface  geometry  coefficient 

base  of  natural  logarithm 

viscous  damping  coefficient;  assumed  to  be  equal 
to  0.03 

Bayleigh  wave  velocity,  m/sec 


<{>  = phase  factor  associated  with  p 

0_  = phase  factor  associated  with  S 
S - n 

3.  The  inputs  required  by  the  seismic  wave  propagation  model 
(eq  Pi)  are  as  follows: 

a.  Source  signal  of  stress  applied  to  the  ground  as  a func- 
tion of  time  for  the  specified  energy  source  (walking  man, 
moving  venicle , etc . ) . 

b.  Terrain  data  on  V0  , V0  , p , and  T for  each  sub- 

— -T  O 

strate  soil  layer. 

c_.  The  viscous  damping  coefficient  £ ; the  coefficient  used 
in  conjunction  with  the  read  Rayleigh  wave  number  k 
(where  k = 2tt/X  , and  X = wavelength  to  account  for 
viscous  damping  in  the  medium),  i.e.  k = k(l  - i£)  , 
where  k is  the  complex  k * 

The  output  from  the  model  is  the  particle  velocity,  particle  acceD.era- 
tion,  or  particle  displacement  as  a function  of  time  for  various  se- 
lected distances  from  the  energy  source.  Predictions  are  normally  ob- 
tained for  each  3-m  distance  to  some  maximum  selected  value.  For  the 
Fort  Huachuca  study,  predictions  of  the  vertical  component  of  particle 
velocity  as  a function  of  time  were  made.  The  particle  velocity-time 
data  were  then  used  as  input  to  the  seismic  sensor  logic  model  for  pre- 
dicting SID  performance  as  discussed  in  the  following  paragraph. 
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Seismic  sensor  logic  model 

4.  The  seismic  sensor  logic  model  is  composed  of  a series  of  sub- 
routines that,  in  general,  utilise  transfer  functions  to  simulate  the 
effects  of  the  frequency  response  of  each  of  the  electronic  components 
of  a SID  on  the  seismic  signal  received  by  the  geophone.  (The  outputs 
from  the  seismic  wave  propagation  model  are  the  inputs  to  this  model.) 

5.  SID  geophone  transfer  function.  The  equation  for  the  transfer 
function  for  the  SXD  geophone  was  developed  from  attenuation-phas^  plots 
(Bode  diagrams).  This  method  is  based  on  the  relation  between  phase  and 
the  rate  of  change  of  gain  with  frequency.  The  equation  used  for  the 
Phase  III  SID  geophones  is 

Vif)2 


rpp  = == 

1 (if)2  + s.6( if)  + 196 


where 


TF.  = transfer  function  for  the  Phase  III  SID  geophone 

= arbitrary  constant  associated  with  geophone  calibration 
i = V^T 

f = frequency,  Hs 

6,  SID  band-pass  amplifier  transfer  function.  The  equation  de- 
veloped for  the  transfer  function  for  the  band-pass  amplifier  is 


= 


Kg  (if)' 


2 [(if)2  + 12.726(if)  + 81] [ (if )2  t G3.63(if)  + 2025) 


where 

TFg  = transfer  function  for  the  Phase  III  SID  band-pass  amplifier 
K2  = arbitrary  constant  associated  with  the  amplifier  gain 

7.  Scientific  geophone  transfer  function.  Because  the  input 
seismic  signal  in  the  particle  velocity-time  relation  was  obtained  with 
a scientific  geophone,  the  response  characteristics  of  the  scientific 
geophone  had  to  be  taken  into  account.  The  transfer  function  for  the 
scientific  geophone  is 

K (if)2 

IV  — 

3 (if)2  + 7.07(if)  + 25 


where 


i.  I 


t • 

j 

ii  * 

1 

\ ! 

t i 


TF^  - transfer  function  for  the  scientific  geophone 

K_,  = arbitrary  constant  associated  with  scientific  geophone 
3 calibration 

8.  The  theoretical  combined  transfer  function  for  the  geophone 
and  band-pass  amplifier  of  the  Phase  III  SID  and  the  scientific  geophone 
is  given  below 


TF  = 


(TF1)(TF2) 

tfT 


! 


■ i 


I 


\ 

* i 
I I 


f 4 

l 


or 


TF  = 


K(if)2[(if)g  + 7-0T(if)  + 25} 


;(il’)2  + 5 ,6(if ) + 196][(if)2  + 12.726 (if ) + «l][(if)2  + b3.63(if)  + 2025] 


where 


K = 


K1K2 


The  value  of  K was  determined  to  be  compatible  with  the  alarm  ampli- 
tude threshold  values  fov  the  Pha.se  J.II  SII>«  The  alarm  levels  are 
38.1  x 10“fa,  95.2  x 10"13,  and  238.1  x io_D  cm/sec  for  high-,  medium-, 
and  low-gain  settings,  respectively  (see  fig.  A3,  Appendix  A).  A 
graphic  presentation  of  the  transfer  functions  describing  the  effects 
of  geophones  and  the  band-pass  amplifier  on  the  seismic  analog  signal 
is  given  in  fig.  Dl. 

9.  The  output  from  the  seismic  sensor  performance  model  is  the 
activation  range  for  low—,  medium—,  and  high— gain  settings  of  the  dID. 
The  activation  distance  obtained  may  nox  be  the  distance  as  which  the 
energy  source  is  first  detected.  This  results  from  the  fact  that  the 
SID  does  not  transmit  a radio-frequency  (RF)  signal  until  the  logic  in- 
tegrates a number  of  signals  from  succeeding  energy  sources,  Thus,  the 
actual  detection  distance  may  be  3 to  5 m less  than  the  value  for  SID 
logic  activation  (that  is,  if  the  source  is  moving  toward  the  SID).  If 
the  energy  source  is  moving  tangential  to  the  SID,  the  distance  at  which 
the  SID  logic  is  activated  will  be  approximately  the  distance  at  which 
the  energy  source  would  be  detected. 
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Analog  Computer  Model 

10.  The  analog  computer  used  to  obtain  predictions  of  SID  detec- 
tion performance  was  a Systron  Donner  Model  SD/80.  A schematic  diagram 
of  the  computer  system  is  shown  in  fig.  D2. 

11.  The  analog  model  was  built  to  be  flexible  and  to  allow  simu- 
lation of  the  wave  shapes  actually  observed  in  a Phase  III  SID.  The 
basic  components  of  the  analog  computer  model  are  as  follows : 

a.  An  input  voltage  comparator,  with  a trigger  level 

corresponding  to  each  respective  SID  gain  setting  (low, 
medium,  or  high). 

b_,  A pulse-shaping  network  , which  forms  a pulse  whose  area 
is  proportional  to  the  seismic  input  signal. 
c_.  A pulse-integrating  circuit,  which  stores  pulses  from  the 
shaping  network  for  a preset  time  interval. 
d_.  A second  voltage  comparator,  with  a trigger  level  cor- 
responding to  the  detection  resulting  from  four  seismic 
signals  in  C sec. 

12  The  analog  computer  was  calibrated  with  a damped  30-Hz,  sine 
wave  whose  peak  positive  v vine  (amplitude)  times  the  appropriate  gain 
setting  (low,  medium,  or  high)  equals  the  threshold  level  for  the  first 
comparator.  This  sine  wave  was  programmed  to  cycle  four  times  in  a 
period  of  6 sec  to  simulate  the  minimum  signal  that  would  trigger  (i.e. 
activate)  a Phase  III  SID  logic  circuit.  The  low-  and  medium-gain 
logics  were  programmed  separately  on  the  computer  to  facilitate  calibra- 
tion and  system  operation  (the  SID  logic  corresponding  to  high  gain  vas 
not  programmed  at  the  time  of  this  report). 

13.  The  computer  model  was  operated  by  inputting  the  M analog 
tape  recordings  of  the  seismic  response  generated  by  a walking  man  or 
moving  vehicle.  The  output  (see  example,  fig.  D3)  was  in  the  form  of 
an  x-y  plot  that  displays  the  position  of  the  energy  source  along  the 
test  line.  Successive  event  markers  on  the  plot  indicated  preselected 
distances,  such  as  1 m,  10  m,  50  m,  and  so  on,  A d-c  voltage  shift 


GEOPHONc  I AMPLITUDE 


Schematic  of  analog  computer  model 


medium- gain  trigger 

POSITION  (ACTIVATION 
DISTANCE  - 72  m) 


NOTE:  PREDICTION  SHOWN  IS  FOR 

M 1 5 I VEHICLE  (SPEED  8 KM/HR) 


Fig.  D3.  Example  of  output  of  analog  computer  model 

(line  displacement)  on  the  plot  indicated  that  the  SID  had  triggered. 
The  detection  range  was  then  determined  by  measuring  the  horizontal 
length  from  the  start  of  the  line  to  the  point  of  displacement. 
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V.  . An y Tv.  t and  Evaluation  Cocra-sand 

V . f . Aivsv  Fleet ronico  I roving  Ground 
Fort*  Huachuca,  Arizona 


1 a AbiTRACT  70  obtain  a quantitative  basis  lor  experiment  design  dec i:* ions  regarding  evaluation  of  seismic 
Intrusion  detector  (.’ID)  ey&toi-is,  test  urea.;  were  selected  (one  for  hand- emplaced  fib’s  and  six  for  air- 
: no  landed  JXb'r)  :jid  special  reiejr.io  and  environmental  Field  experiments  voi\  conducted  at  Fort  Kuachuca  dur- 
i i‘Z  Au.j  19*71  aim  Mar  197- * Data  obtained  were  assumed  to  be  represent  olive  01  vet-  and  ory-tv-i  ruin  conditions 
at  the  tort  sites  . The  M! lection  of  the  0J7j  1 .crt  areas,  their  vet-  .*nd  dry- terrain  cluaracteris’icr , and  the 
orooedure.*  for  collecting  environmental  data  and  conduct. in- ' seismic  response  tests  (i.c.  mun-v all lug,  drop- 
's nrs  icr,  M151  wheeled  and  M:>77  trach.ed  vehicles,  ambient  ur.d  induced  noises)  are  described.  Also  presented 
are  experiment ni  detection  distance  data  on  the  h:uul- emplaced  miniaturized  fir  (:niil.;lh),  :c:d  experimental 
; d n on  the  depth  of  penetration  and  an -le  of  impact  of  the  air- delivered  rcisnic  intrusion  dot ector/.'hort 
( AT? Xv/z ) . A theoretical  system  for  modeling  tlic  quantitative  eflcetr  of  the  terrain  on  L‘l.  detection  per- 
fomioncc  and  methods  by  vliicl:  fib  tert  results  cun  bo  extrapolated  iron  one  sire  condition  to  another  -iro  also 
prerented.  This  system  Lor  been  implemented  on  the  Watcrwayr.  Experiment  station  computer  and  iuvolvvs 
acquisition  of  terrain  and  ceir.mic  data  for  input  to  the  model;  exercising  the  model  for  each  set  of  terrain 
factor  data}  mid  port  ravine  tnc  perlomanco  predictions  for  each  set  of  terrain  arid  tvienic  fw-loi  daiK*  A 
r’/r tei;i  of  equations  for  lieteminiiu;  soil  penetration  and  deceleration  of  air- delivered  fU',r  ir  included.  The 
•yster.  employ r.  u r.ericr,  or  conputcr  routiner  which  treat  the  soil  as  an  elastic- plastic,  struin-hardeninc? 
rompror-dblc  medium.  Input-  for  tlie  penetration  equations  are  dynamic  r t. re r.r- strain  data  obtained  from  lab- 
oratory triaxial  ccriprumiou  an- 1 uni..xial  r-tnin  torts.  The  liond-cj^jliceia-nit  area  I'rovidcd  near-ideal  terrain 
ronditiour  for  enc^ ncerinp  tents  of  hand- emplaced  Gib’s.  Gui-face  and  subsurface  terrain  data  therein  varied 
:orc  from  site  to  site  than  between  vet  and  dry  season-  , dei.nie  ror-ponser  neusui-eU  at  the  1 o test  rites, 
it.hin  the  } 1 and-  eiiip  lac  a .ic  1 it  area  chow  some  rile- to- site  variation,  which  should  be  considered  in  fit’  en£i- 
nceriiiG  testin'.  This  variation  cun  l-e  accounted  for  analyticuUor;  thcrefojx*  it  will  not  necessarily  result 
in  any  erratic  or  inconsistent  ia.»rultr  in  enytneerinc  testin'  of  t]ic  hand- emplaced  flP’c.  Variations  in 
-eis.uic  response  j»araiwter-  portrayed  by  airdrop  areas  1 ami  '.-i,  measured  in  the  dry  season,  generally 
indicate  that  Uic  detectio..  porfon:auico  of  air-implsuitcd  SIl>’s  cai,  be  ovulutiled  for  a large  range  of  terrain 
:ouditlons  . The  variation.-  in  there  condit.ionr.  allow  the  testing  of  air-iiiplojitcd  f-lD’s  under  vaiyiiv 
t.'orr  conditions  (i.e.  hir):  to  lo’-’  rigid-bod,.-  deceleration  valuer)  and  various  soil  strcngll:  conditions. 

"Kc  theoretical  detection  prediction  scheme  and/or  conccx-tr-  provide  an  anoXvtical  f ra; icworh  X'or  predicting 
. II-  performance  for  m:m  valhiiy.  taryclr  in  various  world  (Continued) 
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terrain.';.  Tnu  •.•J'esr.e  requires  specific  terrain  and  seismic  input  data  l'or  execution,  and  hc  a result,  four 
typer,  of  naps  were  prepared : terrain  factor  complex  nape  for  subsurface  layv it  1 an J f?  to  provide  inputs  to 

the  theoretical  seismic  wave  propagation  model j surface  condition  maps  to  be-  used  with  the  fore ing  fine t ion 
model  under  development  J seismic  response  factor  complex  mans  to  show  the  relative  seismic  res  police  of  the 
hand- emplacement  area;  and  SID  oiv-man  walking  detection  performance  map.  The  predicted  depths  of  penctra- 
t.ion  .yl'"'  hr  APSE)/:’  compared  reasonably  veil  with  the  very  limited  experimental  penetration  data;  they  show 
the  depth  oi  penetration:'  to  be  very  dependent  upon  near- surface  moisture  conditions.*  upe rating,  principles 
of  hand- emplaced  miniaturised  ST  and  air-delivered  SIPV.  used  in  special  tests  to  yield  data  lor  aiding 
in  the  evaluation  of  techniques  for  pro  diet  ilk’:  SIP  performance  are  given  in  Appendix  A.  The  test  sites  are 
described  and  supplementary  re  fraction  seismic  data  are  givur.  in  Appendix  L.  The  equations  for  predicting 
penetration  depths  and  deceleration  valuer,  for  air- implanted  SID's  are  presented  in  Appendix  0,  and  the  com- 
puter models  used  to  predict  JID  detection  performance  are  described  in  Appendix  L’ • 
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Factor 

Class 

VS2 
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1 
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4 

1 

1 

5 

1 

1 

b 

1 

2 
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NOTE  ER  DENOTES  ERRATIC  PERFORMANCE  THE  ER  AREAS  ARE  ALONG  THE  ARROYOS. I C.  DRAINAGE 
CHANNELS)  AND  CONTAIN  SLOPING  AND  IRREGULAR  TERRAIN  AND  A GRAVELLY  SANDY  SURFACE 
MATERIAL.  THE  SEISMIC  RESPONSE  IN  THIS  AREA  DEPENDS  TO  A LARGE  EXTENT  UPON  THE 
RELATIVE  LOCATIONS  OT  THE  DROP  HAMMER  AND  THE  SEISMIC  DETECTOR. 
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